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Abstract: We perform a detailed investigation of multiple hard interactions in hadron- 
hadron collisions. We discuss the space-time, spin and color structure of multiple interac- 
tions, classify different contributions according to their power behavior and provide several 
elements required for establishing all-order factorization. This also allows us to analyze 
the structure of Sudakov logarithms in double hard scattering. We show how multiparton 
distributions can be constrained by connecting them with generalized parton distributions 
and by calculating their behavior at large transverse parton momenta. 
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1 Introduction 

When two hadrons collide at high energies, more than one parton in one hadron can have 
a hard interaction with a parton in the other hadron and produce particles with large mass 
or transverse momentum. The effects of such multiparton interactions are suppressed or 
average out in sufficiently inclusive observables, but they have important consequences for 
the details of the hadronic final state. The possible importance of multiparton interactions 
has been realized long ago [1, 2] and phenomenological estimates have been given for 
many final states such as four jets (possibly including b quarks) [3-9], jets associated with 
photons or leptons [10], four leptons produced by the double Drell-Yan process [11-13] or 
from two charmonium states [14-16], as well as a number of channels with electroweak 
gauge bosons [17-26]. Experimental evidence for multiple hard scattering has been found 
in the production of multijets [27-29] and of a photon associated with three jets [30-33]. 
A mini-review of the subject can be found in [34] and an overview of how multiparton 
interactions are modeled in current Monte Carlo event generators is given in [35]. 

At LHC energies, the phase space for having several hard interactions in a proton- 
proton collision is greatly increased compared with previous experiments, and it is expected 
that the effects of multiple interactions will be important in many processes [36-39]. This 
poses a challenge in searches for new physics and at the same time offers the possibility to 
study multiple interactions in much more detail than before. First experimental results on 
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multiple hard scattering at the LHC have already appeared [40] and more can be expected 
in the near future [41]. 

Understanding multiparton interactions is also important for heavy-ion physics, where 
pp or proton-nucleus collisions are used as a baseline for collective effects in nucleus-nucleus 
collisions. Compared with pp collisions, multiparton interactions with nuclei have the 
additional feature that the different scattering partons may come from the same nucleon 
or from different nucleons in the nucleus. Dedicated investigations of multiple interactions 
in pA collisions can be found in [42-46]. 

Phenomenological estimates of multiparton interactions, as well as their implementa- 
tion in event generators, are based on a rather simple and physically intuitive picture, whose 
basic ingredient is the probability to find several partons inside a proton. On the other 
hand, a systematic description of multiparton interactions in QCD has not been achieved 
so far. In the present work, we present a number of steps in this direction. A brief account 
of our main results has been given in [47]. We require all parton-level scatters to have a 
hard scale, so that the concepts of hard-scattering factorization and of parton distributions 
can be applied. Since transverse momenta of final-state particles play a crucial role in the 
characterization of multiple interactions, we fully keep track of this degree of freedom and 
base our discussion on transverse- momentum dependent multiparton distributions. 

In section 2 we give a lowest-order analysis of multiple hard scattering. We find that 
the intuitive picture just mentioned emerges for a subset of all relevant contributions to 
the cross section, but that there are other contributions which may be of comparable size 
and hence call for further investigation. In section 3 we take first steps to extend exist- 
ing factorization theorems for single-hard scattering processes with measured transverse 
momentum [48-51] to the case of multiple hard scattering. While many ingredients for a 
full proof of factorization are still missing (and the possibility that factorization is broken 
cannot be ruled out), we obtain a number of encouraging results that allow us in partic- 
ular to analyze the structure of Sudakov logarithms. Section 4 gives more details about 
the distribution of two quarks or antiquarks in the proton, in particular about the effects 
of spin correlations and the possibility to learn more about multiparton distributions by 
calculating their moments in lattice QCD or by linking them to generalized parton dis- 
tributions. The predictive power of perturbation theory is increased in kinematics where 
all observed transverse momenta (as well as their vectors sums) are large on a pertur- 
bative scale. Complications and simplifications that arise in this regime are discussed in 
section 5, where we will also encounter the conceptual problem of separating single from 
multiple hard-scattering contributions in a systematic and consistent fashion. Section 6 
contains our conclusions. 

2 Lowest order analysis 

2.1 Momentum and position space structure 

In this section we investigate the structure of multiparton interactions in momentum and 
position space, restricting ourselves to graphs with the lowest order in the strong coupling. 
To avoid a clutter of indices we consider scalar partons described by a hermitian field (p, 
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deferring the inclusion of spin and color degrees of freedom to sections 2.2 and 2.3. Our 
derivation of the cross section formula for multiparton interactions uses standard methods. 
For cross sections integrated over transverse momenta in the final state, similar derivations 
can be found in the literature [52, 53]. The extension to cross sections differential in 
transverse momenta is new. For ease of language we refer to the colliding hadrons as 
protons throughout this work, bearing in mind that our results apply without change to 
pp collisions or to any other hadron-hadron collision. 

2.1.1 Definition of multiparton distributions 

We begin by defining the multiparton distributions that appear in the cross section formula 
we will derive shortly. The following definitions need to be completed by a prescription 
to renormalize ultraviolet divergences and by Wilson lines that take into account collinear 
and soft gluons as required to achieve factorization for the cross section. These issues will 
be discussed in section 3. 

The building block from which multiparton distributions can be defined is the n parton 
correlation function 
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where T denotes time-ordering and T anti-time-ordering of the fields. This function de- 
scribes the emission of n partons in a scattering amplitude and in its complex conjugate. 
Throughout this work we assume an unpolarized target: if the target carries spin then an 
average over its polarization is implicit in (2.1) and all subsequent expressions. The parton 
four-momenta in the correlation function are subject to the constraint 
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In (2.1) we have chosen the position of the first field in the matrix element to be ^'^ = 0. 
Taking this position as arbitrary and integrating over it with a factor exp(— i^^Z^) yields 
a delta function for the constraint (2.2). The structure of the cross section will be more 
transparent if we use symmetric variables 
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The constraint (2.2) then turns into 
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(2.4) 
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Figure 1. Assignment of momentum and position arguments in the multiparton correlation func- 
tions and distributions. The dashed line denotes the final-state cut. 



and we can rewrite the correlation function (2.1) as 
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where we have replaced r„ using (2.4). In addition we have used translation invariance 
to shift position arguments in the matrix element by ^n/2. Substituting position variables 
according to 



and Zn = Cn, we obtain 
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The assignment of momentum and position arguments is shown in figure 1. 

We now introduce light-cone coordinates = {v^ it f^)/-v/2 and v = {v^ ^v^') for any 
four- vector f . In a frame where p = we define multiparton distributions 
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This can be written as 
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where we have used the abbreviation 



0{yi, Zi) = (f>{yi - \zi) id^(l){yi + \zi) 
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(2.10) 



for the bihnear parton operators and traded the factors kf for derivatives = •^{d —d)'^ 
acting on the fields. When going from (2.8) to (2.9) we have replaced the time- or anti- 
time-ordered products appearing in (2.7) by usual products, which are understood to be 
normal ordered. To justify this it is crucial that the arguments of all fields in the operators 
(2.10) have a vanishing plus-component. For a generic configuration with all yi and Zi 
different from zero and from each other, all fields in (2.10) have a spacelike separation, so 
that they commute because of causality and can be written in any order. The case where 
fields have a lightlike separation requires special treatment, and different methods for this 
case have been used in the literature for related matrix elements, see [54, 55] and [56]. As 
we shall see in section 5, lightlike field separations in (2.9) also lead to divergences that 
need to be regulated. 

We also introduce distributions that depend partially or entirely on transverse positions 
{yi and Zj) instead of transverse momenta (fcj and r^): 
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and 
F{xi,Zi,yi) 
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In the arguments of (2.11) and (2.12) it is understood that the average transverse position 
of the first two field operators is y^ = 0. The three forms (2.9), (2.11) and (2.12) can be 
used interchangeably, and each of them has advantages in different situations. As we shall 
see, the momentum representation (2.9) naturally appears in Feynman graph calculations, 
the mixed representation (2.11) has a rather simple physical interpretation, and the position 
space representation (2.12) is most convenient for the discussion of Sudakov logarithms. 

The factors of 27r, kf , and 2p+ in (2.8) to (2.12) have been chosen such that the 
coUinear (i.e. transverse-momentum integrated) distribution 
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as well as the distribution 



n— 1 „ 

n J ^'y^ 



F{xi, ki,yi) = F{xi, ki, Vi = 0) 



(2.14) 



admit a probability interpretation. F{xi,yi) is the probability to find n partons with 
plus-momentum fractions Xi and transverse distances yi from parton number n, and 
F{xi,ki,ri = 0) is the probability to find n partons with plus- momentum fractions Xi 
and transverse momenta ki. 

By contrast, F{xi,ki,yi) is not a probability (due to the uncertainty relation one 
cannot simultaneously fix transverse momentum and transverse position) but rather has 
the structure of a Wigner distribution [57] in the transverse variables. Its integral over all 
ki gives the probability to find partons at transverse positions yi, and its integral over all yi 
gives the probability to find partons with transverse momenta ki. A related interpretation 
for generalized parton distributions can be found in [58]. In figure 1 we can identify ki 
as the "average" transverse momenta of the partons and yi as their "average" transverse 
position, where the "average" is taken between the partons to the left and to the right 
of the final-state cut in the figure. In a physical process, this corresponds to an average 
between partons in the scattering amplitude and its complex conjugate. 

The interpretation of multiparton distributions becomes more explicit if one represents 
them in terms of the light-cone wave functions of the target, see [59]. Most conveniently 
derived in the framework of light-cone quantization, this representation is analogous to 
the wave function representation for single-parton densities [60] and generalized parton 
distributions [61, 62]. The distributions in (2.13) and (2.14) can be written in terms of 
squared wave functions in impact parameter or transverse-momentum space, which makes 
their probability interpretation manifest. The wave function representation also offers a 
way to model multiparton distributions in the region of large momentum fractions, where 
one can expect a small number of partonic Fock states to be dominant. We shall not pursue 
this avenue in the present work. 

In later chapters we will also need collinear distributions that depend on the momentum 
transfer variables r^. 
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We will see in the following section that F[xi,yi) or equivalently F[xi,ri) appear in 



multiple-scattering cross sections. This is not the case for the distributions 



F{xi)=F{xi,ri = Q) 



j{ d^k. n 



=1 
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F{xi, ki, yi) , 



(2.16) 



which give the probability to find n partons with momentum fractions Xi and unspeci- 
fied transverse positions or transverse momenta. We note that the integrals over ki in 
(2.13), (2.15) and (2.16) are logarithmically divergent and require appropriate regulariza- 
tion, which will be discussed in sections 5.1.2 and 5.3.2. 
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Figure 2. Graph for the cross section of a collision with n hard scatters at parton level. The 
dashed line denotes the final-state cut. Here and in the following, the lower blob is associated with 
the right-moving proton and the upper blob with the left-moving proton in the collision. 

The definitions in this section are given for right-moving partons, with Xi being plus- 
momentum fractions. Analogous definitions for left-moving partons are obtained by ex- 
changing the plus- and minus-components of all position and momentum vectors. 

2.1.2 Cross section for n hard scatters 

We now evaluate the cross section for a process with n scatters at parton level, as sketched 
in figure 2. We work in a reference frame with p = p = and consider kinematics where 
the squared cm. energy qf of each scatter is large and where each transverse momentum 
IgJ is much smaller than qf' and q^ . Defining 

X, = qr/p- , (2.17) 
we can then approximate 

qf^2qtq-^XiXiS, (2.18) 

where s = {p + p)^ is the squared overall cm. energy. We neglect the target mass through- 
out, so that s ~ 2pp K, 2p^p~ and the flux factor in the cross section is l/(4pp). One can 
trade the momentum fractions Xi and Xi for qf and the rapidities 

Y^ = llog^ (2.19) 
2 g. 

with 

dxi dxi = - diqh dYi , (2.20) 
s 

where we have again used (2.18). We note that for the very high s achieved at the LHC, 
both Xi and Xi are rather small, except if or qf is very large. 
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The cross section for n hard scatters is given by 
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where the combinatorial factor C contains a factor k\ for each set of k identical hard- 
scattering final states."^ The remnant of proton p (p) consists of m (m) spectators with 
momenta px,j (Px j)- denotes the squared matrix element for the ith hard scatter, 
with truncated propagators of the incoming parton lines. Hi includes integration over the 
internal phase space of the final state produced by the hard scatter, with only the four- 
momentum Qi kept fixed. If such a final state is the decay product of a single particle with 
mass M and width T (e.g. a or a Higgs boson) then Hi includes a factor 
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which in the limit of narrow width constrains qi to be on the mass shell. If the final state 
is a stable single particle with mass M, then Hi includes a delta function 



27T6{qf - M2) 



(2.23) 



so that together with the integration element d^gj/(27r)'* in (2.21) one obtains the correct 
one-particle integration measure d^q^/ \2qf {27:)^Y We now rewrite the cross section in 
terms of the correlation functions (2.1). To this end we use 
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^An often used notation for two hard scatters is to write m/2 in the place of 1/C, with m = 1 if the 
hard-scattering final states are identical and m = 2 if they are distinct. 
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Using the analogous relation for the matrix element between X and p and rewriting the 
momentum conservation constraint in (2.21) as 



(n m m \ 

i=l j=l j=l ^ 

(n n n \ 

i=l i=l i=l ^ 



d Ifi d In 
(2^ (2^ 



Eh + Ypx,j-p) (2vr)V^) + , (2.25) 

1=1 7 = 1 ^ ^1=1 1=1 ^ 



we can express the cross section as 
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where in the last step we have switched to the set of symmetric variables (2.3). They 
have the important property that the kinematic constraints on rj and fj do not involve the 
final-state momenta Qi, which will lead to a great simplification below. 
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Hard-scattering approximation. The parton-level scattering processes involve a hard 
scale, which we collectively denote by ~ without assuming a particular hierarchy 
among the individual squared momenta g?. The case where one of them is much larger 
than the others is of particular relevance for the description of the underlying event, but we 
shall not investigate the consequences of such a hierarchy in the present work. In the graph 
of figure 2 it is understood that partons emerging from the shaded blobs have virtualities 
much smaller than . The components of the various four- momenta thus scale like 



and 
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where A denotes the size of the transverse momenta or the scale of non-perturbative 
interactions, whichever is larger. The momentum conservation constraint 5^^^(rj + fj) 
enforces that the components 

AVQ (2.29) 



are small, although by general scaling arguments they could be of order Q. The constraint 
(5^^)(gj — ki — ki) leads to 



kt 



k7 



(2.30) 



up to relative corrections of order /Q"^ . We make these approximations in the correlation 
functions $ and ^> and see that the longitudinal momenta of the partons entering the hard 
scattering are fixed by the final-state kinematics. In the squared hard-scattering matrix 
element Hi{qi,ki,ki,ri,fi) we can neglect all transverse momenta and all components of 
order A'^/Q. With (2.30) this only leaves a dependence on the independent variables qf 
and q~ . Since Hi is invariant under a boost along the z axis, it can then only depend on 
2qfq~ qf. Altogether we then have 
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Inserting this into the cross section (2.26) and using the definition (2.8) of the multiparton 
distributions gives 
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Rewriting d'^qi = p'^p dxi dxi d^qi, we obtain our final result for the cross section in mo- 
mentum representation, 
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for the ith parton-level subprocess and used the approximation (2.18). We have carried out 
the integrations over f j using the constraints 6^'^\ri + ri), so that the distributions for the 
two protons are evaluated at opposite values of their last arguments. Fourier transforming 
these to position space, we have 
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and the distributions are evaluated at equal values of yi . Transforming also the arguments 
ki and fc,-, we have 



da 



nr=i dxi dxi d^qi C 



ai{xiXiS 

i=l 
n-l 

'2 



n 



dh 



n 



i=l 



\J (2vr)2 
F(^Xi, Zi, yj) F(^Xi, Zi, yi) , 



(2.36) 



- 11 - 



where all position arguments in the two distributions coincide. 

The interpretation of the distributions F{xi,ki,yi) discussed in section 2.1.1 extends 
to the cross section formula (2.35). In each individual hard subprocess, two partons with 
average transverse momenta ki and ki produce a final state with transverse momentum g^. 
The ith scatter occurs at an average transverse distance yi from the nth scatter. The hard 
scatters are approximated to be local in transverse space, so that their average distance is 
equal to the average distance between the colliding partons in each proton. We thus find a 
rather intuitive interpretation of the variables in our cross section formula, provided that we 
"average" the transverse momenta and positions between the amplitude and its conjugate. 
Let us however emphasize that we have obtained (2.35) from calculating Feynman graphs 
using standard hard-scattering approximations, without any appeal to classical or semi- 
classical arguments. 

Integrating the cross section over all transverse momenta Qj we obtain a simple result 
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n"=i dxi dxi C 



ai{xiXis) / d^y. 
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in terms of collinear multiparton distributions. This formula has long been known and 
provides the basis of most phenomenological analyses of multiple interactions in the liter- 
ature. It was derived in [52] for scalar partons in a way very similar to the one we have 
employed here. 

2.1.3 Single vs. multiple hard scattering 

The approximations we have made in the previous section give the leading term of an 
expansion in powers of K/Q. Let us investigate how the resulting cross section (2.35) 
scales with Q. As can readily be seen from its definition (2.11), the mass dimension of 
F{xi,ki,yi) is —2 and one has F{xi,ki,yi) A~^. To obtain this power behavior, it is 
essential that the distribution is invariant under a boost along the z axis. For instance, 
a hadronic matrix element that transforms like the plus-component of a vector would 
be proportional to or another large plus component and thus scale like Q times the 
appropriate power of A. Note that the dependence of F{xi,ki,yi) on the large scale Q 
via renormalization group or Sudakov logarithms (see section 3) is neglected at the level 
of power counting. The hard-scattering cross sections have a power behavior ai ~ Q~^, 
and the integrations over transverse momenta count as d^kid'^ki 6^'^\qi — ki — ki) ~ A^. 
Finally, the distances yi in (2.35) are generically of size 1/A so that d'^yi ~ A~^. Putting 
all ingredients together, one finds 



da 



nr=i dxi dxi d'^Qi 



^ (2.38) 

multiple 



A2 Q2n 



for the cross section of n hard scatters. One obtains of course the same result if the power 
counting is done for the representations (2.33) in momentum space or (2.36) in position 
space, using d'^Vi ~ A^ or d'^Zi ~ l/A^. 
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Let us compare this with the cross section for producing the final states with momenta 
Qi in a single hard scattering. With 



j=i j=i 
the factorization formula for this case reads 
da 

single 



E 

i=l 



p- 



(2.39) 



dx dx dq 



a{xxs) j d^k d^k 6^^^ {q-k-k) f{x, k) f{x, k) , (2.40) 



where a is the appropriate hard-scattering cross section and f{x, k) and f{x, k) are trans- 
verse-momentum dependent single-parton densities. The definition of f{x,k) can be ob- 
tained from (2.11) by setting n = 1, which gives a power behavior f{x, k) ~ A~^. We now 
make (2.40) differential in the internal momentum variables of the final state, which we 
choose as 



Ui = xjx = qf/q^, 
and with i = 1, . . . , n — 1. We then have 



Xi/x = Qi /q" 



(2.41) 



da 



nr=i dxi dxi d'^qi 



da 



single nr=i dUidUid'^qi 



d-kd^-k6^^\q-k--k)l^ 



7.n~l 



(2.42) 



The differential hard-scattering cross section on the r.h.s. behaves as Q , so that we have 



da 



nr=i dxi d^qi 



single 



A2Q 



2n 



(2.43) 



We obtain the important result that if one leaves the cross section differential in the trans- 
verse momenta Qj, the contributions from single and from multiple hard scattering have 
the same power behavior in the large scale Q, so that multiple hard scattering is not power 
suppressed. It is easy to see that the power behavior in (2.38) and (2.43) holds for any 
combination of single and multiple hard scatters, e.g. when producing the final states with 
momenta qi and q2 in a single hard scatter and each final state with momentum ^3, 54, 
etc. in a hard scatter of its own. 

Let us now see what happens if we integrate over the q^. Li the multiple-scattering 
mechanism, each transverse momentum q^ is the sum ki + ki of two parton momenta and 
thus limited to be of size A, so that the phase space volume is HILi d'^^i ~ A^". With a 
single hard scattering, however, the individual momenta q^ can be as large as the hard scale 
Q, and only their sum q is limited to be of order A by the constraint q = k + k in (2.42). 
The phase space volume in this case is therefore U7=i d^Qi = d^q UZi d^q^ ~ A^Q^"""^, 
and we have 



da 



nr=i dxi dxi 



A' 



2n-2 



do- 



multiple 



Q 



2n 



nr=i dxi dxi 



single 



1 



(2.44) 
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for the cross sections integrated over all transverse momenta. Multiple-scattering contribu- 
tions are now suppressed by at least one power of A^/Q^ and are hence power corrections 
to the contribution from a single hard scattering, as has been known for a long time [63]. 
This is indeed necessary for the validity of the familiar collinear factorization theorems, 
which only take into account single hard scatters. 

Power counting in the hard scale Q provides an essential criterion for determining which 
contributions to the cross section are important. There are, however, other important 
factors to keep in mind. We already mentioned Sudakov logarithms in q^/Q"^, which 
appear in the cross section differential in and are different for single and multiple hard 
scattering. They will be discussed in section 3.4. Another aspect in which single and 
multiple scattering contributions differ is the dependence on the momentum fractions Xi 
and Xj, which can be rather small as we remarked after (2.20). We will return to this point 
in section 2.4. 

2.1.4 Impact parameter representation 

The cross section in (2.36) involves distributions F{xi, Zi,yi) that depend on the transverse 
positions of the scattering partons but still refer to proton states with definite (zero) trans- 
verse momenta. In this section we give a formulation completely in transverse position 
space, closely following the construction of impact-parameter dependent parton distribu- 
tions in [62, 64, 65]. 

To begin with, we define a non-forward correlation function ^{l^,l^;p,p') exactly as 
in (2.1) but with a state {p'\ having a different momentum than the state \p). Using the 
same arguments as in section 2.1.1 we can derive a representation of the form (2.7) for 
<I>(/j, /^;p,p'), with {p\ replaced by {p'\. The constraints on the parton momenta read 

n n n 

p-J2h=P - Y.''^=P-P- (2-45) 

i=l 1=1 1=1 

in this case. In the same manner we define multiparton distributions F(xi,ki,ri;p,p'), 
F{xi,ki,yi;p,p') and F{xi, Zi,yi;p,p') as in (2.8) to (2.12), but taken between states 
{p^iP'l Ip^jP)- Note that we take the same plus-momentum in the bra and ket state, 
even if their transverse momenta are different. 

We now consider a transverse boost, i.e. a Lorentz transformation that changes the 
transverse components of a four-vector v as 

v^v-v"^ ?-^LE_ (2.46) 
2p+ 

and leaves plus-components unchanged. Invariance under this transformation implies 

F{xi,ki,ri]p,p') = F[xi,ki - XiP,ri;-\^,\/^) (2.47) 

with 

P = \{p + p'), A=p'-p. (2.48) 
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In impact parameter space we then have 



F{xi,Zi,yi;p,p') 



n/d2fc.e-'^^ n/ 



'-1=1 



(2vr)2 



F{xi,ki,ri;p,p') 



(2vr) 



X F(xi,fci-XiP,ri;-iA,iA) 

= e^PEr=i-.-.F(x„^,?/,;-iA,iA). 

We now introduce proton states with definite impact parameter: 

d^p 



One readily obtains their normahzation 



(2vr)^ 



-ibp 



P ,P) 



{p'+, b'\p+, b) = 47rp+6{p'+ - p+) 5(2) {b' - b) 



(2.49) 
(2.50) 
(2.51) 



from the usual relativistic normalization = (2vr)'^ 2p'^ Sip''^ — P^) 5^'^\p' —p) 

of momentum eigenstates (recall that at fixed p one has dp^/p^ = dp'^ /p^ in the invariant 
integration element). For later use we also give the projector on one-particle states, 

dp+d'^b 



1 



4:TTp^ 



\p , , b\ 



(2.52) 



which is readily checked by taking the matrix element between the one-particle states in 
(2.50) and using (2.51). We finally define the center of momentum of m particles with 
plus-momenta p^ and transverse positions b^ as 

m , m 

i=l ' i=l 

By virtue of Lorentz invariance, this is a conserved quantity. Note the analogy between 
(2.46) and non-relativistic boosts if d is a momentum and if one replaces plus-momenta 
by masses. The center of momentum is thus the analog of the center of mass in the 
non-relativistic case, which is of course conserved. 

Let us consider the matrix element of the same operator as in (2.12), but taken between 
impact parameter instead of transverse-momentum eigenstates. We have 

■n— 1 „ -1 n—1 

n 2p+ dyr {p+^-h-\d\0{Q,Zn)J{0{yuZi)\p+,-b+\d) 



i=l 



dz 



i=l 



i=l 



I 



d^p' d^p i(p'-p)b-i(p'+p)d/2 r?( /N 

(27r)2 (2vr)2 r [x^, z„y„p,p ) 

d^A d^P ib^-iPd+iPY.'l^^xiZi pf t,.--lA iAl 
(27r)2 (27r)2 ^i, Ui-, 2^) 

5i^)[d-Y^^x,z^ y l_^e-''^F(x„^„y,;-iA,iA). 



(2.54) 
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The delta function in the last line reflects the conservation of the center of momentum, 
which equals 

n-1 



i=l 
n-1 



(2.55) 



i=l 



for the bra and ket states in the matrix element, respectively. Here Xs = 1 — XlILi -^i 
and Zg is the center of momentum of the spectator partons. We define impact-parameter 
dependent multiparton distributions by 



/d A. 
__ (,-^b^F{xi, z,,yi; -i A, i A) . 



(2.56) 



If we set Zj = then the matrix element in (2.54) is taken at d = and hence becomes 
diagonal. We can interpret F(xj,Zj = 0,yi;b) as the probability to find n partons with 
plus-momentum fractions Xi in a target that is localized in impact parameter space, with 
parton number n at a transverse distance b from the center of the target and partons 1 to 
n — 1 at relative transverse distances yi from parton n. 
Inverting (2.56) and setting A = we get 



F{xi,Zi,yi) = j d^bF{xi,Zi,yi;b) 



(2.57) 



and can therefore represent the multiple-scattering cross section (2.36) as 



da 



n"=i dxi dxi d^qi 



1 

C 



n 

i=l 
n-1 

n 



(JiiXiXiS 



d^yi 



n 



(2vr)2 



d^bd^b Fixi, z„ yf, b) F{x^,z„yi; b) . (2.58) 



Integration over leads to Zi = as in (2.37). The resulting cross section formula was 
already derived in [66], and it has a very intuitive geometric interpretation shown in figure 3. 
As already noted after (2.36), the approximations we have made for the hard-scattering 
subprocesses imply that each pair of colliding partons in the hadrons p and p must be at the 
same position in impact parameter space. The relative distances y^ between the partons 
are hence the same in both hadrons, but the distance of the partons from the center of 
their parent hadron is in general different in p and p. The relative transverse distance b — b 
between the hadrons is integrated over in the cross section. 

Our result (2.58) shows that the representation of the cross section in terms of impact- 
parameter dependent distributions remains simple even if the transverse momenta q^ are 
kept fixed. In the geometric interpretation just described, we then have to replace "dis- 
tances" by "average distances", with the average taken between the amplitude and its 
conjugate. What is lost in this case is a probability interpretation of the multiparton dis- 
tributions. The two fields associated with a parton in the target are now taken at a relative 
transverse distance Zi, whose typical size is \zi\ ~ 
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Figure 3. Visualization of the cross section formula (2.58) for n = 2 when Qi and are integrated 
over. Each hard scatter produces a heavy gauge boson in this example. 



2.1.5 Reduction to single-parton distributions 

In order to build a phenomenology of multiple interactions, one needs a simple ansatz for 
multiparton distributions as a starting point. It is natural to approximate those distribu- 
tions that have a probability interpretation by the product of single-parton densities. In 
this section we show how one can formally implement this approximation and generalize it 
to the distributions F{xi, Zi,yi) or F{xi,ki,yi), which do not represent probabilities. 

To this end we insert complete sets of intermediate hadron states in the operator 
product appearing in the multiparton distributions: 

71—1 r rt— 1 

O{0,zn)ll O{y„Zi) = O{0,Zn) Yl Y,\Xi){X^\0{yi,Zi) 

1=1 L i=l Xi 

rn-1 



O{0,Zn)\Xn^l) ll{Xi\0{y,,Zi)\Xi^i] 

X„-i,...,Xi 



1=2 



[xAO{yi,zi). 



(2.59) 



Note that the two parton fields in each operator 0{yi,Zi) are associated with the same 
plus- momentum fraction Xi in the multiparton distributions. 

The approximation that gives a product of single-parton distributions is to assume 
that among all intermediate states \Xi) the dominant ones are single-proton states. This 
reduces the complete sets of intermediate states to the projection operators (2.52), and one 
obtains 



6^''^(d-Y,x^Zij F{x,,z,,yi;b)^ 11 
dpf£h 



1=1 

71-1 

n 

1=1 

71-1 



1 = 1 



dz- 



n-l 



n dy, 



1=1 



p+,-fo-id|o(o,z„)|p+_i,b„_i; 



II {pf ,b^\0 [yi ,Zi)\pt_^,h^^^) {p+,h^\0{yi,z{)\p+,-h + \d] 



(2.60) 



1=2 



Translation invariance and the definition (2.50) of impact-parameter states imply 
{ptA I 0{y,,z.,)\pU,h,_,) = ^Tipt-pt-d {p+,b,-y, \O{0, z,)\pt_„b,_, - y.,_,) (2.61) 



and hence 

n 

(5^^^ (rf - ^ XizA F{xi, Zi, yi] b) 



i=i 



n / d\ 



i=l 

n-1 



dz, 



2tt 



(p+,-6- id|O(0,z„)|p+,6„„i) 



i=2 



27r 
dz^ 



27r 

Using (2.54) for n = 1, we have 
dz 



n / ^eW+(p+,6i-yi|O(0,z,)|p+,&.-i-y^-i 



(2.62) 



27r 



e"^ P (p+,-6-id|O(0,z)|p+,-6+id) = 5(2)(rf_2;z)/(a;,z,fo), (2.63) 



where f{x,z;b) can be written as 

/(x,^;fe) = J ^e-''^/(x,z;A) 

with 



/(a;,^; A) 



2tt 



^ixz p+ (p+ iA|O(0,z)|p^ 



' 2' 



A reader famihar with generahzed parton distributions wih recognize that 

dh 



f{x,k;A) = I A£_e— ^/(x,2;A) 



(2.64) 



(2.65) 



(2.66) 



is a transverse-momentum dependent generahzed parton distribution at zero skewness. We 
wih shortly need the cohinear distributions 



f{x;b) = f{x,z = 0-b) 



f{x-A) = f{x,z = 0;A) 



(2.67) 



as weh. Introduced long ago in [64, 65], the impact parameter density f{x;b) gives the 
probability to find a parton with momentum fraction j; at a transverse distance b from the 
center of the proton. 

The delta function on the r.h.s. of (2.63) implies that 

n-l 

bn-1 = - ^ XiZi + ^XnZn bi 1 = bi + XiZi for 1 < f < n - 1 (2.68) 

1=1 

in (2.62), so that we obtain the desired approximation 

F{xi, Zi, yi] b) f[xn, Z^, b + ^{xiZi + . . . + Xn-lZn^l)) 
n-l 

Y[ f{xi, Zi] b + y,i + \{xizi + . . . + Xi_iZi_i) - i(xj+iZi+i - . . . x„z„)) 

X /(xi,Zi;fo + 2/i - \{X2Z2 + . ■■XnZn)) (2.69) 



1=2 
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x'2 



Vi 



b + Vi 



Figure 4. Illustration of the approximate relation (2.70) for n = 2. 



of a multiparton distribution. Setting Zj = and integrating over 6, we obtain in particular 
the collinear multiparton distribution F{xi,yi) in terms of impact-parameter dependent 
single-parton densities, 



/n—l 
(fh f{xn;b) Yl f{xi;b + yi) 

i=l 



(2.70) 



This relation is illustrated in figure 4, which uses the representation of parton distributions 
as squared light-cone wave functions we mentioned briefly before (2.15). 

Let us now insert (2.69) into the cross section (2.58). For measured transverse momenta 
Qj, the different single-parton distributions are entangled by their Zi dependence. By 
contrast, the integrated cross section simplifies to 



da 



1 



nr=i c 



XiXiSj 



-i=l 
n-1 



d^bd^bf{xn;b)f{xn;b) 



n / d^Vi f{xub + yi) f{xi;b + yi) 
i=i '' 

^ j ^"^P Yl^iixiXis) j d'^yi f{xi;yi - p) f{xi;yi) 



(2.71) 



where the only integration variable linking the different factors is the relative distance 
p = b — b, and where we have renamed the integration variable b to y^ in the second step. 
In different forms, this relation (or more precisely its analog for quarks and gluons instead 
of scalar partons) has long been used as a starting point of phenomenological studies, see 
e.g. [67-71] and [8, 72]. ^ 

As observed in [59] for the case of collinear distributions, the reduction of multiparton 
to single-parton distributions also takes a simple form in the transverse-momentum repre- 
sentation. This remains true if one keeps the transverse parton momenta unintegrated. To 
see this, we integrate (2.69) over b and Fourier transform w.r.t. yi and Zi as specified by 
(2.11) and (2.12). Changing integration variables from b and y^ to the impact parameter 



^ We note that in [8, 72] the impact parameter arguments of / are p — yt and yi instead of — p and yi 
(if we translate to our notation). This is equivalent in the spin independent sector, where the single-parton 
distributions are independent of the direction of the impact parameter. 
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arguments of the distributions on the r.h.s. of (2.69), we obtain 



F(Xi ^ki^Ti) ~ / (yXji^ kji 2 (r\ -\- . . . -\- Tn—\)i fj-. 



r-n-l 



Y[ f{xi,ki - ^Xi{ri + ... + ri_i) + ^Xi{ri+i + . . .r„);rj) 



i=2 



f{xi,ki + ^xi{r2 + ... + Vn); ri) , 



(2.72) 



n-l 



where we recall that Vr 



^ Vi. Integrated over the momenta ki this simply reads 



i=l 



so that the cross section (2.71) becomes 

■n-l , ,2 



(2.73) 



da 



1 



nr=i dxi dxi c 



n 

1=1 



(2vr)2 



JJ ai{xiXis) f{xi; n) /(x^; -r^ 



i=l 



(2.74) 



The arguments in (2.73) can easily be anticipated from figure 1. 

We emphasize that the relations (2.69) to (2.73) have been obtained by restricting 
a sum over all intermediate states to a single proton. We do not have a motivation for 
this restriction other than observing that it results in neglecting correlations between dif- 
ferent partons in the proton. It seems plausible to assume that this is a reasonable first 
approximation, at least in a certain region of variables, but one should not expect it to 
be very precise. Possible deviations from this approximation and their phenomenological 
consequences have recently been discussed in [8, 73-77]. 

2.2 Parton spin 

Let us now see how the scattering formulae (2.33) to (2.37) are modified in QCD, where 
partons have nonzero spin. In (2.21) to (2.31) the squared amplitude Hi of the ith hard 
scattering and the hadronic matrix elements of parton field operators acquire spinor indices 
in the case of quarks and Lorentz indices in the case of gluons. These indices can be treated 
as in the case of a single hard scattering. For the time being we still omit color degrees of 
freedom, which will be discussed in section 2.3. 

2.2.1 Quarks 

The correlation function for n quarks entering the hard scattering is 

-n-l 



$ai/3i...a„/3„(^i,'^i) 



n 

4 = 1 



d Zj 



n 

i=l 



d Vi 
(2^)4 



-Wiri 



{p\t qp^{yi - \zi) ■■■qpS-\zn)\ T\^qaShn) ■ ■ ■ qa^{yi + \zi) 



(2.75) 
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When this is integrated over the parton minus-momenta, the anti-time and time ordering 
can be omitted and one can reorder the fields as 



n-l 

n 



(2.76) 



by an even permutation. For antiquarks entering the scattering, one has an operator 
product QaAVi - Qp.iVi + instead of qp.ivi - \zi) Qa^iVi + ^Zi). 

Consider the case of a quark entering the hard scattering. We wish to rearrange the 
spinor indices in the product ^apHi^j^a-, where for brevity we write a,/3 instead of aj,/3i 
and leave out all other indices on which Hi and $ depend. The rearrangement is achieved 
by the Fierz transform 



Hi,pa = ^S(Satr{^Hi) + l(75)^^tr(i75/7i) + i(7^)^„ tr(i7^ifi) 
+ i(7^75)/3atr(i757^ifi) + liia'"'j5)pati{liau^,j5Hi) . 



(2.77) 



The Dirac matrices with open indices on the r.h.s. multiply fields q[s{yi — \zi) qa{yi + ^ 



Zi) 

in the correlation function <&q,/3 for the right-moving proton. The dominant terms in the 
cross section are those where that matrix is F = ^7"*", ^7"^ 75 or ^ia~^^j5 with j = 1,2, 
because F^q,<I>q^ is then proportional to the large momentum component p"*" ~ Q by virtue 
of Lorentz invariance. The traces over the hard scattering matrix Hi on the r.h.s. of (2.77) 
have both large plus and minus components since Hi depends on the boson momentum qi. 
One thus has 

^apHi^pa = tr(i7+^) tr(i7~^^*) + tr(l7+75$) tr(l757"i^i) 

-|- tr(^icj-'''"75<I>) tr(^i(7-'~ 75/^4) -|- {power suppressed terms} , (2.78) 

where a sum over the transverse index j = 1, 2 is understood. 

When defining distributions for scalar partons in (2.8), we included a factor for 
each parton i = 1, . . . , n. For quarks we do not do this, but instead include this factor 
in the definition of the parton- level cross section dj from the squared matrix element Hi. 



Writing k- ^ for the collinear approximation of fcj (i.e. k^ 
recognize in 



and ki 



kttr{l^-H,) = ltr{ 



"i,cHi^ — 2 ^ ^ '^t 



{ki,c) Hi Usiki^cj 



0) we 



(2.79) 



the spin averaged squared amplitude for an incoming on-shell quark. The corresponding 
terms with 7~75 and icr"-' 75 = 75 7''7~ are respectively associated with scattering on a 
longitudinally and transversely polarized quark. 

Integrating (2.75) over the minus components of the parton momenta, one obtains 
multi-parton distributions as in (2.9) with the scalar field operators (2.10) replaced by 
quark bilinears 



Oa{yi,Zi) = qivi 



^Zi)raq{yi + \zi\ 



(2.80) 
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Figure 5. The numbering of fields in two-parton distributions specified in (2.85). The color indices 
j,j',k and k' will be discussed in section 2.3.1. 



where a = q, Aq, Sq labels the polarization and 



2 ' 



j7 75 



1 ■ i+ 



(2.81) 



We recognize the operators that appear in the definition of single-parton densities for un- 
polarized, longitudinally polarized and transversely polarized quarks, see e.g. [78, 79]. For 
antiquarks entering the hard scattering one proceeds in an analogous way. The correspond- 
ing operators are 



OaiVi^Zi) = -qivi + ^Zi)Taq{yi 



=y+=0 



(2.82) 



with 



q 5 



(2.83) 



The overall minus sign in (2.82) reflects a change in the order of fleld operators from 
QaiUi -\zi)qi3{yi + \zi) to qj3{yi + \zi) qa{yi - \zi) , cf. our remark after (2.76). In the case 
of O^q a further minus sign is included in T^q, so that the operator corresponds to the 
diff'erence of antiquarks with positive and negative helicity. 

From now on we concentrate on two-parton distributions. The formalism can be ex- 
tended without conceptual difficulties to higher multiple interactions, but the resulting 
expressions become rather unwieldy. As one encounters nontrivial features already for 
double hard scattering, it is natural to elaborate this case first. To simplify the discussion, 
we introduce a compact notation 



n 

i=l 



dz- (Pi 
(27r)3 



X 2p 



dy {p\Lp{y -\zi)if{-\z2)^p{\z2)<~p{y + \zi)\p) 



=y+=0 



(2.84) 
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for the Fourier transformed matrix element of a product of field operators ip. Their indices 
are assigned according to 



1 




y + 






2 










3 










4 




y - 


\zi 





amentum fraction xi in conjugate amplitude (2.85) 

as shown in figure 5. Throughout this paper we consider unpolarized incident hadrons, so 
that an average over the proton spin is understood in (2.84). A two-quark distribution is 
then given by 

Fai,a2{xi,ki,y) = {{{q3'f^a2Q2){q4'f^aiqi))) , (2.86) 

and if the parton with momentum fraction X2 is an antiquark one has instead 

Fai,a2{xi,ki,y) = {{{q2'i^a2Q3) {qiTaiqi))) ■ (2.87) 

In straightforward extension of the case of single-parton distributions [79], the matrix 
elements defining distributions for quarks and antiquarks are thus connected as 

Faj,a2ixi,X2,ki,k2,y) = -^ai ,a2 (^^l , -3^2 , ^1 , "^2, ?/ ) (2.88) 

with sign factors cjg = asg = +1 and cja^ = —1- Definitions and relations analogous to 
(2.86), (2.87) and (2.88) hold for the case where the parton with momentum fraction xi is 
an antiquark. 

The previous arguments can be repeated for the partons in the left-moving proton, 
with the roles of plus and minus components interchanged. We define the hard-scattering 
cross section for a right-moving quark and a left-moving antiquark as 

a^,aa = ^ [Pa{k,)] [P-a{h)\ H^^p^^-^ (2.89) 

with spin projectors 



Piaik) = PL-ikc) = Hhl' (2.90) 



constructed from the collinear momenta ki^c introduced before (2.79), i.e. kf^ = kf for 
right-moving partons and k'^^ = k~ for left-moving ones, with all other components equal 
to zero. The spin projectors match the Fierz decomposition (2.78) and the operators in 
(2.80) and (2.82), and they can be expressed in terms of quark or antiquark spinors as in 
(2.79). It is understood that for each label 6q or 6q the cross section (2.89) depends on a 
transverse Lorentz index, which has not been explicitly displayed. In most reactions the 
partonic subprocess involves only chirality conserving interactions. Since incoming quarks 
and antiquarks are approximated as massless in the hard scattering, only the combinations 



-23- 




Figure 6. Graphs for double hard scattering initiated by different combinations of quarks and 
antiquarks in the amphtude and the conjugate amphtude. 



^q,qi o'Aq,Aq, ^q,Aq^ ^^q,q ^^'^ ^Sq,5q then nonzero. For parity conserving processes 
such as the production of a virtual photon, one is left with only dq^q, <5"Ag,A(7 and asq.Sq. 
Hard-scattering cross sections 

^i^aa foi" right-moving antiquarks and left-moving quarks are 
defined as in (2.89) with an appropriate change of spinor indices. 

We now have everything at hand to write down the expression for the double-scattering 
graphs of figure 6a and b. For a single quark flavor, one has 



da 



Wi=i dxi dxi d^q. 



11 ,<i2=9i^9i'5<? «=1 

Hi ,B,2=q,'^q,Sq 



I 



X / d^y 



fig. 6a,b 

0"l,aiai (xiXis) a2,a2a2 {x2X2s) Fa^^ai i^i, ki, y) Fa^^ai i^i, ki,y) 



,aiaiixiXis) (5"2,a2a2 {x2X2s) Fai ,a2 {xi, ki, y^ Fai,a2 {Xi, ki, y) 



(2.91) 



where S* = 2 if the final states of the two hard scatters are identical and S = 1 otherwise. It 
is straightforward to Fourier transform the previous expressions either from the interparton 
distance y to the relative transverse momentum r, or from average transverse momenta 
ki,ki to transverse positions Zi, as we did in (2.9), (2.12) and (2.33), (2.36) for scalar 
partons. 
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Notice that (2.91) involves a polarization dependence in the multiparton distributions 
and hard-scattering cross section. This is because, even for unpolarized hadron beams, 
the polarization of the two partons with momentum fractions xi and X2 can be correlated 
among themselves. We will discuss this in more detail in section 4.1.1. 

The two-quark and quark-antiquark distributions considered so far have the form 
((0102)), where the Oi are bilinear operators from (2.80) or (2.82). As we discussed 
after (2.14), these distributions can be interpreted as probabilities or pseudo-probabilities 
in the sense of Wigner distributions for two partons in the proton that carry momentum 
fractions xi and X2, respectively. 

There are further double-scattering graphs that contribute to the cross section and in- 
volve distributions which represent interference terms rather than probabilities. In figure 6c 
we show the case where the parton with momentum fraction xi is a quark in the scattering 
amplitude and an antiquark in the conjugate scattering amplitude. Such interference terms 
in fermion number have no equivalent in single hard-scattering processes, where they are 
forbidden by fermion number conservation. For their description we introduce interference 
distributions 

Iai,a2{Xi,ki,y) = (((^2ra2g4) (gsTaigi))) , 

hi,a2{xi,ki,y) = (((^4ra2g2) (giTaigs))) • (2-92) 

In the absence of a probability interpretation, the choice of quark vs. antiquark labels in 
the Dirac matrices is pure convention. We assign labels such that a indicates a quark and 
a an antiquark in the amplitude, i.e. for the parton indices 1 and 2 in figure 5. The graph 
in figure 6c contributes to the cross section as 



da 



= ^ Yl -f^l,ai/3iQi/3i(^l'^l) [PaAh)]a,f}2 [^'aa (^2)] /^laj 

fig. 6c ai,a2=q,A.q,Sq 
5-1 ,a2=q,^q,Sq 



Yli=i dxi dxi d'^q. 

X -f^2,a2/32a2/32(^2,^2) [Pai{h)] p^a, [^a2 (^2)] ^2^1 

X '[ljd'^kid'^ki6'''^\qi-ki-ki) j d^y Ia^fl^{xi,ki,y) Ia^^a2{xi,ki,y) ■ (2-93) 

i=l 

We see that the contraction of Dirac indices ties together the two hard-scattering kernels 
and the spin projectors P, so that one cannot define separate partonic cross sections ai 
and (72. The power behavior of this contribution is the same as in (2.91). 

Taking different quark flavors into account, we obtain further interference terms. The 
contributions in figure 7a and b involve the interference of different quark flavors, and those 
in figure 7c and d the combined interference in fermion number and flavor. The relevant 
matrix elements are easily written down, reading e.g. (( {uz^a2 d2) {d/i^ax ^i) )) for the lower 
part of figure 7a. 

Which interference distributions are of appreciable size is interesting from the point of 
view of nucleon structure and important for phenomenology. One may for instance imagine 
that diquark-like correlations in the nucleon wave function play an important role in this 
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Figure 7. Graphs for double hard scattering with interference in quark flavor. For simphcity, the 
blobs indicating the hadronic matrix elements are not shown. 

context. In section 2.5 we will argue that for small values of Xj both fermion number and 
quark flavor interference distributions should become relatively small. 

2.2.2 Gluons 

If gluons enter a hard-scattering subprocess, special attention needs to be paid to their 
polarization. In covariant gauges such as Feynman gauge, an unlimited number of right- 
moving gluons with polarization in the plus direction can be attached to a hard scattering 
graph without any power suppression. The effect of these gluons is resummed into Wilson 
lines, which we will discuss in detail in section 3.2.1. Alternatively, one may work in light- 
cone gauge = 0, where the corresponding gluon polarization is absent. One then has 
to be careful about subtle effects from Wilson lines at infinity, see our remark at the end 
of section 3.2.1. 

Once the right-moving gluons with plus polarization (and the left-moving gluons with 
minus polarization) are taken into account, the leading contribution to the cross section 
comes from gluons with transverse polarization, corresponding to field operators with 
J = 1, 2. It is for these gluons that one introduces parton distributions similar to those of 
quarks. To decompose the product of two gluon potentials with transverse polarization, 
we use the relations 

gii'gfefc' _ gjkgj'k' _ ^jk' gj'k ^ 

^jj'M' = 1 (jifc^i'fc' + s^k'^j'k _ ^jj'^kk'-^ ^ (2.94) 



-26- 



where the indices j,j',k,k' = 1,2 are restricted to be transverse and where e^^' is the 
two-dimensional antisymmetric tensor with e^^ = 1. The tensor t^^''^^' is symmetric and 
traceless in each of the index pairs (jj') and {kk'). As an analog of the decomposition 
(2.77) for fermions, we can thus write 



H: 



for the squared hard-scattering matrix element, where in the last step we have used the 



relation r^fM'^n'M' 



The tensors depending on in (2.95) are to be contracted with a product A^' of 
gluon potentials in the multigluon correlation function 



n 



(2vr) 



n-1 

n 

i=l 



X {p\T 



n-l 



n-l 



{-hn) n (y^ - T ^'Hhn) n (y^ + ^0 



i=l 



i=l 



\p) ■ (2.96) 



In analogy to the definition (2.8) for scalar partons (and in contrast to the one for quarks) 
we include a factor kf for each gluon i when defining multi- gluon distributions F from <I>. 
One then obtains 



n 

1=1 



1 



27r 

n-l 



(2vr) 



n-l 



n 2p+/ '^^^^ 



i=l 



(p|a„(0,z„) J]a,(2/„2i)|p>, 



j=i 



where 



0,(?/„ z,) = n^^'G+^^'(y. - \z,) G+\y, + i^,) 
with polarization labels a = g, Ag, 5g and 



a 



(2.97) 



(2.98) 



(2.99) 



The operators Og and 0/\g appear in the usual densities for unpolarized and longitudinally 
polarized gluons. By contrast, describes the interference of two gluons whose helicities 
differ by two units, or equivalently the difference between linear gluon polarization in 
two orthogonal directions. Such distributions have previously been discussed in different 
contexts, see [80, 81] and [82-87]. 

In going from (2.96) to (2.97) we have traded gluon potentials for field strengths 
using the relation G^^ = d'^ A^ valid in the light-cone gauge A'^ = 0. Under the Fourier 
transform this turns kf A^' A^ into {kf)~^G^^'G^^ and explains the factor 1/ixip'^) for 
each parton in (2.97). It is plausible that gluon field strengths rather than potentials should 



-27- 



appear in the definition of gluon distributions, since G^^ lias a simple behavior under gauge 
transformations and can be used to construct gauge invariant operators. How a definition 
with G^^ emerges in Feynman gauge is rather involved and has been shown explicitly for 
the case of a single hard scattering in [88]. 

Parton-level cross sections for gluons are defined as in (2.89) with spin projectors 

pkk'^i^kk'^ pkk^ ^ _i-^kk' ^ ^pu^^kk' ^^ii',kk' (2.100) 

following from (2.95). In Pg one readily recognizes the average over the two transverse 
gluon polarization. The expressions (2.98) to (2.100) are for right-moving gluons. For 
left-moving gluons, one has to change -|- into — coordinates in (2.98) and reverse the sign 
of the e tensor in HAg and P/\g . This is because in a covariant decomposition of the matrix 
elements the two-dimensional e tensor arises from the four-dimensional one as e^^ = e"' 
and thus changes sign when -|- and — coordinates are interchanged. 

Using our shorthand notation (2.84) we can write two-gluon distributions as 

Fa,,a,{x.,k,,y) = {x,p+)-\x2P+)-' {{{u','^' Gf"' Gf") {Ui^;GfGt'))) • (2.101) 

Of course, there are also multiparton distributions involving both quarks and gluons. When 
discussing the mixing of two-quark and two-gluon distributions in section 5.1.3 we shall 
need quark-gluon distributions of the type 

Fa,,aA^^,h,y) = (xip+)"l (((gsra.^s) (n^{G+^' G+^) )) (2.102) 

with fli = g, Ag and 02 = q, Aq. 
2.3 Color 

In contrast to single-parton densities, where two parton fields are always coupled to a color 
singlet, multiparton distributions have a nontrivial color structure. We limit ourselves 
to two-parton distributions here, i.e. to correlation functions with four parton fields. In 
this section we give general decompositions of their color structure. Dynamical aspects 
where color plays an essential role will be encountered throughout section 3, as well as in 
sections 5.1.3 and 5.2.2. 

2.3.1 Quarks 

For two-quark distributions we write 



Fjj',kk' = {{{q2,,k'^a2q2,k) fej-'Taigij))) = 



ON 

'F6,rSkk' + ^^^^'Ft^^,tt,, 



(2.103) 



where j,j' and k, k' are color indices and A'^ is the number of colors. The indices 1 and 2 on 
the quark fields are shorthand for the position space arguments associated with momentum 
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fractions xi and X2, as given in (2.86). For ease of writing we do not display the polarization 
labels oi, 02 of F when not necessary. The functions and ^F can be projected out as 



6j'j 6k'k Fjj',kk' 

2N 



^/iV2 



1 



j.a j.a TP 
^j'j ^k'k ^jj',kk' 



((feTaa 92) (g4rai ^l))) , 

2N 

^^^^ {{{q3Vq2) {q^ra^t^,))) . 



(2.104) 



We see that for = 3 the quark lines carrying the same longitudinal momentum are 
coupled to color singlets and color octets in ^F and ^F, respectively.^ Obviously, only ^F 
admits an interpretation as the joint density of quarks with momentum fractions xi and 
X2, summed over their respective colors. The prefactor of ^F in (2.103) has been chosen 
such that it also appears in the projection (2.104). For this choice the color singlet and 
color octet distributions enter with equal weight 



{^F'^F + ^F^F)/N 



(2.105) 



in the cross section of processes where hard scatters produce color-singlet systems. In this 
sense, the size of ^F relative to ^F directly indicates its relevance to phenomenology. 

For parameterizing the color structure of Fjji^/^y one can alternatively use ^F and the 
matrix element 



6jik Sk'j F, 



j ^jj',kk' 



{{{q3,j^a2q2,k) iq4,k i ai qi, 



ViV2 - 1 

N 



'F + 



1 

N 



'F. 



(2.106) 



in which quark lines carrying different longitudinal momentum couple to color singlets. We 
note that this combination becomes equal to ^F in the limit of large N. It can be rewritten 
in terms of matrix elements 



= {{{qiTa2q2) (qs'^aiqi) 



(2.107) 



that involve bilinear quark operators with no uncontracted color or spinor indices. This is 
achieved by a Fierz transform of Fqj Fq^ w.r.t. the spinor indices of q-sj and qij, followed 
by a Fierz transform w.r.t. the other two indices. Writing 



Cai,a2 = iq3,j^a2q2,k) (^4,^ Faj §1 j) , 

one has 



ai ,a2 



{qA,k^a2q2,k) {q3,j'^a^ qi,j) (2.108) 



^Oq^q + OAq,Aq\ 
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^q,Aq ~ ^Ag.g 


\ 0%M ) 
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(2.109) 



For convenience we use the notation F for general A*'. 
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and 



q,oq ' oq,q 



q,Sq 



Sq,q 



^Aq,Sq ^ ^5q,Aq 



5q 



SqAqJ 



\ 













jk 





le 



\ 









fnk 

^Aq,5q ^ ^5q,Aq 



Sq,q 
(■nk 



k 

Aq,Sq 



ot 



(2.110) 



AqJ 



where r is defined in (2.94) and the transverse indices k, k' of the tensor operators on the 
r.h.s. are summed over as appropriate. The global minus sign in both equations comes from 
the reordering of fermion fields between O and O. The inverse transformation goes with 
the same matrices. Of course, the distributions do not have a probability interpretation 
since the quark fields coupled to color singlets carry different momentum fractions. 

To illustrate that the color octet combination need not be small let us consider a 
three-quark system, as is done in constituent quark models. Irrespective of the details in 
the model, the color part of the three-quark wave function is ej^i- The color structure of a 
two-quark distribution is thus given by 



Fjj',kk' OC €jkl €j>k'l = 6kk' — Sjk' 5k 



(2.111) 



where I is the color index of the spectator quark and therefore summed over. With (2.104) 
one readily finds ^F = —^/2 (^F). The combination in (2.106), where the quark lines {13} 
and {24} are coupled to color singlets is then ^ (-v/8 ^F + ^F) = — (^F) and thus as large as 
^F itself. 

The preceding expressions can easily be adapted for the quark-antiquark distributions 
Fai,a2 defined in (2.87). With color indices labeled as in figure 5, the corresponding matrix 
element reads {{{q2,k^a2'l3,k') (Qij'^aiQij))) and is decomposed as on the r.h.s. of (2.103) 
with interchanged indices k and k' . An extra minus sign appears in the transformation 
laws (2.109) and (2.110) whenever a label A.q is changed to Aq, because T^g = — Fa^. 

An analogous color decomposition can finally be made for the interference distributions 



defined in (2.92), 



Ijj',kk' 

with 



(( {q2,kTa2 q^J') {q3,k' Tai qi,j) )) = 



1 6jk'Sj'k + 



2N 



St j.a j.a 
- ^ ^jk'T^j'k 



(2.112) 



I — 5k' j 5kji Ijj',kk' 

2N 



(((92ra2g4) (gsTai ^l) )) , 

2N 

"{q2Va2tyi){qzra,t''qi)\ 



In analogy to (2.106) one can alternatively use together with 
5k' j' 5kj Ijj',kk' = {{iq2,j^a2q^,j') (^S^'Tai 91,, 



1 

N 



(2.113) 



(2.114) 



By the same transformation as in (2.109) and (2.110), with appropriate sign changes for 
the antiquark matrices , one can rewrite this as a linear combination of matrix elements 



= (( (93 ^3,2 94) {q2 Tai 91 ) » 



(2.115) 
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where the quark bihnears have no uncontracted spin or color indices. Using the relation 
^'jk'^fk ~ k^jk^k'j' — jN^jk'^j'k One Can also rewrite (2.112) as 



Ijj',kk' = — 



1 

N2 



_ I{^jk'5j'k - Sjk6j'k') + y 2{N + 1) ^^^i^'^i'^ + SjkSj'k') 

(2.116) 



with 

The transformation between ®/ and ^I, is orthogonal. For = 3 we can rewrite 
^jk'^j'k ~ ^jk^j'k' = ^jj'i^k'ki ^nd recognize that describes the case where the quarks 
with momentum fraction xi are coupled to a color antitriplet, whereas ®/ describes the 
case where they form a sextet. 

2.3.2 Gluons 

The color structure for multi-gluon distributions requires the coupling of several color octets 
and is hence more involved than for quarks. For a two-gluon distribution we proceed by 
first coupling each of the gluon pairs {14} and {23} to irreducible representations and then 
coupling these pairs to an overall color singlet. For the color structures that can mix with 
quarks we write 

F-''^^' = (xip+)-i(x2P+)-i {{{Gtu^.Gl) (G2'n,,G?))> 
1 



(iV2 - 1)2 



■^F 5°-°-' 5^^' - ^^"^ ~ ^ /""'^ f^''''^ + ^^-^"^ ~ ^ Sp ^aa'c^bb'c _|_ 



(2.118) 



with a shorthand notation C^'lIa^G" = Hii^^ Qo-' Qo-^+i for the contractions of gluon 
polarization indices. As is readily seen from 

(xip+)(x2P+) 'F = {{{G\lia,G\) (G2n,,G?))> , 



{xrp+){x2P^)^F = - ^ (( {r'^'Gi Ha, Gl) (/-'Gf n,, G^) )) , 

{x,p+){x2P+)'F = ^]^r/ «K^'''g|n,,G^2) {d''''''GiUa-,Gt))) , (2.119) 

each of the pairs {14} and {23} in ^F, ^F and ^F is respectively coupled to a singlet, 
an antisymmetric and a symmetric octet. For hard-scattering processes producing color 
singlet states, these distributions enter the cross section as 

[^F'^F + ^F^F + ^F^F + ■■■] /{N^ - 1)2 . (2.120) 
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The ellipsis in (2.118) and (2.120) stands for terms where the gluon pairs are in higher 
color representations. For SU{3) these are 10, 10, 27, and the full decomposition reads 



aa'bb' 



1 

7V=3 64 L 



§0.0.' gbb' 



V8 



Ap raa'c rbb'c r 

3 5 



aa'c jbb'c 



+ 



10j7 ^.aa',bb' _|_ 2 



Jp f 



10 771 (,aa' 



M'\* 



) + 



27p ,aa' ,bb' 



(2.121) 



with tensors [89] 



,aa',bb' 
''10 

aa'M' 



gabga'b' 
jrabca'b 



27 



gab'ga'b 
ab' ra' 6 



2 jaa'c jbb'c 
l^aa'gbb' _ 



•abc fa'b'c _|_ jabc^a't'c-^ 



G^aa'c^bb'c 



(2.122) 



In ^'^F the indices {aa') are coupled to 10 and {bb') to 10, whereas in ^^F the opposite is 
the case. The normalization factors in (2.121) are such that the production of color singlet 
particles involves the combination ['^F'^F + ^F^F+^F^F + '^^F^F+^F'^'^F+'^'^F^'^F]/6A. 
Useful relations between the / and d tensors can be found in [90] . 

We conclude this section with the color decomposition of the quark-gluon distributions 
introduced in (2.102). The quark lines can only couple to a color singlet or octet, which 
has to be matched by the gluon lines in order to obtain an overall singlet. A complete 
decomposition is thus given by 



jpaa 



^ ^^F 5""' -^fV2 if^"'^ fjj, + 



iV(iV2 _ 1) 



2iV2 



Ar2 



(2.123) 



with 



(xip+) 'F = (( (g-3r„, q2) {G2Ua, G?) » , 



2A^2 



Ar2 



{{{q3Vl2)(.d''''''G'':Ua,G'l) 



(2.124) 



The factor i in (2.123) has been chosen so that "^F is real valued (since if^o'c jg Hermitian 
w.r.t. the indices a and a'). The normalization factors multiplying "^F and ^F are the 
geometric means of their counterparts in (2.103) and (2.118). 

2.4 Power counting and dominant graphs 

In section 2.1.3 we have already compared the power behavior in A/Q of single and multiple 
hard scattering cross sections. We now take a closer look at this issue and extend our 
analysis to the interference of single and multiple scattering. 

As building blocks for establishing the power behavior of the cross section we take cor- 
relation functions involving n parton fields and amplitudes T^^m for hard-scattering 



-32- 



processes with k incoming partons and m final-state particles. The relevant parton corre- 
lation functions are obtained by replacing the scalar parton fields in (2.1) by quark fields 
q and q, or by the transverse components of the gluon potential as in (2.96). To treat 
quarks and gluons on a common footing, it is convenient to use modified correlation func- 
tions that are divided by for each quark or antiquark line with momentum I, and 
modified hard-scattering amplitudes T'j^^^ that are multiplied with the corresponding fac- 
tor Furthermore, pairs qji and of quark fields in are contracted with one of 
the Dirac matrices Fq in (2.81) that give the dominant contributions to the cross section.^ 
The products T'^^^ T'j^i^^ modified hard-scattering amplitudes with their complex con- 
jugates are to be contracted with the corresponding Dirac matrices specified in (2.78). 
Since 

^ - ^ - " ^ - +^759, A^A^ (2.125) 



have the same mass dimension and are invariant under boosts along the z axis, the power 
behavior of the modified correlation functions is 

~ (2.126) 

regardless of the parton species. The power on the r.h.s. is just the mass dimension of 
By definition, all internal lines of the hard-scattering subgraphs are off shell by order Q^, 
so that the power behavior of the amplitudes Tl_^^ (where the propagators of external 
particles are truncated) is also determined by their mass dimension. For the processes 
considered in the following, one has 

n^^-Q'-'-^, (2.127) 

as can readily be checked for the example graphs in figure 8. 

For definiteness we consider the production of two particles with large masses and 
respective four-momenta qi,q2- Examples are the weak gauge bosons W, Z or a Higgs 
boson. The power behavior of the cross section is the same if we replace one or both of the 
heavy particles by a set of light particles such as a lepton pair or a pair of jets, provided 
that we integrate over the internal phase space of the final-state particles while keeping qi 
fixed. Replacing for instance a particle with momentum qi and mass Mj by two massless 
particles with momenta pi and p2, we have to change 

d% 27r6iqf - Mf) T^, Tj^U, = dx, dxi d^q, TrdlxiXi - ^ T^^U, (2.128) 



>2 ' 



into 

(2^ J '2^ -Pi- P2) Tk^2 = dxi dxi d J T^^a ^fc'- 



(2.129) 



*For the purpose of power counting, it is not important which of the matrices Ta is taken and which 
pairs of quark fields are contracted together if there are more than two of them. We will not specify these 
details in the present section and use <3?^ in a generic sense. Likewise, color indices are not relevant for 
power counting and will be omitted. 
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Figure 8. Example graphs and power behavior for different combinations of single and double hard 
scattering contributions to gauge boson pair production. It is understood that internal lines of the 
hard-scattering subgraphs are off shell by order Q^, whereas partons emerging from the proton 
matrix elements are off shell by order A^. 

where dOi is the solid angle of pi in the rest frame of qi. We recall that x^ = qf /p'^ 
and = qX Ip~ are defined in terms of final-state momenta and thus directly observable. 
According to (2.127) the scaling behavior of the squared hard-scattering amplitudes changes 
by l/Q^ when going from (2.128) to (2.129), which is compensated by the phase space 
volume sdVL\ ~ . One may put restrictions on the phase space integration, such as a 
minimum transverse momentum of pi, as long as d^\ remains of order 1. For each further 
final-state particle, the squared amplitude acquires an extra l/Q"^-, which is compensated 
by an extra phase space integration with volume of order Q^. 

After these preliminaries we can establish the power behavior of the conventional 
mechanism with a single hard scattering, shown in figure 8a. The cross section formula 
can be obtained in exactly the same way as in section 2.1.2. Omitting all factors that are 
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not relevant for power counting (including the 6 functions constraining (2.128)) we 

have 
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where for simplicity we have not displayed the momentum arguments of T', and ^' 
(which can readily be inferred from figure 8a). It is understood that in the second step we 
have made the collinear approximation and neglected I and I in the hard scattering, as well 
as l~ compared with l~ , and compared with l'^. The power behavior of the integration 
regions is d'^l ~ and dl^ ~ dl+ ~ A^/Q, so that together with (2.126) and (2.127) we 
obtain 
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For the double hard-scattering contribution in figure 8b one has 
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(2.131) 
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Note that we have used the constraint d{qf — l'^ — l'^) to fix the large component l'^ 
at its value qf in the collinear approximation, thus leaving the integral over the small 
component . If instead one uses the constraint to fix l'^ = Qi — one would have to 
count the integration element df^ as order A^/Q since l'^ can only vary by that amount. 
An analogous remark applies to the constraint J(gj~ — — l'^). The power behavior of 
(2.132) is 



sda 



Y\i=i dxi dxi d'^qi 



fig. 8b 



A2Q2 



(2.133) 



and hence the same as for single hard scattering, in agreement with the result we obtained 
for scalar partons in section 2.1.3. 

Let us now see how the power behavior changes if on one side of the final-state cut the 
two quark-antiquark annihilation graphs are connected by a hard gluon. We then have an 
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interference between double hard scattering and a single hard-scattering process as shown 
in figure 8c, 
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(2.134) 



This is power suppressed compared with the contributions in figures 8a and b and may 
therefore be neglected. It is instructive to see why the power counting changes between 
(2.132) and (2.134). Compared with T!^%iT!^_^^, the hard-scattering amplitude T'^_^2 is 
down by a factor of l/Q^ , which in the example of figure 8c is due to two additional quark 
propagators and one additional gluon propagator relative to figure 8b. The additional 
loop integrations over the large components l'^ and in (2.134) each scale like Q, but 
for the transverse momentum integrations one now has d%dH\ ~ A^ compared with 
dH\ d% (5(2) (q^ - l[ - l[) ~ A2 before. Altogether one has thus lost a factor of A^/Q^. 

By an analogous argument one finds that the differential cross section for the pure 
single hard-scattering mechanism in figure 8d is power suppressed by a factor of A^/Q^ 
compared with the one in figure 8c. 

The graphs in figure 8c and d involve single hard scatters with four incoming partons. 
There is, however, also an interference between double hard scattering and single hard 
scattering with two incoming partons. This involves correlation functions for three partons, 
of which at least one must be a gluon due to fermion number conservation. An example is 
shown in figure 9a, which gives 
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This is the same power behavior as the squared single and double hard scattering contribu- 
tions in figures 8a and b, so that interference terms of this type are not power suppressed. 
The example graph at hand has a suppression by Og since two-gluon fusion into gauge 
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Figure 9. As figure 8, but with parton correlation functions tliat involve gluons. 

bosons only starts at one- loop level, but for other final states like jets there is no such sup- 
pression. We will encounter these interference terms again in section 5.2.1 (see figure 38). 

Adding a hard gluon between the two single scatters on the left of figure 9a leads to 
the interference between different single hard-scattering processes in figure 9b. In the same 
way as above one finds that it is power suppressed by A^/Q^ compared with the leading 
contributions to the cross section. 

The contributions discussed so far have hard-scattering subprocesses with the same 
number of incoming partons from one and the other proton. One can, however, also have 
a parton in one proton scatter on two partons in the other proton. Examples for this 
are shown in figures 9c and d, and one finds that their power behavior is the same as for 
figure 9b. 

The pattern emerging from the preceding examples is clear: leading-power contri- 
butions are obtained as long as all hard-scattering processes involve only two incoming 
partons. This includes contributions from single scattering, double scattering and their 
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interference. For each hard scattering initiated by four partons one has a suppression 
by h? /Q"^^ and each hard scattering initiated by three partons comes with a suppression 
factor A/Q. 

2.4.1 Cross section integrated over transverse momenta 

So far we have considered the cross section differential in and q2- We now discuss how 
the power counting is changed when the cross section is integrated over these transverse 
momenta. As we aheady observed in section 2.1.3, the integration measure (Pqid?q2 
counts differently depending on whether Qi and q2 are both restricted to be of order A, 
or whether they can individually be of order Q and only their sum Qi + ^2 is restricted to 
size A. The latter requires a single hard-scattering process in both the amplitude and its 
conjugate. For our examples we thus find an integration volume d?qi (fiq2 of order A^Q^ 
for graphs 8a, d and 9b, c, d, whereas in the other cases the integration volume is of order 
A^. The resulting power behavior of the cross section is given in the figures. 

We see that the pattern of power suppression is different from the one we found for 
the cross section differential in and particular, the leading-power contribution 

now comes only from the standard single hard-scattering in graph 8a. 

The power behavior of the other contributions can be made more transparent by taking 
a closer look at the correlation functions they involve. As is evident from (2.130), the single- 
hard-scattering contribution of graph 8a goes with the transverse-momentum integrated 
correlation function f dl" d'^l and its counterpart for which are proportional to the 
usual collinear quark or antiquark densities. By contrast, integration of (2.134) over 
and ^2 gives a four-parton correlation function 



and its counterpart for $4. In these correlation functions all independent transverse parton 
momenta are integrated over, and correspondingly all field operators have the same trans- 
verse position. In physical terms, the single hard scattering in the conjugate amplitude 
has forced all hard scatters in figure 8c to take place at the same transverse position.^ By 
contrast, the double hard scattering contribution in figure 8b has two pairs of fields with a 
relative transverse distance y as we have seen in section 2.1.2, corresponding to two hard 
scatters taking place at positions that can be separated by a typical hadron size. This dif- 
ference has recently been pointed out in [66]. One obtains the same twist-four correlators 
(2.136) when integrating the contribution of graphs 8d and 9d over and ^2- 

^As is well known, integrals over transverse parton momenta in the correlation functions are logarith- 
mically divergent. If these divergences are avoided by a transverse- momentum cutoff of order Q (which is 
the largest scale in the process) then the relative transverse positions of the partons are of order 1/Q. 




(2.136) 
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Similarly, one finds that the transverse-momentum integrated cross sections from 
graphs 9a, b and c involve correlation functions 



J dq dq d'^q d^q $3 

^ [ dC^ d^^ e^^^'t+i^Pt (^p\Aqo) q{C2)raq{^i)\p) . , (2.137) 

which are proportional to collinear twist-three distributions. Again, a single hard scattering 
in the amplitude or its conjugate is enough to put all fields at the same transverse position. 

The power behavior like A^/Q^ of the integrated cross sections for graphs 8c and d is 
now readily understood, as it involves a collinear twist-four distribution for both colliding 
protons, each of which is responsible for a power suppression by A^/Q^. Likewise, graphs 9a, 
b and c involve the product of two collinear twist-three distributions, and graph 9d the 
product of a twist-two with a twist-four distribution. In both cases the integrated cross 
section is therefore suppressed by A^/Q^ (which happens to be the same suppression factor 
as for the double-hard-scattering contribution of graph 8b). 

In the transverse-momentum integrated cross section, graphs 8c and 9a with a double 
hard scattering in the amplitude play no particular role compared with their respective 
counterparts, graphs 8d and 9b, which involve the same correlation functions and have the 
same power behavior. Indeed, one may regard graphs 8c and 9a simply as higher-twist 
contributions with disconnected hard-scattering graphs on one side of the final-state cut, 
rather than associating them with multiple hard scattering. This was recently advocated 
in [66]. 

We emphasize that such a view is appropriate only if the cross section is integrated over 
transverse momenta. For observed transverse momenta Qi and 92 we have a different power 
behavior for graphs 8c and d, as well as for graphs 9a and b. In particular, the interference 
contribution from graph 9a then has the same leading-power behavior as graphs 8a and b. 
Let us also note that for graph 9a the quark and antiquark in each proton are not at the 
same transverse position for fixed Qi and 92 • If we express the correlation functions and 
$3 through matrix elements {p\AqO) q{£,2)^aQ{Ci)\p) and {p\A'^{0) q{^i)Taq{C2)\p) then 
the transverse- momentum integrations in (2.135) can be carried out and give 

j d% d% 5(2) {qi -h- U) e-'^'l^-^^^ = (27r)2 (^. - e-^^'"?' (2.138) 

for i = 1,2. With IqJ ~ A we thus have a typical quark-antiquark distance |^;^ —^2! ~ 1/A. 
2.5 Effects at small x 

Typical values of Xj and Xi at the LHC can be quite small, as we already noted after 
(2.20). At ^/s = 7TeV and qf = m?^ one has for instance \J XiXi = 1.3 x 10"^. Although 
phenomena at small x are not the main focus of this work, we wish to make a few comments 
on them in the present section. 

We begin by recalling that the usual densities for quarks, antiquarks and gluons rise 
with small x. This rise can be approximately described by power laws q{x) ~ q{x) ~ 
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x~^~"^'' and g{x) ~ with exponents Xg and between and 1 that depend 

on the factorization scale n. The abundance of small-x partons can be understood as a 
consequence of repeated radiation, which is essentially described by ladder graphs. Such 
graphs are in particular resummed by the DGLAP evolution equations, which make the 
rise at small x steeper when fi is increased. 

In the simple approximation where correlations between partons are neglected, multi- 
parton distributions are the product of single-parton densities as discussed in section 2.1.5. 
The distribution of n quarks or antiquarks then approximately behaves like F{xi,ki,yi) ~ 
{xiX2 ■ ■ ■ if all momentum fractions are sufficiently small. If all momentum frac- 

tions are of similar size, Xj ~ ~ x, this gives a factor x^^^^^^'**'?) in the cross section 
(2.35). If the same final state is produced by a single quark-antiquark annihilation, the 
corresponding factor is only according to (2.42).^ The multiple scattering mech- 

anism is thus enhanced for small momentum fractions, both for observed and integrated 
transverse momenta q^. In terms of graphs, this enhancement can be traced back to mul- 
tiple ladders, one for each pair of partons with the same momentum fraction Xi in figure 2. 
We expect that such an enhancement exists, although the above estimate based on com- 
pletely uncorrelated partons is likely too simplistic. 

Note that a strong rise at small x is only observed for parton densities that mix with 
gluons under evolution, but not for combinations like q{x) — q{x) or u{x) — d{x), which 
rise more slowly than x~^. A corresponding pattern is expected for multi-parton distribu- 
tions. Since they cannot mix with multigluon distributions, the interference distributions 
in fermion number or quark flavor discussed at the end of section 2.2.1 are not enhanced 
at small x. We hence expect them to play a minor role in small-x kinematics. 

The preceding arguments apply to both quark and gluon distributions in the framework 
of hard-scattering factorization, and based on the experience with single-parton densities 
one expects them to be relevant for momentum fractions Xj ~ Xj of order 10~^ or smaller. 
At very small x the gluon is by far the dominant parton species in the proton, and one 
may use high-energy factorization and the BFKL approach to describe the dynamics of 
gluon ladders. The primary expansion variable of this approach is log^, rather than the 
ratio Q/A used in the power counting arguments on which hard-scattering factorization 
is based. Basic quantities in high-energy factorization are Green functions depending on 
transverse gluon momenta, which bear close resemblance with the transverse-momentum 
dependent gluon distributions discussed in this work and naturally allow one to keep track 
of transverse momenta in the final state. 

Investigations of multiparton scattering in the BFKL approach can be found in [91-95]. 
In agreement with the arguments given above, one finds that the two-gluon distribution 
receives a contribution from two independent BFKL ladders, with a small-x exponent twice 
as large as for a single BFKL ladder [92]. More complicated graphs with four gluons in 
the t channel have been analyzed in [92, 94, 95]. As to the high-energy behavior of three t 

®For this comparison it is important that the hard-scattering cross section on the r.h.s. of (2.42) depends 
only on the momenta qi and not on p or p. It is hence proportional to Q^^" (without any further factors 
of Xi or Xi) and thus of the same order as the product (Ti(T2 ■ ■ o"n in (2.35). 
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channel gluons, all solutions found so far for the corresponding evolution equations have a 
weaker small-x growth than a single BFKL pomeron [96, 97]. 

Let us see how small-x dynamics affects the different graphs investigated in section 2.4. 
In line with our above discussion, we assume that correlation functions for four partons give 
a faster growth of the cross section with energy than correlation functions for two partons. 
We have no definite expectation for the corresponding behavior of three-parton correlation 
functions; if the results just mentioned for gluons in the small-x limit are a guide, then 
contributions with three t channel partons are not favored for small momentum fractions. 
For the cross section differential in we then find that the multiple-scattering mechanism 
of figure 8b is actually favored over the single-scattering graph 8a, which has the same 
power behavior in A/Q but a weaker rise at small x. With the caveat just mentioned, the 
interference contribution in figure 9a is expected to be less important. 

The cross section integrated over is dominated by the conventional single-scattering 
mechanism in figure 8a by power counting in A/Q. Among the contributions that are 
suppressed by A^/Q^ the double-scattering graph 8b is enhanced at small x. To a lesser 
extent the same is true for graph 9d, which involves four t channel partons in only one of 
the two protons. There may be situations where the small-x enhancement overcompensates 
the power suppression by A^/Q^, for instance in the high-energy production of minijets, 
where the hard scale Q is not too large. In such cases the BFKL approach may be more 
adequate than the one using hard-scattering factorization. 

2.6 The "effective cross section" 

The cross section for double hard scattering is often written as ads = cr2/{Cacs), where 
(7i and a"2 are single hard scattering cross sections, C is the combinatorial factor introduced 
below (2.21) and a^s is an "effective cross section" characterizing the strength of multiple 
interactions. Let us see to which extent such a formula holds true in the light of the results 
we have derived so far. 

Under the assumption that there are no correlations between different partons in the 
target hadron we derived the factorized form (2.70) for multiparton distributions in a model 
theory with scalar partons. This derivation carries over to the color singlet distributions 
of two unpolarized quarks, antiquarks or gluons, i.e. to ^Fq^^g^, ^Fg^^^j, ^Fq^^g^, ^-^91,92 ^^'^ 
^Fg^g, where the two quark fiavors qi and q2 may be different. If one further assumes that 
the impact-parameter dependent distributions of a single quark, antiquark and gluon have 
the form fc{x; h) ~ F{h) fc{x) with a common impact parameter profile F{b) for all parton 
species c, then the cross section (2.71) for double hard scattering takes the form 
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and 



foff 
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(2.141) 



The second form of (2.141) has recently been given in [98] and uses that the Fourier 
transform F{r) of F{b) depends only on because of rotation invariance. 

It is natural to ask whether (2.139) extends to the cross section differential in q^. If 
one has fc{x, z; b) ~ F{b) fc{x, z) for all parton types c then the factorized approximation 
in (2.69) gives Fc-,,c2{xi, Zi,y;b) ^ F{b+ ^xiZi) F{b + y - \x2Z2) fc2{x2, Z2) fc^{xi, zi). 
Inserting this into the cross section formula (2.36) does not lead to a factorized form, 
because XiZi appears in the arguments of the impact parameter profile F. A simplification 
occurs however if the measured transverse momenta are large compared to a hadronic 
scale A (while being much smaller than Q)7 The typical values of Zj in the cross section 
are then small compared with 1/A, whereas typical values of b and b + y are of hadronic 
size. We can thus approximate F{b + ^xiyi) F{b + y— ^X2Z2) ~ F{b) F{b + y) and obtain 



and (Jeff as in (2.141). Both (2.139) and (2.142) can be made differential in further variables 
describing the sets of particles produced by the two hard scatters. If one integrates these 
relations over kinematic variables in the presence of cuts, they only retain their validity if 
each cut refers to particles in one of the two sets but not in both. 

The assumptions that allow one to relate the cross sections for double and single hard 
scattering by a single process independent constant cJeff are quite strong, and a number of 
effects can invalidate (2.139) and (2.142): 

• an impact parameter profile Fc {b) that is not the same for different parton distribu- 
tions. The effect of this was estimated for a specific model in [75]. 

• a correlation between the x and b dependence in the single-parton distributions 
fc{x; b) or fc{x, z; b). Evidence that such a correlation is appreciable for x above 0.1 
comes from the calculation of the Mellin moments J dx x^~^ [fqix] b) + (— l)"/q(a;; 6)] 
with n = 1,2,3 in lattice QCD, see [99, section 4.4.5] and references therein. The 
interpretation of HERA measurements [100, 101] for jp — )■ J/^ p in terms of gener- 
alized parton distributions shows that the average squared impact parameter (b^) of 
small-x gluons in the proton has a weak logarithmic dependence on x [72, 102, 103]. 
An estimate of how a correlation between x and b in fc{x;b) affects multiparton 
interactions has been made in [104]. 

^This kinematic region is examined in detail in section 5. Notice the change of notation compared with 
the previous sections, where A denotes a hadronic scale or the size of |qj, whichever is larger. 
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• correlations between different partons in the proton, which invahdate the relations 
(2.69) and (2.70) between two-parton and single-parton distributions. In [105] it was 
argued that such correlations are significant. 

• an appreciable size of multiparton distributions that describe spin correlations be- 
tween two partons (section 2.2), of distributions where partons with the same mo- 
mentum fraction Xi are not coupled to color singlets (section 2.3), or of interference 
distributions in fermion number or quark flavor (section 2.2.1). 

Finally, the assumption that the observed cross section is given by contributions from 
either single or double hard scattering is invalid if their interference (see figure 9a) is 
important. All in all, we feel that (2.139) or (2.142) may be useful for order-of-magnitude 
estimates but should be used with great caution. Of course, one may define a^s as the ratio 
{dai/ dVi) {da2/ dV^) / {S da^isl dVidV2) of differential cross sections for single and double 
scattering. Since this ratio can depend on the process and on all kinematic variables, a^s 
is then not a universal constant. 

3 Beyond lowest order: factorization and Sudakov logarithms 

So far we have analyzed the lowest-order graphs that contribute to multiple scattering 
processes. For a systematic treatment in QCD we need to go beyond this approximation 
and in particular take into account graphs where additional gluons are exchanged. A 
complete analysis should eventually establish whether an all-order factorization formula 
can be written down for a given observable. We will not attempt to do this here, but 
provide some building blocks for such an analysis. We use the framework of hard-scattering 
factorization, which essentially organizes the dynamics according to virtualities (as opposed 
to high-energy or small-x factorization, where the organizing principle is based on rapidity). 
We focus on the cross section differential in small transverse momenta and in particular 
investigate the structure of Sudakov logarithms. In section 3.5 we will make some remarks 
on transverse-momentum integrated cross sections, described by collinear factorization. 
For reasons given in section 3.2 we will concentrate on the double Drell-Yan process, i.e. 
on the production of two electroweak gauge bosons, which for definiteness we take to be 
virtual photons. Likewise, we will use the single Drell-Yan process as an example when we 
recall the ingredients for factorization with a single hard scattering. 

3.1 Dominant graphs 

One of the first tasks when establishing factorization for a given process is to identify the 
dominant graphs in the kinematic limit one is interested in. The appropriate tool for hard- 
scattering factorization is the method of Libby and Sterman [106, 107], which we briefly 
recapitulate. The first step is to trade the limit of large kinematic invariants (which we 
collectively denoted by Q earlier) for the limit of vanishing masses of all partons. In doing 
so, one uses that up to an overall normalization the quantities of interest depend on the 
ratio of Q and the masses. If we keep small transverse momenta in the differential cross 
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section, those must be sent to zero as well. One is thus led to examine which graphs and 
which phase space regions give rise to mass divergences. Such divergences come from the 
poles of Feynman propagators, but only if for a suitable loop integration variable there are 
poles on both sides of the real axis, which "pinch" the contour of the loop integration. If 
there is no pinch, the poles can be avoided by deforming the integration contour. One finds 
that lines that give pinch singularities are either soft (i.e. all their momentum components 
are close to zero) or collinear to one of the incoming or outgoing particles of the process. 
All other lines are far off-shell (possibly after complex contour deformation). The leading 
contribution in the large Q limit of a given graph comes from regions of phase space in 
the vicinity of the pinch singular configurations just described. To obtain a factorization 
formula, one has to express subgraphs with collinear or soft lines in terms of matrix elements 
that make sense beyond perturbation theory. Parton densities and related quantities are 
constructed from these matrix elements. Off-shell lines are organized into hard subgraphs, 
which can be calculated perturbatively. 

A physically intuitive interpretation of the previous construction is provided by the 
Coleman-Norton theorem [108]. The pinch singular configurations of a graph correspond 
to a scattering process where the lines with collinear momenta correspond to classical 
trajectories in space-time. The trajectory associated with each line is proportional to 
its four-momentum, so that it shrinks to a point for soft lines. In the "reduced graph" 
that represents the corresponding classical process, off-shell lines in the original graph are 
likewise contracted to points. 

The preceding analysis is based on the denominators of Feynman propagators and 
gives only a necessary condition for the occurrence of mass singularities. A power counting 
analysis taking into account the numerators of Feynman graphs (similar to the one we gave 
in section 2.4) provides further restrictions on the contributions that actually dominate a 
given observable. At this level, the polarization of gluon lines is found to play a crucial 
role. 

For single Drell-Yan production at fixed small transverse photon momentum, one finds 
that the dominant graphs have the structure shown in figure 10a. For each of the colliding 
protons there is a collinear subgraph. On either side of the final-state cut there is one 
hard subgraph producing the final state boson, connected with each collinear subgraph 
by exactly one fermion line and an arbitrary number of gluon lines, which must have 
polarization in the plus direction for right-moving and in the minus direction for left- 
moving collinear gluons. Finally, there is a soft subgraph with soft gluons attaching to 
either of the collinear subgraphs. There are no soft gluons coupling to the hard subgraphs. 

The dominant graphs for double Drell-Yan production are easily identified and just 
have an additional hard subgraph for the second produced gauge boson on either side 
of the final-state cut. As we have already seen in section 2.4, hard subgraphs that are 
connected to each collinear graph by a single parton line have leading power behavior. The 
power counting for the soft graph is not affected by having one or two hard subprocesses. 

*This was not stated in the original work by Libby and Sterman, who considered transverse-momentum 
integrated quantities. 
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Figure 10. Leading graphs for the single (a) and double (b) Drell-Yan process at measured qr ^ Q- 
The upper and lower blobs denote coUinear subgraphs, the blob crossing the final-state cut (dashed 
line) denotes a soft subgraph, and the blobs with a final-state gauge boson denote hard subgraphs. 
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Figure 11. Space-time representation of the leading graphs for the single (a) and double (b) Drell- 
Yan process at measured qT ^ Q. The parton lines move along light-like paths and have been 
drawn with a slight curvature only for clarity. Likewise, the two bosons in figure (b) are meant to 
be produced at the same point in t and z. 

Finally, the absence of soft gluons coupling to a hard subgraph has the same reason as 
in the single Drell-Yan case, namely that such soft gluons increase the number of hard 
propagator denominators in the hard subgraph, without providing a compensating large 
numerator factor or phase space volume. 

The space-time representation in the sense of the Coleman-Norton theorem is shown in 
figure 11 for the graphs in figure 10. Parton lines from one and the other proton meet at one 
point in the t-z plane and annihilate into a gauge boson. For double Drell-Yan production, 
the two bosons are produced at the same point in t and z. The transverse momenta of 
partons and the produced bosons are neglected in this interpretation (see above), so that 
the classical scattering process takes place at fixed transverse coordinates (x and y). 

3.1.1 "Rescattering" contributions 

Before discussing in detail the leading graphs in figure 10, we wish to comment on graphs 
of the type shown in figure 12a. They have been associated with "rescattering" in the 
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Figure 12. (a) A graph associated with, "rcscattering" in the literature, (b) Space-time represen- 
tation for the production of two high-p^ partons in a single hard 2 — > 2 scattering process. The 
dashed lines have momentum components in the transverse x-y plane; all other lines move on the 
light-cone z — ±.t. 

literature [109, 110] and were calculated in terms of two hard 2 — )• 2 QCD processes, where 
the parton with momentum k in the figure is treated as an outgoing parton in the first 
scattering and as an incoming parton in the second one. It is understood that the transverse 
momenta P2 and p are all large (we denote their order by pT below). 

We argue here that this is not a correct way to calculate the graph, at least not 
within the usual hard-scattering factorization framework used in [109, 110]. According 
to our discussion in the previous section, the lines that enter or exit a hard-scattering 
subgraph must correspond to pinched singularities and thus admit a classical space-time 
interpretation in the sense of the Coleman-Norton theorem. This is not possible for the line 
with momentum k in figure 12a. As illustrated in figure 12b, the two partons emerging from 
a hard 2—7-2 process have large transverse momenta and, being on shell, thus have finite 
rapidities. In other words, their velocity in the z direction is smaller than the speed of light. 
As soon as such a parton has propagated over a finite distance, it can no longer scatter on 
another parton from one of the two initial protons, since those partons do move with the 
speed of light along z. The proper treatment of the parton with momentum k is thus to 
regard it as an internal line in a single hard-scattering process with three incoming partons 
(hjhj) and three outgoing ones {pi,P2,P)- As we saw in section 2.4 such a contribution 
is power suppressed (if Pi + P2 + p is integrated over, it involves a parton distribution of 
higher twist) and can hence be neglected. 

Put differently, the parton with momentum k is generically far off-shell in the leading 
region of the graph in figure 12a. A kinematical analysis readily shows that the final-state 
momenta Pi,P2,P fix the sum if -|- to a large value of order px, up to small corrections 
of order l/px- The value of if is however integrated over a large interval of order px- For 
a particular value of if in this interval, the propagator of k does have a pole, but this pole 
is not a pinch singularity (the gluons adjacent to k are far off-shell when /c^ = and their 
propagator poles are a distance of order px away in the complex if plane). One can thus 
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deform the integration contour of such that k"^ 



is always of order p'^ and thus large. ^ 



3.2 CoUinear and soft gluons 

We now return to the graphs in figure 10. They contain an arbitrary number of collinear 
and soft gluons, and further simplifications are required to obtain a useful factorization 
formula that involves a limited number of nonperturbative quantities. 

In existing factorization theorems, the effects of collinear and of soft gluons are de- 
scribed by Wilson line operators, to all orders in the strong coupling. The possibility to 
obtain such a simple structure is crucial for establishing factorization. Detailed analyses of 
this issue can be found in [48-51] for single Drell-Yan production or for its crossed-channel 
analogs, the production of back-to-back hadrons in e~^e~ annihilation or semi-inclusive deep 
inelastic scattering (SIDIS). By contrast, for hadron-hadron collisions producing back-to- 
back jets or hadrons with measured transverse momenta, serious obstacles to establishing 
factorization have been identified in [111] and in previous work cited therein. A systematic 
treatment of transverse-momentum dependent factorization for jet or hadron production 
in multiple hard scattering will probably need to wait until a suitable formulation for single 
hard scattering has been found. 

We therefore limit our considerations in this section to the double Drell-Yan process. 
Extending our arguments to the production of other colorless particles is trivial if the hard 
scattering is initiated by quarks or antiquarks and should be possible if it is initiated by 
gluons. We shall not attempt to give a full proof of factorization even for double Drell-Yan 
production. Instead, we will analyze the lowest-order graphs with an additional exchanged 
collinear or soft gluon. To a large part this will be a recapitulation of the corresponding 
analysis for the single Drell-Yan process. We nevertheless give the necessary steps in some 
detail, in order to see how the arguments generalize to double hard scattering. We will pay 
particular attention to the color indices for quarks and antiquarks, since the color structure 
of two-parton distributions is nontrivial compared with the single-parton case. Finally, we 
will point out which further issues need to be settled to obtain a full proof of factorization. 

3.2.1 Prom collinear gluons to Wilson lines in parton distributions 

Figure 13 shows an example where several gluons collinear to the right-moving proton 
p couple to a left-moving quark or antiquark. The quark or antiquark is thus taken far 
off shell, so that its propagator and its coupling to the gluon belong to one of the hard- 
scattering subprocesses. 

We now recapitulate the analysis of one such coupling, which is well-known from single 
Drell-Yan production, taking particular care of color indices and of the distinction between 
quarks and antiquarks. The relevant part of the graph in figure 14a can be written as 



^When actually calculating the hard scattering, one can nevertheless integrate along the real axis; 
the pole of + ie) then provides an absorptive part to the hard-scattering amplitude. The possibility 

to deform the integration contour of if justifies the perturbative treatment of the propagator for k. 




I 



(3.1) 
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Figure 13. Example graph for the double Drell-Yan process with coUinear gluons coupling to 
left-moving quarks or antiquaries before those undergo a hard scattering. 




a j a f 



Figure 14. Collinear gluons in the Drell-Yan process. Top row: subgraphs with a right-moving 
gluon coupling to a left-moving quark or antiquark before it annihilates. Bottom row: corresponding 
graphs after the off-shell propagators have been replaced by eikonal lines. 

where in a shorthand notation we write {■ ■ ■ qj A°'''^ . . .) and {. . .qj/ . . .) for the hadronic 
matrix elements of the right and left moving proton, respectively. The subscript c on I 
indicates the collinear approximation specified after (2.90), i.e. l~ = l~ , = and Ic = 0. 
Instead of the spinor u{lc) for the incoming quark we could also use the projection operator 
P{1), see the discussion after (2.79) and (2.90). The vertex with the produced photon and 
the spinor for the incoming right-moving quark are not needed for our argument and have 
been omitted. Our sign convention for the strong coupling g is such that the covariant 
derivative reads D'^ = d'^ + igA^ . 

The expression in (3.1) has the structure W^Ha, where R is the matrix element of 
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the right-moving proton and H a hard-scattering ampHtude. One therefore has \R~^\ ^ 
\R~\, \R\, whereas ah components of H°' are genericaUy of the same size. To leading-power 
accuracy we therefore have R'^Ha ~ R^H~. We now introduce an auxihary spacehke 
vector V with 

t;->0, v+<0, v = (3.2) 

and either |t;"'"| ~ or \v~^\ <^v~ . We can then write 

RH ^ R+H- = R+v^—^ e+H-^Rv-^ — (H , (3.3) 

where in the last step we have used the conditions on the components of -R", if" and 
just stated, as well as l^"*"! ^ \£\ for the momentum i of the right-moving gluon. For 
reasons given in the next section, we have provided an ie prescription to the factor 1 /iv in 
(3.3) such that the pole in is on the same side of the real axis as in the propagator of 
the off-shell quark that couples to the photon in figure 14a: 

£v + ie = i+v~ -r\v+\+ie, {£ + lcf + ie = 2£+l- + ie , (3.4) 

where it is important that l~ > 0. With (3.3) we can rewrite (3.1) as 

1 



Ta = {...qjA'''\..){-igt'^^,v^) 



iv + ie 



7(£ + /",) 



(7^)n(/, 



{...q,,...). (3.5) 



With = J {£ + Ic) — jlc and (jlc) u{lc) = we finally obtain 



Ta = {... qj . . .){-tgt';^,v^) -j^^ u{l) (. . . q,, . . .) . (3.6) 

In the hard-scattering amplitude we have thus traded the coupling —igf^^a of the gluon to 
the quark and the adjacent quark propagator ij^il + /c) for the coupling —igf^Va of the 
gluon to a so-called eikonal line and the eikonal propagator i/{iv + ie). 
Repeating the same steps for the graph in figure 14b gives 

n={...q,.. .) v{l) 7. ~\ i-ig) t-^, {. . . A-''^ qy . . .) 

= {...q,...) vile) {-tgty v^) {. . . A-''^ qy ...). (3.7) 

The change from an incoming quark to an incoming antiquark in the hard scattering has 
changed the overall sign of the propagator ij^il + /c)) which is reflected in an overall sign 
change of the eikonal propagator ijiiv -|- ie). On the other hand, the momentum flow in 
the graph and the resulting ie prescriptions have remained the same. 

It is instructive in this context to compare Drell-Yan production with SIDIS, where 
one has an outgoing quark or antiquark in the hard scattering. The corresponding graphs 
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3 « i' 



Figure 15. As figure 14, but for semi-inclusive deep inelastic scattering. 

are shown in figure 15a and b. Taking the same vector v as before, we have 

= (. . . A^-'^ . . .) i-ig) t^^n^ vile) (• • • Qf • • •) , 

= (. . . qj A"''' . . .){-igt';^,Vo^) v{Q (. . . qr ■■■). 

U = {■■■<!,■■■) m 7. i-^9) ty {. . . A-'^ q,. . . .) 

= {...qj...) u{Q ^ (-i<7i •j't'a) (• • • ^"•'^ qr ...). (3.8) 

Compared with graphs 14a and b, the relative flow of the momenta £ and Ic in the off-shell 
quark or antiquark has changed. Hence the corresponding propagator has a denominator 

{I - kf + ie = -2£+r + ie (3.9) 

instead of the one in the second equation of (3.4). As a result, the sign of ie in the eikonal 
propagator is now reversed. 

A graphical notation for eikonal lines needs to specify the flow of the momentum i 
relative to 

1. the color flow (and hence the fermion number flow in the quark line which is repre- 
sented by the eikonal line). This determines the overall sign of the eikonal propagator. 
We denote the color flow by an arrow on the eikonal line, which points in the same 
direction as the arrow on the original fermion line. 

2. the flow of the large momentum Ic in the original fermion line, which is either an 
incoming or an outgoing line in the hard-scattering subprocess. This determines the 
sign of ie in the eikonal propagator. We indicate this graphically by a full or an empty 
circle at the end of the eikonal hne, such that the large momentum flows from the 
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> O + ie 



< O iv + ie 



a, a 



iv — ie O > 



iv — ie O < 



igfv" 



Figure 16. Feynman rules for cikonal lines representing quarks or antiquarks. The rules for lines 
to the left of the final-state cut (denoted by the dashed line) directly follow from (3.6) to (3.8), and 
those for lines to the right of the final-state cut are obtained as usual by complex conjugation. 



full to the empty circle. Since incoming and outgoing partons in the hard scattering 
can be associated with a classical path in space-time according to section 3.1, the full 
circle represents the past and the empty circle the future time direction. 

The corresponding Feynman rules are given in figure 16, and the graphs corresponding 
to the eikonal representation in (3.6), (3.7) and (3.8) are shown in the bottom rows of 
figures 14 and 15.^" 

We now briefly review how eikonal lines are generated by Wilson line operators in 
the hadronic matrix elements that appear in a factorization formula. The relevant part of 
the expression (3.6), together with the relevant integrations over momentum and position 
variables reads 



Using the representation 



(2vr)4 



iv + ie 



iv + ie 



we can rewrite this as 



dX e*(^^+f-«)^ vA'^iOi-igf 



(3.10) 



(3.11) 



(2vr)4 
d\vA''{i-\v)ej^, 



3f 



(3.12) 



Our graphical notation differs from that in the literature. In [48] for instance, an arrow on the eikonal 
line was associated with the flow of the large momentum, and the overall sign due to the color flow was 
indicated by explicit color indices and taken into account in the vertex between a gluon and an eikonal line, 
rather than in the eikonal propagator. 
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Introducing the Wilson line 



W{C;v) = Pexp 



ig / dXvA^iC- >^v)t'' 







(3.13) 



where P denotes path ordering, we recognize the term in square brackets in (3.12) as the 
term of order g in the expansion of VF^(^; v). In a full factorization proof, one has to show 
that the coupling of two or more collinear gluons to the incoming quark line in figure 14a 
exponentiates, so that their combined effect is the replacement 

qAO^QjiO[w\^;v)]jf (3.14) 

in the operator defining the parton distribution. Likewise, the expression in (3.7) corre- 
sponds to the one-gluon term in the replacement 

qj{0^[WiC;v)],,,q,,iC). (3.15) 

The conditions we imposed on v after (3.2) hold in a frame where p moves fast to the 
right. One readily finds that in the rest frame of p one has > 0, so that the Wilson line 
(3.13) relevant for Drell-Yan production has a path pointing into the past. By contrast, 
the reversed sign of ie in the eikonal propagators for SIDIS corresponds to 

qriO gj^jg QjiO [W'HC; , qM) g^^jg t^'^^' ^^'^^^ ^^"^^^ 

with a future-pointing Wilson line 



W'{^;v) = Pexp 



/•CO 

/ dXvA''{C + Xv)^ 
Jo 



(3.17) 



The preceding discussion was for right-moving collinear gluons and generalizes trivially to 
left-moving collinear gluons in the proton with momentum p. The corresponding Wilson 
lines are to be defined with an auxiliary vector w that satisfies 

w^>0, w-<0, w = (3.18) 

and either \'w~ \ ~ or \w~ \ <^ in a frame where p moves fast to the left. In the rest 
frame of p one then has > 0. 

The manipulations in the preceding arguments are all concerned with a single hard- 
scattering subprocess at a time, so that they readily apply to double Drell-Yan graphs such 
as in figure 13, where they give the order g part of a Wilson line for each quark or antiquark 
operator in the multiparton distributions. The full operator for a two-quark distribution 
then reads for instance 

[q{-hz2)W\-^Z2;v)]k'ra, [W{^Z2;v)qC^Z2)]k 

x[q{y-^zi)W\y-^zr,v)]j>Ta,[Wiy + lzi;v)qiy + ^zi)]j , , _^ . (3.19) 

zJ=z{=y+=0 

The open color indices k,k' , which were carried by quark fields in the lowest-order for- 
mula, are now carried by the "ends" of the four past-pointing Wilson lines. The projection 
on color singlet and color octet distributions is done as in (2.103). 

Let us now mention how the previous arguments need to be generalized to obtain a 
complete factorization proof for double Drell-Yan production. 
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• The step from (3.5) to (3.6), which ehminates an internal fermion propagator in the 
hard-scattering graph, is elementary when applied to the lowest-order hard scattering. 
For more complicated graphs (with loop corrections or further external gluons) one 
needs a Ward identity to achieve this simplification. In a model theory with Abelian 
gluons, this is quite simple to establish, see e.g. [50, chapter 10.8]. The formulation 
for QCD is more complicated and involves external ghost lines in addition to external 
gluons in the hard scattering (see [50, chapters 11.3 and 11.9]). 

• We have considered only one gluon coupling to each hard-scattering subgraph. One 
needs to show that the coupling of an arbitrary number of gluons exponentiates and 
gives a full Wilson line W{y; v) or its complex conjugate. Again, this is simple to 
show for Abelian gluons (see [50, chapter 10.8]). To the best of our knowledge, an 
explicit proof for transverse-momentum dependent distributions in QCD has not yet 
been given. 

We note that the present and the previous point only concern one hard-scattering 
subprocess at a time. It should therefore be straightforward to extend arguments 
valid for the single Drell-Yan process to the case of double Drell-Yan production. 

• The two Wilson lines W{^Z2;v) and W{y + ^zi;v) in (3.19) correspond to gluons 
in the scattering amplitude, where all gluons fields should be time ordered. With 
v'^ < the gluon operators in one Wilson line have a spacelike separation, so that 
they commute and can readily be brought into the order required by path ordering. 
Two gluon operators in different Wilson lines do not necessarily have this property, 
and the possibility to reorder the fields needs to be investigated. A similar statement 
holds for the two Wilson lines W^—^Z2;v) and W^{y — ^zi;v) that correspond to 
gluons in the conjugate scattering amplitude. 

• The operator in (3.19) is not explicitly gauge invariant, because the Wilson lines 
end at different positions at infinity, namely at ai — oov with finite spacelike a, for 
i = 1,2,3,4. The same issue already arises for single-parton distributions and has 
been discussed in [113, 114] for lightlike Wilson lines, i.e. for v'^ = 0. In a gauge 
where the gluon potential (and any product of gluon potentials) has zero expectation 
value at a — oov, one can trivially complement the operator (3.19) with Wilson 
lines that go in the transverse direction and connect the lightlike Wilson lines to a 
common reference point, e.g. to —oov. After projecting the open color indices at this 
reference point onto color-singlet or color-octet combinations, the resulting operator 
is explicitly invariant under local gauge transformations. The extra Wilson lines in 
the transverse direction are essential in the gauge vA = 0, where the Wilson lines in 
(3.19) reduce to unity, see the discussion in [113]. 

As we will see in section 3.2.3, the choice = is not suitable for transverse- 
momentum dependent factorization. To obtain a gauge invariant definition of the 
relevant parton distributions, one needs to extend the procedure just described to 
the case where < 0. This holds both for single and multiple hard scattering. 
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Figure 17. Coupling of a soft gluon to a collinear partoii tliat (a) enters tlie liard scattering or (b) 
is a spectator. 

3.2.2 Soft gluons and the soft factor 

We now turn to the exchange of soft gluons between right- and left-moving partons, i.e. to 
the soft subgraph in figure 10, and show how it can be described in terms of a soft factor 
that is defined as the vacuum expectation value of Wilson lines. 

For definiteness we consider one soft gluon with momentum I, exchanged between the 
soft subgraph 5" and the collinear subgraph i?" of the right-moving partons. Here a is the 
polarization index of the gluon considered, and the indices for any other external gluons 
are omitted for simplicity. We assume that the components of i are of comparable size, 
l^+l ~ \£~\ ~ \£\, as well as the momentum components of all other soft gluons attached 
to S. The components of are then also comparable to each other. Since ji?"*"! ^ 
we then have 

£R w rR+ , SR « S-R+ . (3.20) 

Introducing an auxiliary spacelike vector w as in (3.18) with <C , we furthermore 
have Sw ~ S~w~^ , so that we can write 

+ a 

SaR^'^S-- — £-R+^Sa-„ £R (3.21) 

£-w+ + ^e £w + te ^ ^ 

with a factor w°' / (£w + ie) that will eventually turn into a Wilson line. The ie prescription 
for the pole at £w = is adequate for £ flowing from S into R in the scattering amplitude, 
i.e. on the left of the final-state cut in figure 10. We note that this prescription corresponds 
to the one for collinear gluons in the previous section, cf. figures 14 and 17a. 

The approximations in (3.20) and (3.21) break down in the so-called Glauber region, 
i.e. for soft momenta dominated by their transverse components, \£.\ ^ \£~^\, \£~\- A major 
part of a factorization proof for hadron-hadron collisions is to establish that this momentum 
region does not contribute to the final factorization formula. With the ie prescription we 
have chosen, the pole of l/{£w + ie) is on the same side of the real £~ axis as the propagator 
pole of the quark with momentum I + £ in figure 17a, which is readily seen by adapting 
(3.4). In the graph of figure 17a one can avoid the Glauber region by a contour deformation 
to complex i~ . With the same contour deformation one can however not avoid propagator 
poles in graphs where the gluon couples to a spectator parton (rather than to the parton 
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Figure 18. Graphical illustration of the Ward identity (3.22). Shown is only the part of the 
quark-antiquark distribution to the left of the final-state cut. The Feynman rules for eikonal lines 
are given in figure 16. 

entering the hard subprocess). An example is shown in figure 17b, where the pole in i~ of 
the propagator for the hne p — I — i is on the opposite side of the real axis than the pole 
of the propagator for the line I + i in figure 17a. To apply the Ward identities discussed 
below, one has to make the same contour deformation for i~ in both graphs of figure 17 and 
will thus pick up a residue contribution from the propagator of the spectator parton. For 
single Drell-Yan production one can show that the sum over all such residue contributions 
cancels due to unitarity, see the discussion in [50, chapters 14.3 and 14.4] and in the original 
literature cited therein. We do not know whether and how such arguments can be extended 
to the case of double hard scattering and leave this issue as an important task for further 
investigation. We will proceed under the assumption that such an extension can be made. 

Following the procedure for single Drell-Yan production, the next step in our argument 
is to use a Ward identity to relate the collinear subgraph with a gluon attachment to the 
same subgraph without a gluon. For the correlation function describing quark-antiquark 
emission and an additional gluon in the amplitude, this identity reads 

^lf,kA^-' h.kAA) = sw i-igt^j ^^rMh - ^, kAA) 

+ Sw{-igCk)j-^^<^,j,,rnk'ik,k-^,A,l2) (3-22) 

and is depicted in figure 18. Analogous identities can be written down for the emission of 
two quarks or two antiquarks, with a factor i/{iw + ie) for each quark line and —i/{iw + ie) 
for each antiquark line in the amplitude. We leave it to future work to give a general proof 
of these identities, but verify them here for two simple examples. 

Our first example is a quark-antiquark pair with a pointlike coupling to a target. The 
corresponding two-parton distribution is then proportional to 6jki{'yh)~^ ^ {—i){'yl2)~^, 
where li and I2 are the respective momenta of the quark and antiquark, and j and k are 
their respective color indices. The tensor product (8) refers to the spinor indices, whose 
coupling at the vertex with the target we need not specify for our argument. Attaching a 
soft gluon in the amplitude, we have to add the graphs in figure 19a and b. Contracting 
the gluon polarization index a with i" and using the same trick as in the step from (3.5) 
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Figure 20. Graphs for a gluon coupling to a three-quark system tliat originates from a colorless 
fermion via a pointlike vertex, j, k, I and a are color indices. 

to (3.6), we obtain 

ih -f{h-i) 7/2 7*1 7(^2 -«) 7*2 

„ ( i i } —i „ i ( —i —i 



l{h-i) 7^1 J 7^2 ^"ih \l{h-^) ih 

7(/i-£j 7/2 7/1 7('2-«j 

Multiplication with Sw/{iw + ie) gives (3.22) for this particular case. 

As a second example we take a colorless fermion target coupled to three quarks by a 
pointlike vertex. The two-quark distribution is then proportional to 

ejki ^ ^ u(k) , (3.24) 
7*1 7*2 

where ^3 is the momentum of the spectator quark. Coupling a gluon to this system, we get 
the three graphs shown in figure 20, which after contraction with give 

1 % % 

7/1 7(/i - 1) 7/2 

% % % 

7/1 7/2 7(^2 --t) 

+ (.-ig) ttm^ikm ^ ® ^ ® U{h) (7^) 



^'^'"7^1 7^2 ' l{h-^) 
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Figure 21. Diagram witli one soft gluon exclianged between tlie left- and riglit-moving partons to 
tlie left of the final-state cut. 



i 



£) 7/2 



1 1 

— n(/3) + i {-ig) ttm ^iml — ® 



7^1 



(3.25) 



i'j'm^niki + tkm^jmi + ^Im^jkm' ^his tensOT is zero, because it is 

= 0. 



with the tensor V^-^ 

completely antisymmetric and hence proportional to V^ji^^ij^ = ^(^jj + ^kk + 
Multiplication with Sw/{iw + ie) finally gives the equivalent of (3.22) for a two-quark 
distribution. 

The preceding arguments can readily be adapted for a soft gluon attached to left- 
moving collinear partons by exchanging -|- and — components of the relevant vectors. The 
auxiliary vector w is then replaced by v as in (3.2), with \v'^\ <C v~ . Likewise, one can 
repeat all arguments for soft gluons in the conjugate amplitude, i.e. to the right of the 
final-state cut in figure 10. In the corresponding Ward identities one then has to use the 
Feynman rules on the r.h.s. of figure 16. 

Consider now the diagram in figure 21, where in the amplitude one gluon is exchanged 
between the left- and right-moving partons. Its contribution to the cross section is propor- 
tional to 



(27r)4 (27r)' 



(2T:f5^'^\i + l)S%{i,I) 



i=i '' 



d%d%{2^)H^^\qi-l\-l[) 



X [*^ai,a2] '21 ^'l> ^2) [^ai,a2] jj',fcfc'(^; ^li ^2) ^'li ^2, 



d^l dH 



(27r)4 (27r)4 



(27r)V^)(£ + l) 



IW 



£w + ie 



IV 



/3 



£v + ie 
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2 

n / d% d% (27r)45(^) -k- k) J d% d% (27r)^j(4) (gi - l[ - l[) 

4 = 1 

( — igt'j^) [^ai,a2]mj' ,kk'{h "^1^21^11^2) ~ ("^fi'^mfc) [^ai ,02] jj'.mfc' (^1 ' ^2 "^1^15^2) 



(-Wtkn) [*^'ai,a2]ii',nfc'(^l, ^2 -^^1^2) " (-*5*ni) [*^*ai ,02] nj'.fcfc' (^1 " ^2 ' ^1 ' ^2) 



(27r)V4)(£ + l)(-i54 



X y dHi d% (27r)45(^)(gi - /i - /"i -£) J d% d% (27r)^<5(^)(g2 -h-h-l) 
X y d^/; d^/"; (27r)45(4)(gi - l[ - l[) 

,a2] mj' ,kk' (^ij ^2) ^'i) ^2) [^ai,a2]ij',nA;'(^i) ^2) ^'i! ^2) + {three more terms}, (3.26) 

where in the last step we have shifted the integration variables li and ^2- For simplicity 
we have omitted a global factor, as well as the expressions for gg — 7*, which in the hard- 
scattering approximation only depend on the external momenta qi and q2 and thus do not 
appear under the loop integrals (see section 2.1.2). 

To provide a representation beyond perturbation theory, we represent the soft subgraph 
(which for two external gluons is just the gluon propagator) as a matrix element, 



{27r)^6'^^\i + i)SU^,i)= / AA"e*^^+*«^0|^^(e)4(O|0) 



(3.27) 



Here we have omitted the time ordering between the fields, which requires justification 
when and £^ do not have a spacelike separation. We gloss over this point here (see also 
our discussion at the end of section 3.2.1) but return to it briefly at the end of section 3.3.1. 
Using (3.27) and (3.11) we then have for the first term in (3.26) 



00 />oo 



JO 



«^+^«^ / dX dXe'^^'"+'^^'' {-igt';j{-igtlJ{0\wA''{C)vA\C)\0) 



(27r)4 J (27r)4 

X / d% d% (27r)2 6{qt - /+) 6{q^ -li) [ d^^i e-^^i^^i-^i-fi-^) 



X / d% d% {2-nf 6{q+ - /+) 6{q^ -q) / d^^2 e-^^^fe-'a"'" 



X J d% d% (27r)2 8{qt - /'+) 5{q-{ - j d^$,\ e-'^iin-'i-^i) 

X [*^ai,a2]mj',fcfc'(^li ^2' 'l' ^2) [^ai ,02] (^1 ' ^2 ' ^1 ' ^2) 

J d^ii d^ii cP^'i e-^^i''i-^^292-i4'iqi 

X j d\ d% d% e^^l'l+*^^'2+■*^'l^^ (2vr)3 5(g+ - /+) 5iqt - it) Siqt - I'l^) 



ai,a2j mj 



'jkk'ihi ^2) ^1) ^2) 



-58- 



d% d% d% e*«i'i+^^2'2+*€'j; (2vr)3 _ 1-) 5{q- _ /"-) 5{q- _ l'{ 



X ['^'01,02] ji',nfc'(^l) ^2) ^1' ^2) 

fOO 



-i5 y d\wA''{iiT - \w)t 



a 

jm 



-ig / d\vA\i2T-\v)An 



|o>, 



(3.28) 



where S.^j- denotes the four- vector with = = and transverse components ^j. 
The corresponding expression for the diagram without soft gluon exchange is obtained 
by replacing the last line in (3.28) by 6jmSkn- Using (3.13), we recognize the factors in 
square brackets in that line as the order g terms in conjugate Wilson lines W^^(^it; w) and 
{^2T]v). In the transverse plane, the paths of these Wilson lines are at the positions that 
are Fourier conjugate to the transverse quark momenta li and I2 in (3.28). The three other 
terms in (3.26) give analogous contributions, with Wilson lines W{£,2T]w) and W{^,it,v) 
at the positions that are Fourier conjugate to the transverse antiquark momenta I2 and Zi, 
respectively. 

After a change to symmetric momentum and position variables as specified between 
(2.1) and (2.7), and after restoration of global kinematic factors, the second to fifth lines 
on the r.h.s. of (3.28) turn into the product Fa^^a2{^ii ^iiy)Pax,a2{xi-, z,i,y) of two-parton 
distributions in transverse position space, and the Wilson lines are to be evaluated at the 
appropriate transverse positions of the quark or antiquark fields in the definition of these 
distributions. 

It is straightforward to repeat the preceding derivation for a soft gluon exchanged to 
the right of the final-state cut, as well as for the case where the gluon crosses this cut. For 
an model theory with Abelian gluons, it is not difficult to see how soft subgraphs with an 
arbitrary number of external gluons add up to full Wilson lines, in close analogy to the 
case of single Drell-Yan production. We do not attempt here to give a corresponding proof 
for the nonabelian theory, given that even for the single Drell-Yan process this is quite 
involved. The structure suggested by our analysis of one-gluon exchange is however clear: 
the effect of all soft subgraphs is to multiply the Born-level cross section (2.36) in position 
space representation by a soft factor. This factor is the vacuum expectation value of a 
product of Wilson lines, with one Wilson line for each external quark or antiquark in the 
multiparton distributions. We thus have 



da 



Y{i=i dxi dxi d^Qi C 



XiXiSj 



i=l 



d'^Zi 



-IZiQi 



d^y 



a\,a2\ mm 



,nn' i^iT 



qq\ mm' ,nn';jj' ,kk' 



,a2'\jj',kk' 



+ {further terms} - 



(3.29) 



where the "further terms" describe the remaining combinations of quarks or antiquarks in 
the two-parton distributions, as discussed in section 2.2.1. The soft factor reads 

['S'(j(j] mm' ,nn' ijj' ,kk' i^i: y) 
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{0\[W{yT + ^ziT;v)W^yT + ^zit; [W{yT - lziT;w)W^yT -\zit]- 

X [W{\z2T-,w)W\\z2T]v)]kn[W{-\z2T-,v)W\-\z2T-,w)\n'k'\^). 



\] nv 



(3.30) 



We notice that Wilson lines W{i; v) T4^t(^; y;) and W{i] w) W^{i] v) are contracted pairwise 
in their color indices. From our derivation we see that this color contraction follows from 
the fact that the hard scatters produce color-singlet particles, so that the color indices of 
annihilating quarks and antiquarks are directly contracted with each other. The "further 
terms" in (3.29) have a soft factor Sqq multiplying -F^i ,02 -^ai ,a2 ^^'^ ^ soft factor Sj mul- 
tiplying the product of interference distributions Iai,a2lai,a2- These factors are defined in 
analogy to (3.30) with an appropriate interchange of arguments and indices for W and . 

In analogy to two-parton distributions, we can represent Sgg in a singlet-octet basis 
for index pairs jj' , kk', etc. 



mm' ,nn' -Jj' ,kk' 

2N 

+ 



1 

iV2 



2N 



Defining the matrix 



^qq ^mm' ^n'n^j'j^kk' + 
Sic fb fb X X 



IS*? - ^ A fa fa 
'-'qq "mm' "n'n ^j'j ''kk' 



■'qq '' mm' n' j' j kk' ~^ j\^2 



880 _ j.b j.b j.a j.a 
'-'qq '^mm''^n'n''j'j'^kk' 



(3.31) 



'qq 



'lie _ 180 

^qq '-'qq 

810 _ 880 _ 

'^qq -Jqq 



(3.32) 



we then have 
da 



Yli=i dxi dxi d^qi C 

^ I -^oi ,02 (^ij -^i) y) 
\ -^ai ,02 (-^ii -^t) y) 




n 

i=l 



d^Zi 
(2^ 



d'y 



Fai ,0,2 {■^i J ^i^y) 



+ {further terms} . (3.33) 



One can of course rewrite the cross section in terms of distributions F{xi, ki,y) or F{xi, ki,r) 
depending on transverse momenta. The result involves a Fourier transformed soft factor 
and is a convolution in transverse- momentum variables. 

The soft factor (3.30) for double Drell-Yan production generalizes the corresponding 
factor appearing in the single Drell-Yan process, which reads 

Sg{z) = 1 {o\[w{^ZT;v)W^{^ZT;w)]mj [W{-lzT;w)wH-^ZT;v)]jm\0) (3.34) 

for the annihilation of a right-moving quark with a left-moving antiquark. The color indices 
are now contracted to an overall singlet, as they are in ^^Sqq. The analog of (3.33) is 



da 



dx dx d^q 



ai{xxs) 



d^z 
(2^ 



e^*^*? fq{x, z) Sq{z) fq{x, z) ^ {further term} , (3.35) 



where the "further term" corresponds to a right-moving antiquark and a left-moving quark. 
For the discussion in subsequent sections we note that at z = the product of Wilson lines 
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in (3.34) reduces to the trace of the unit matrix, so that 5*^(0) = 1. Similarly, one finds 
from (3.30) and (3.31) that 



To close this section let us collect the issues in the soft-gluon sector that need to be 
worked out for a full factorization proof. Some of them we have already mentioned. 

• One needs to show that the exchange of gluons in the Glauber region cancels in 
the cross section. Such a cancellation requires a specific choice of ie prescription in 
the eikonal propagators. For the prescription in (3.21), which corresponds to past- 
pointing Wilson lines, one can show that Glauber gluons do cancel in single Drell-Yan 
production. It is natural to expect that the same prescription is appropriate for the 
double Drell-Yan process, if there is any choice for which Glauber gluons decouple 
in that case. 

• The Ward identity (3.22) for attaching one gluon to a collinear subgraph needs to be 
proven, and it needs to be extended to the case where additional gluons are attached 
to the subgraph. One then needs to show that the attachment of an arbitrary number 
of gluons exponentiates to the Wilson lines in the soft factor (3.30). 

In a model theory with Abelian gluons, a corresponding proof should be a rather sim- 
ple extension of the corresponding arguments for single-parton distributions, which 
can be found in [50, chapter 10.8]. An explicit proof for transverse- momentum de- 
pendent factorization in QCD is still lacking even for single hard scattering, as far as 
we know. 

• It must be shown that explicit time ordering of the gluon operators in the soft factor 
(3.30) can be omitted. It must also be established that one can complement the 
Wilson lines along v and w in the soft factor in such a way that one has an explicitly 
gauge invariant definition. We expect that for both issues it should be possible to 
extend a proof for single Drell-Yan production to the double Drell-Yan process, but 
we are not aware of an explicit proof for the single Drell-Yan case. 

The second and third bullet items are closely connected with the corresponding points for 
collinear gluons, which we discussed at the end of the previous section. 

3.2.3 Towards a factorization formula 

In section 3.2.1 we have seen how collinear gluons give rise to the Wilson line operators 
(3.19) in the matrix elements defining multiparton distributions. However, these Wilson 
line operators contain not only collinear but also soft gluons, which are already taken into 
account in the soft factor (3.30). At the level of graphs, this is refiected in the fact that 
the gluon momentum i in figures 14 and 17a can be either collinear or soft. To prevent 
double counting of soft gluon contributions, the factorization formula for the cross section 
requires appropriate subtractions. 




(3.36) 
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Let us briefly recapitulate how this problem can be solved for single Drell-Yan pro- 
duction. The necessary subtractions can be performed by dividing out vacuum matrix 
elements of the form (3.34). There is a certain freedom of whether to absorb these matrix 
elements into the soft factor or into the parton distributions that appear in the final fac- 
torization formula . The former choice was made in the original work [48] of Collins and 
Soper/^ whereas both [49] and [51] have made the latter choice. Finally, in recent work 
by Collins [50] (see [115] for a brief summary) all matrix elements of the form (3.34) have 
been absorbed into the parton distributions, which gives a factorization formula without an 
explicit soft factor. Whichever choice is made, a consistent formulation requires one to take 
matching ie prescriptions in eikonal lines when treating collinear or soft gluon attachments, 
as we did in (3.3) and (3.21). 

Another detail that admits several choices is the direction of the path in Wilson lines. 
In section 3.2.1 we have seen that the approximations needed for right-moving collinear 
gluons require a vector v that corresponds either to large negative or to central rapidity, 
where we define the rapidity of a spacelike vector as 

(3.37) 

One cannot take the limit yy — )• — oo, i.e. one cannot take v lightlike in the parton density, 
since this would give divergences from the region where gluons coupling to eikonal lines have 
small i'^ but large i~ , i.e. from the region of large negative gluon rapidities [48, 50, 116]. 
The approximations for soft gluons in section 3.2.2 require a vector v with large negative 
rapidity and a vector w with large positive rapidity. Again one cannot take the limit where 
yy —7- — oo and — t- -|-oo, as we shall see explicitly in section 3.3.1. However, this limit 
can be taken for appropriate combinations of matrix elements, which leads to important 
simplifications, see [49] and [50]. 

A further technical point is that the matrix elements discussed so far include contribu- 
tions from self energy graphs of Wilson lines and from graphs where gluons are exchanged 
between different Wilson lines pointing in the same direction (see e.g. figure 24 below). For 
spacelike vectors v and w such graphs give infinite results already at tree level, as shown 
in appendix A of [117]. Such graphs do not appear in the derivation of the factorization 
formula: as we have seen in the two previous sections, Wilson lines appear when treating 
gluon exchange between partons that have a large rapidity difference. The offending graphs 
cancel in the combination of matrix elements that appears in the final factorization for- 
mula, but in the individual factors they must be explicitly excluded. (Only in the scheme 
of [50] do these graphs already cancel in the parton distributions.) 

Finally, the hard-scattering subgraphs have radiative corrections themselves. Since 
we require the produced bosons to have small transverse momenta, there are only virtual 
corrections: radiation into the final state can only be collinear or soft and is included in 
the collinear or soft factors. For Drell-Yan production at one-loop accuracy, one thus only 
has the vertex correction to the quark-antiquark-photon three-point function. The regions 
of soft and collinear gluon momenta in the virtual graphs have to be explicitly subtracted 

^^This may not be quite obvious in [48] but has been clearly pointed out in section X.A of [49]. 
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in the definition of tlie liard-scattering cross section, in order to ensure tliat this factor is 
dominated by large virtuahties. This removes in particular the well-known soft divergences 
of the virtual graphs (which in the more familiar case of inclusive observables cancel when 
real emission graphs are added). 

We expect that the above procedure can be generalized to the case of double Drell-Yan 
production. The division by vacuum expectation values of the form (3.34) will be replaced 
by multiplication with the inverse of the matrix (3.30) in color space. We leave it to future 
work to show that this can actually be done. In the remainder of this section, we will take 
a closer look at the elementary building blocks of factorization, namely at the soft factor 
in (3.30) and at the dependence of the proton matrix elements of the operator (3.19) on 
the direction v of the Wilson lines. 

To conclude this section we note that a factorization theorem for the double Drell- 
Yan process also needs to provide a proper separation between the production of the two 
gauge bosons by one or two hard-scattering processes. We will discuss this problem in 
section 5.2.3. 

3.3 The soft factor at small transverse distances 

If all three transverse distances y,zi,Z2 are small compared with a hadronic scale A~^, 
the soft factor in (3.30) is dominated by perturbative dynamics and can be evaluated in 
perturbation theory. In this section we compute the short-distance form of this factor to 
leading order in the strong coupling. 

From (3.11) and (3.12) one obtains the representation 

ig d\vA\i-\v)e = -igev^ j ^ e^?^ j A e*^'^S(C) (3.38) 

for the exponent of the Wilson line W{^; v). Together with the Feynman rules for eikonal 
lines and their coupling to gluons in figure 16 this gives a Feynman rule for the order g 
term of the Wilson line to the left of the final-state cut, as shown in figure 22. 

3.3.1 The basic graphs 

The expansion at 0{as) of the soft factor (3.30) involves the three types of graphs shown in 
figure 23, which we now calculate. Graphs a and b already appear in the soft factor (3.34) 
for the single Drell-Yan process, whereas graph c is specific for multiparton interactions. 
As discussed in the previous section, we discard graphs as in figure 24, where gluons are 
exchanged between Wilson lines pointing both along v or both along w. 

To regulate ultraviolet divergences, we work in 4— 2e dimensions, and to exhibit infrared 
divergences in individual graphs we use a small gluon mass A. For the time being we omit 
the color matrices in the Feynman rules. One finds that each pair of graphs in figure 23 
gives the same result, so that in the following we only list the expressions of the left-hand 
graphs and multiply by two. 
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Figure 22. Feynman rules for the term of order g in the expansion of the Wilson lines defined in 
(3.13). The rules for lines to the left of the final-state cut (indicated by the dashed line) follow from 
(3.38), and those for lines to the right of the final-state cut are obtained by complex conjugation. 
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Figure 23. Basic graphs contributing to the soft factor at order as 



Before renormalization, the vertex correction graphs in figure 23a give 



(27r)4-2e -^w + ie'^ 'i^-X^ + ie' ' iv + ie 

vw 2. f f°° di- 1 



V+W+ ^ J (27r)2-2e J^^ 2tt 2e+i- -f-X'^ + ie 

^ ^ (3.39) 



-64- 



Figure 24. Graphs with gluon exchange between eikonal lines having the same rapidity. These 
graphs are excluded in the evaluation of soft factors, as discussed in the text. 



where we recall that the vectors v and w have zero transverse momenta and satisfy v~ > 0, 
< 0, > 0, w~ < 0. For < all poles in i~ are on the same side of the real axis, 
so that this region gives a zero contribution to the integral over i~ . For > we close 
the integration contour in i~ around the pole of the gluon propagator and obtain 

Ua = -as / ,^ ,o / dii+y 



V w J {2tt) 



(3.40) 



As both poles in {i'^)'^ are on the same side of the real axis, one can deform the integration 
contour to — oo < (£^)^ < and obtains 



via 

Ua = -2as , log 



v-w+-v+v~ ^\v+w-J'^ J (27r)2-2^ £2 _^ A2 



v~ 






v 




v~ 











cp-^'£ 1 



= + logf V^") r(.) (^Y. (3.41) 

We see that both vectors v and lu must be chosen away from the light cone, since taking 
even one of them lightlike gives an infinite result for Ua- With the definition (3.37) for the 
rapidity of a spacelike vector, we have 

, = 2(y^ - Vv) , , = tanh(y^ - y^) . (3.42) 

V^W J V — v~^v 

We require a large rapidity difference, \yw — Uvl ^ ^, i-e. v~'w~^ » v~^w^ . This is satisfied 
by the choice yt, <C 1 and 3> 1 made in section 3.2.2, but it also allows one of v or w to 
have central rapidity, as long as the rapidity of the other one is large (positive or negative) . 
We can then cLpproximate ta,nh.(?/|y; — yy) ^ 1. The expression (3.41) is ultraviolet divergent, 
and using MS subtraction we obtain the renormalized result 

Sai^t) = [C/a(^*)]"° = (y«, - y.) log ^ . (3.43) 



- 65 - 



For the graphs in figure 23b, which describe gluon emission into the final state, we have 
= 2/^'^ / T^TO (^i^^) ^(^' - A') e^^T^-^^'r' Hgw^ {-g^p) {tgv^-^ ^— 



26 /) /* oo /* oo 



1 1 



vw 



j2-2e p , /-oo 



v-w-"^ J (27r)2-2^ 

1 



X , , ^ r-r^-^^ • (3-44) 



Comparing with (3.40) we observe that Ua + Ub = at ^ = ^' . For $ ^' the integral over 
i in (3.44) is ultraviolet finite and no subtraction is needed before we set e = 0. We thus 
obtain 

Sbiz) =^ [Ub{z)Y=' = 4a. {y^ - y.) | ^ ^"'^2^ = V ^^"^ " ^''^ ^MM^l) 

a A2z2 

= iVw - Vv) log , (3.45) 

A->o vr Oq 

where 6o = 2e~"' and 7 is the Euler number. 

The graphs in figure 23a and b give the full 0(0.) contribution for the vacuum matrix 
element Sq defined in (3.34), which appears in single Drell-Yan production. Let us evaluate 
this contribution as a side result. The color factor for all graphs is A^^^ tr(t"i'^) = Cp in 
this case, so that we have Sg{z,fi) = 1 + S{z,fj.) + 0(02) ^[^^ 

S{z, /i) = Cf [Saifi) + Sbiz)] = -Crivn^- Vv) log ^ (3.46) 

for z 7^ 0. Notice that the infrared divergences regulated by a gluon mass A have cancelled 
in the sum over all graphs, as is required for a perturbative evaluation of Sq. For z = 
the relation Ua + = ensures that the condition 5^(0) = 1 does not receive radiative 
corrections, in agreement with the general result discussed below (3.34). This complete 
cancellation between real and virtual corrections plays a crucial role for collinear factor- 
ization, see section 3.5. The fact that the limit z — )■ of (3.46) is infinite rather than zero 
is due to our use of modified minimal subtraction. For z = the ultraviolet divergences 
in Ua and Ub cancel each other (as do the finite parts of the graphs), so that there is no 
ultraviolet subtraction for their sum. At z 7^ 0, however, Ua + Ub does require ultraviolet 
subtraction since the first term is ultraviolet divergent whereas the second term is not. 
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Finally, the graphs in figure 23c give 
Uc = 2/.- / (-5-") ~''f, i-^9v^) ' 



^ J (27r)2-2e 27r 2i+i- - ^ - \^ + ie 

^ ^ (3.47) 



+ ^1+ -ie i- + 



For < all poles in i~ are on the same side of the real axis and give a zero net result, 
whereas for > the integral over i~ can be obtained from the residue of the gluon 
propagator. This gives the same expression as in the second step of (3.44). The graphs in 
figure 23c hence give the same result as those in figure 23b: 

S,{z) = Sb{z) . (3.48) 

Note that the graphs in figure 23b and c only differ by the position of the final-state 
cut. The result (3.48) shows that this does not matter and provides a consistency check 
for the omission of explicit time ordering in (3.27). Namely, time ordering of the two 
fields ^(0 and A{^) associated with the gluon propagator is required when the gluon does 
not cross the final-state cut. When it does cross the cut, no time ordering prescription 
arises, however, since the two fields then have a complete set of final states between them: 
Ex Mm) {X\A{0 = A{OA{0. 

The identity (3.48) also holds before setting e = and can then be used also at z = 0. 
This provides a simple explanation of the relation Uh{0) = —Ua discussed above. Indeed, 
the equality Uc{0) = —Ua is already evident from the starting expressions (3.39) and 
(3.47), which only differ by the sign of the eikonal propagator ^i/{£v + ie). This is a direct 
consequence of the Feynman rules in figures 16, since the upper eikonal lines to which the 
gluon couples in figures 23a and c only differ by the flow of color charge. 

For later use we note that the graphs in figures 23b and c change sign when the flow 
of color charge in one of their eikonal lines is reversed, whereas those in figure 23a remain 
the same. 



3.3.2 Soft factor associated v^ith two-quark distributions 

With these building blocks at hand we can construct the soft factor associated with two- 
parton distributions in the cross section formula (3.33). We first discuss the factor Sqq in 
some detail and give the results for Sqq and 5/ in the end. 

The graphs contributing to Sqq are those in figure 25. In the graphs of figure 25a and 
bi there is either no gluon coupling to the pair {14} or no gluon coupling to the pair {23} 
of eikonal lines, so that the lines have the same color coupling at the bottom and the top of 
the graph. The color factor for the vertex correction graphs 25a \s Cp, independent of the 
color flow along the eikonal lines and of the way in which their color indices are coupled. 



-67- 



For figure 25bi we have color factors 



■^trif^n trl = Cf 



for ^^Sgg , 



for , 



(3.49) 



where the prefactors and 4/(A^^ — 1) come from the anafog of the color decomposition 

(3.31) for Sqq. We note that these color factors are the same as for the ladder graphs we 
shall discuss in section 5.1.3. 

By contrast, the graphs in figure 25b2 and c can change the color coupling of the eikonal 
lines. As an illustration, consider the first graph in figure 25b2 with the color indices at 

jj' kk' 



the bottom multiplied by octet matrices t''-,t'l,,. The color structure is then 
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2N 



"mm' nn' 



(3.50) 



and we see that this graph contributes to both ^^S and ^^S. If the indices jj' and kk' are 
coupled to singlets, then the graph is proportional to tj^^/ and thus contributes to ^^S. 
For the first graph in figure 25c the color factors are analogous, with the difference that 
one has 



/ mm' \ / nn' o 



nn' 2 

iV2 



0mm' Onn' + 2^ 
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^mm' ^nn' 
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(3.51) 



instead of (3.50). Putting everything together and using the results of the previous section, 
we have 
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(3.52) 
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Figure 25. Graphs for the soft factor associated with two-quark distributions. The "mirror graphs" 
not shown are as in figure 23. The numbers, position arguments and color labels on the eikonal 
lines correspond to those of the quark lines in the distribution Fa^^a2i see figure 5. 



The terms Sc come with an overall minus sign, because in graph 25c one of the eikonal lines 
coupled to a gluon has its color flow reversed compared with graph 23c, which reverses the 
sign in the eikonal propagator according to the Feynman rules in figure 16. 

The matrix elements in (3.52) are such that we can replace Sa and = Sc by their 
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Figure 26. A higher-order grapli for the soft factor, which connects more than two eikonal lines. 



sum, which is infrared finite according to (3.46). Using S = CpiSa + Sb) we have 

Sqq{Zi, y, fl)=[l + S{ZI,H) + S{Z2, /U)] 



1 
1 



cSd 
cSd -{l + c^)Sy-2c'Sd. 



with a color factor 



and linear combinations 



1 1 



(3.54) 



Sdizi,y) = s(y+ ^^-^^,/^^ +s(y- ^^-^^,fi 

( — Z2 \ ( Z\ — Z2 

-S\y^ — 2 — ' -^\y 2 — ' ^ 

Sy{zi,y) = S{zi,fi) + S{z2,n) - s(^y + ^^^ '^) ~ ^{v ~ ' ^'^'^^^ 

Note that the ^ dependence of S has canceled in Sd and Sy. 

At one-loop level the soft factor for two-quark distributions can thus be expressed in 
terms of its analog S for single-parton distributions. This simplification will most likely no 
longer hold at higher orders in a^, since one then has graphs like in figure 26, where more 
than two eikonal lines are connected. 

It is easy to see that the soft factor Sqq{zi,y) for quark-antiquark distributions is 
obtained from Sqq{zi,y) by replacing Z2 — )• —Z2- This corresponds to the interchange of 
the labels 2 and 3 in the definitions of -Fai,a2 and -Fq^^qj (see (2.85) to (2.87)) and uses the 
result that the kernels Sb and Sc for gluons that cross or do not cross the final-state cut 
are identical. 

The soft factor multiplying the interference distributions is obtained from the same 
graphs as those in figure 25, apart from changes in the color labels and in the color fiow 
of the eikonal lines. Making the appropriate replacements of position arguments in Sqq we 
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get 



11 



Si = l + Cf 



2S.{^) + SAy+'-^)+sJy-'-^ 



8% = l + 2Ci.5,(A.)- — 



Sh[y^ 7, — ]+Sh[y — 



+ [Cf 



2N 



Sciy^ — +Sc(y — 



1 



+ - [S,{z,) + S,iz2)] , 



2N 



which can be rewritten as 



Sc(y + ^^iT^) +Sc(y- ) - Stizi) - Skiz2) 



(3.56) 



Si{zi,Z2,y) = [1 + S{zi,fi) + S{z2,fi)] 




(3.57) 



3.4 CoUins-Soper equation 



As is well-known, cross sections with measured transverse momenta \q^\ ~ (/t much smaller 
than the hardest scale Q in the process contain large logarithms in qr/Q, which need to 
be summed to all orders in as in order to have a perturbatively stable result. A powerful 
method to resum these Sudakov logarithms is due to Collins, Soper and Sterman (CSS) 
[118]. This method uses transverse-momentum-dependent factorization and is therefore, 
up to now, limited to the production of color singlet particles, such as a Drell-Yan lepton 
pair or a Higgs boson. In this section we show how this formalism extends to double Drell- 
Yan production, and we calculate the corresponding Sudakov factor to next-to-leading 
logarithmic accuracy. 

We begin by a brief account of the CSS method for processes with a single hard 
scattering initiated by quark-antiquark annihilation. As we mentioned in section 3.2.3, 
the Wilson lines in the definition of transverse-momentum dependent parton distributions 
must be taken along a direction v with a finite rapidity. Their dependence on this rapidity 
is governed by the Collins-Soper (CS) equation [48], whose solution resums Sudakov log- 
arithms to all orders. By Lorentz invariance, the distributions depend on v via the scalar 
parameter 



12 



(3.58) 



To avoid the appearance of square roots, our definition (3.58) follows the convention of Ji et al. [51] 



and differs from the one of Collins and Soper [48], with ^ 



C[48] ■ 
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Figure 27. Graphs describing tlie dependence of a single parton distribution on the Wilson line 
rapidity at order a^. The complex conjugate graphs are not shown. Note that graph b includes 
the case where the gluon couples to the left quark line, i.e. graph a is included in b. 

As discussed earlier, we require v to be spacelike (although the manipulations in the fol- 
lowing would also work for timelike v, which is the choice adopted in [51]). In terms of the 
rapidities of v and yp of the proton, we have 

C = 2M sinh(yp - y,) ^ M 6^'^"^'" , — ^ = - tanh(yp - y,) A ~ (3.59) 

91ogC dy^ dyy 

for a right-moving proton, where M is the proton mass. The approximations are valid for 
Uv Up, which we have seen to be necessary in section 3.2.1. This implies ^ M. 

Following the original paper [48] we work here with "unsubtracted parton distribu- 
tions" in the parlance of [51, 119] and [50], i.e. with matrix elements of operators con- 
structed from quark fields and Wilson lines along v as in (3.19). Subtractions for the soft 
momentum region are not made in the parton distributions but in the soft factor that 
appears in the cross section (see [49]). The fields in the operator are renormalized, so 
that the distributions depend on an ultraviolet renormalization scale fi. Fourier transform- 
ing from transverse momentum to position space, we then have single-quark distributions 
f{x,z;(^,fi). The fi dependence is given by a homogeneous renormalization group equa- 
tion [48] 

^ fix,z;C,fi) = 2-fg{asin)) f{x,z;C,IJ-), (3.60) 



dlog fi 



where 7^ = SCpcts/i'i^'^) + C(ck^) can be identified with the anomalous dimension of the 
quark field in the axial gauge vA = (where the Wilson lines reduce to unity). In a 
covariant gauge it corresponds to renormalization of the composite operator W{^)q{^) or 
g-(e)VFt(e), see (3.19). 

In terms of graphs, the dependence of / on v arises from the propagators of eikonal 
lines and from their coupling to gluons, as is obvious from the Feynman rules in figure 16. 
A power counting analysis shows that dominant contributions to df /d( come from regions 
where the momentum i flowing through eikonal lines is either hard or soft; contributions 
with i collinear to the proton are power suppressed. The only important graphs where ^ is a 
hard momentum have the form of a vertex correction shown in figure 27a (at higher orders 
this vertex graph becomes dressed with further gluon and quark lines) . If the momentum i 
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flows through spectator partons (see figure 27b and c), the hard region is power suppressed 
and only the region of soft i is important. 

The contribution from the region of hard momenta i reads 

—^—f{x,z;Q,ii) = G{xC/fJ-) f{x,z;C,n) . (3.61) 

O log Q hard 

Since £ is hard, the kernel G can be calculated perturbatively. G is extracted from sub- 
graphs whose external lines are the quark on the left of the final-state cut and the adjacent 
eikonal line, or the corresponding two lines on the right of the cut. As a result, G is inde- 
pendent of the transverse distance z between the two quark field operators. Furthermore, 
it depends only on the longitudinal quark momentum xp rather than on the proton mo- 
mentum p, so that the dependence on ^ is via the combination = 2(xp) v / ^/\v^ . For 
dimensional reasons this parameter must be divided by fi, since G is dominated by hard 
momenta and hence independent of nonperturbative scales. At leading order in a^, one 
obtains G from the graph of figure 27a and the complex conjugate graph, with subtractions 
made for the region of soft £. Using MS subtraction for the ultraviolet divergences, one 
finds [48, 119] 



G{xC/fi) = -^GF 



+ O(oj), (3.62) 



If the momentum carried by the eikonal lines and hence by the gluons they couple 
to is soft, we have graphs with soft gluons coupling to parton lines that move fast to the 
right. We can then use the same procedure as in section 3.2.2, i.e. approximate the gluon 
coupling to the right-moving particles and then use a Ward identity. The gluons then 
couple to eikonal lines with large positive rapidity, associated with a vector w, with one 
eikonal line for each parton line attached to the collinear subgraph. This gives 

f{x,z;C,fi) = _ '^'^'? (^'/^) __}_ f{x,z;C,fJ') , (3.63) 

dlogC soft dy^, bg' {z,ii) 

where we used 9/(5 log C) = —d/{dyy) and where we have explicitly displayed the depen- 
dence of Sq defined in (3.34) on the two vectors v and w. As shown in [50], one can take 
the limit of lightlike w in (3.63) and thus has 

/(x, z; C, /u) = [G(xC//i) + K{z,p.)\ /{x, z; C, fi) (3.64) 

with 



K{z,t,) = - hm (3.65) 

Having taken the limit of infinite the dependence of K on y.y has disappeared as well, 
since by Lorentz invar iance K could only depend on — Vv For large z the kernel K 
is dominated by nonperturbative dynamics, just as Sq. For sufficiently small z we can 
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however use the perturbative expression of Sq derived in the previous section. From (3.46) 
we readily obtain 

K{z,f,) = -^CFlog^, (3.66) 

verifying that at lowest order K is independent of v and w. 

An analysis of the ultraviolet divergences in G and K shows that they satisfy renor- 
malization group equations [48] 

,/ G{xC/fi,as{n)) =--pr — K{z,^) = -iK{as{lj)) , (3.67) 
dlogn ' alog/i 

so that the sum G + is independent of /i. At lowest order in this is readily verified 
from (3.62) and (3.66), which also give the leading term 

lK{as) = —CF + 0{al) (3.68) 
vr 

of the anomalous dimension. The 0{a1) term is also known [118]. 

We are now ready to generalize the CS equation to (unsubtracted) two-quark distri- 
butions. The contribution from hard eikonal momenta £ is again given by vertex graphs 
as in figure 27a, with one graph for each of the four quark legs. This gives just the sum of 
kernels G{xiQ/ ^i) + G{x2C/ ^j)- As for the soft momentum region, the argument leading to 
(3.63) can be repeated. In section 3.2.2 we have seen that the Ward identity for soft gluons 
coupling to collinear lines in a two-parton distribution has a nontrivial color structure. We 
therefore now have a matrix equation in color space. 



d I -^ai,a2 (-^jj -^i) yi C) /^) 



91ogC yFa^^a^ixi.z^y-Cli)^ 
with a kernel 



soft 



K,,(..,y;^) {'P^.^^i-^^-^^y-^C,,)\ (3gg^ 



K,,(z„y;/.) = - hm ^?I!ifll^[S^^-(z,,y;/.)]-\ (3.70) 

Analogous equations hold for quark-antiquark and for interference distributions, with ker- 
nels K constructed from the appropriate soft factors S. Note that the kernels are sensitive 
to the color charge of partons (i.e. to the difference between quarks and antiquarks) but 
not to their polarization. For K this follows from the analogous property of S, whereas 
for G it follows from parity invariance applied to the relevant subgraph with one external 
quark and one eikonal line. Putting everything together, we can write the CS equation for 
a two-quark distribution into the form 

^ (1^) = [G{x,c,^^) + G{x2C^^) + K{z,,^,) + K{z,,^,)] (^f^ 

+ M{zuZ2,y)Ly\ (3.71) 
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with 



M(zi,Z2,y) = K{zi,Z2,y;n) - [K{zi,fi) + K{z2,id)] 




(3.72) 



According to (3.67) the ^ dependence cancels between the kernels G and K in (3.71). The 
matrix M is independent of /i as well. This can be traced back to the fact that the only 
ultraviolet divergent graphs for S (and hence for K) have the form of a vertex correction as 
in figure 25a. As discussed in the previous section, these graphs only contribute to the color 
diagonal elements ^^S and ^^S and are the same in both cases, since they are insensitive to 
the color of the different eikonal lines. 

Finally, the dependence of -Fa^^aj, -Fai,a2 etc. on the renormalization scale is given by 



d 



dlogjj, 



F{xi,Zi,y;C,f^) = 47g(a5(/x)) F{xi, Zi,y;C, fJ-) 



(3.73) 



for both and ^F. This is because ultraviolet renormalization in F is performed for 
individual operators W{^)q{^) and q{C)W^^). 

3.4.1 General solution 

Before solving (3.71) let us first consider the simpler equation 
d 



dlogC 



F{xi, Zi,y; C, /i) = G{xiC/fi, as{n)) + G(x2C//^, as{fi)) 



+ 2Ki2izi,Z2,fi) F{xi,Zi,y;C,i^i) 



where we have abbreviated 



1 r. 



Ki2{zi,Z2,fi) = - \ K{zi,ij) + K{z2,fi) 



(3.74) 



(3.75) 



We have included the argument of the running coupling in G since this will be needed 
shortly. For the moment we do not assume that K is given by a perturbative expansion. 
The solution of (3.74) can be obtained by adapting the well-known solution of the CS 
equation for single-parton distributions [48, 118, 120]. We have 

F{xi,Zi,y;C,f^) = exp[-S{xiC, zi, Z2, Ho) - 5(2:2 C, ^i, ^2, Mo)] F^"^ {xi, Zi,y; fi) , (3.76) 

where F^" specifies the initial condition of evolution in The scale fiQ should be chosen 
such that F'^o does not depend on widely disparate scales. If this is not possible because 
Zi and y widely differ in size, further steps may be required in order to resum all large 
logarithms. We note that since 5 is /x independent (see below), the ^ dependence of F^" 
is given by the same renormalization group equation as in (3.73). The Sudakov exponent 
in (3.76) reads 



S{xC,Zi,Z2,Ho) 



dC_ 



G(xC7/",as(/^)) + Kl2{zl,Z2,^l) 



(3.77) 
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with X equal to xi or X2- It is well-known from the solution 



f{x,z;C,fJ-) = exp[-S{xC,z,z;no)] f''°{x,z;ti) 



(3.78) 



of the CS equation (3.64) for single-quark distributions. We note that a more general set 
of solutions can be obtained by multiplying x(^ with a constant C2 in the upper integration 
limit of (3.77); the initial condition in (3.76) then depends on that constant. One can 
easily restore the C2 dependence of the expressions to follow, but for simplicity we limit 
ourselves to the choice C2 = 1 here. 

The integrand of (3.77) contains functions that depend on the ratio of two large scales. 
To make this dependence explicit, one uses the renormalization group equation (3.67) for 
G and K. Obviously, K12 has the same fi dependence as K, so that the /i dependence 
cancels between G and K12 in (3.74) and (3.77). Using (3.67) one can rewrite 



G{xC /fi,as{fj.)) + Ki2{zi,Z2,n) 

= G{l,asixC')) + Ki2{zi,Z2,Ho) - / 

J u 



(3.79) 



Inserting this into (3.77) one can perform the integration over for the terms containing 
or K12 and obtains 



S{xC,Zi,Z2,fio) = — 



JK{as{^J')) log G(l,Q^(/i)) 

A* 



Ki2{zi, Z2, Ho) log 



xC 



(3.80) 



The term with 7^- in (3.80) gives rise to the leading double logarithm in xC/fio, whereas 
the other terms give only single logarithms. Using (3.62) and (3.68) and neglecting the 
running of Os one has 



S{xC,Zi,Z2,fJ.o) = — Gplog^ — 



a 



— Gp + Ki2{zi,Z2,IJ.o) 



vr 



xC 



log ^ + 0{ai) (3.81) 



at leading order in a^. A more precise expression is obtained by rewriting J dfi/fi = 
das/ p{as), where (3 = das{^i) / d\og fj? . After expanding \/j3{as) in a^, the integral 

in (3.80) is straightforward to evaluate for the one-loop expression (3.62) of G and the 

two-loop expression of 7x (given e.g. in [118]). 

The form (3.80) is valid even if Ki2{zi, Z2, h) cannot be evaluated perturbatively 

because one or both of zi and Z2 are large. If both distances are small, K{zi, /i) and thus 

Ki2izi, Z2, IJ^) is given by a power series in as(/i). An alternative form of the Sudakov 

exponent [118] is then obtained by rewriting (3.79) as 



G{xC' /^i,asifJ.)) + Ki2{zi,Z2,n,a{n)) 

= G(l,Q,(xC')) +i^i2 (21,^2, /io,as«')) - / ^A{zl,Z2,^lo,as{^Ji')) (3.82) 

J an fJ- 
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with 



d 



A{zi,Z2,iiQ,as) = 7i^(as) + 2/3(as) ^— Ki2(2i, ^2, /^o, a^) , 



(3.83) 



where we now distinguish between the exphcit ^ dependence of K12 and the imphcit de- 
pendence via the running couphng. One then obtains 



S{xC,Zi,Z2,Hq) = I 
J u 



A{zi,Z2,i^o,as{^i)) log G{l,asifJ,)) 



Ki2{zi,Z2,Ho,asin)) 



(3.84) 



where all perturbative functions are evaluated with as at the same scale. With the one-loop 
expression of K in (3.66) we have 



Ki2{zi,Z2,fJ.,as 



Cirlog -2 h ©(a J 



A natural choice for the starting scale of evolution in is thus 



Mo 



^1 ^2 



(3.85) 



(3.86) 



with a constant Ci of order 1. If one takes Ci = 6q then K12 vanishes. 

It is now easy to write down the solution of the full CS equation (3.71) for a two-parton 
distribution. It is given by 

(^F{xi,Zi,yX,fJ-)\ r , ,-, 

[s-p/^ _ ^ ,,^] = ^^V[-S{XiC,Zi,Z2,Ho) - S{X2C„Zi,Z2,Ho)\ 



^F{xi,Zi,y;C,l^), 



X exp 



M{zi,Z2,y) log 



VxiX2 C 

Mo 



'^F>'^'{xi,Zi,y;fi)' 
^^F''"{x^,Zi,y■,n)^ 



(3.87) 



with 5 given by (3.80) for arbitrary values of Zi and by (3.84) if both zi and Z2 are 
small. The logarithm in the second line has been chosen such that it coincides with the 
one that multiplies 2Ki2 when one evaluates —S{xi() —S{x2C) from (3.80). Other choices 
are possible and lead to different initial conditions ^F^^° and ^F^o. 

Unless all distances zi, Z2 and y are small, the matrix M cannot be calculated pertur- 
batively and we cannot further simplify the exponentiated matrix in (3.87). Nevertheless, 
(3.87) contains some important information, since it gives the explicit form of the de- 
pendence on the large scales xi( and X2C- In particular, we see that to leading double 
logarithmic accuracy, where only squared logarithms of xiC/mo and X2C/fJ'0 are retained, 
the Sudakov factor for two-quark distributions is the same for ^F and for ^F and given by 
the product of the corresponding Sudakov factors for single-quark densities with momen- 
tum fractions xi and X2. At next-to- leading logarithmic accuracy, ^F and ^F mix under 
evolution in with the amount of mixing depending on the transverse distances zi, Z2 
and y. 



-77- 



It should be possible to generalize the CS equation (3.71) and its general solution (3.87) 
to the case of multiparton distributions for more than two partons. The same holds for 
multi-gluon distributions, where the general structure will remain the same but the kernels 
G and K will be different. 

3.4.2 Small transverse distances 

Let us now consider the situation when zi, Z2 and y are all small enough to calculate M 
in perturbation theory. The kernels K{zi,fi) in (3.71) are then given by (3.66) at leading 
order in a^, and the Sudakov exponent S in the solution (3.87) of the CS equation can 
be evaluated from (3.84). It remains to investigate the matrix in (3.87), where we 
abbreviate 



L = log^^. (3.88) 

We treat the kernel Mgg for two-quark distributions -Fai.a2 i^i detail and discuss its analogs 
for quark-antiquark distributions Fa^ a^ and interference distributions Iai,a2 later. 

Using the definitions (3.70), (3.72) and our perturbative result (3.53) for Sqq, we readily 

find 

with c given in (3.54) and 

Z1 + Z2 \ Z1+Z2 



Kd{z,,y) = K(y + ^^^,f?j +K(y 



2 



/ Zi — Z2 \ ( Z\ — Z2 
-K[y + - K{y ,fJ. 

Ky{zi,y) = K{zi,f,)+K{z2,fi) - K(y + ^^^,/") - K(y - ^^^,/^) (3-90) 

in analogy to (3.55). Using the explicit form (3.66) or the renormalization group equation 
(3.67) for K, we see that M is /i independent, as we anticipated earlier. We must of course 
choose a scale in when using the one-loop result (3.66) for evaluating and Ky, which 
gives rise to a residual scale dependence of order a^. The situation is the same as for a 
physical (and hence formally ^ independent) quantity evaluated in fixed-order perturbation 
theory. An appropriate scale of in and Ky will be constructed from zi and y. 

Let Dqg be the diagonal matrix with the eigenvalues of M^^ and let Egg be the matrix 
whose columns are the corresponding eigenvectors, i.e. 




MgqEgg = EggDq, with = ( ^ ) . (3.91) 

One then has 
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The matrix (3.89) has eigenvalues 



d± 



(3.93) 



and a matrix of eigenvectors 




(3.94) 



so that 



1 



d+ — d 

Let us see how this matrix behaves for 

\cKd\^ \Ky\ . 

One can then Taylor expand the square root in (3.93) and obtains 



d.e^'^- 



d. 



1 + c2 K, 



+ c 



Hi. 

y 



-{l + ^)Ky + 0{^Kd 



(3.95) 



(3.96) 



(3.97) 



if Ky > 0, whereas the role of (i+ and in (3.97) is interchanged if Ky < 0. In both cases 
one gets 



exp 



L 



iV2 K„ 



( 

J_Kd 
\Nb Ky 



Nb K. 



) 



\Nb Ky) J 



(3.98) 



where we have traded the color factor c = l/V N"^ — 1 for 

= : . 



(3.99) 



Since c ~ the condition (3.96) holds in the large- limit. Inserting (3.98) into 

(3.87), we see that ^-F(C) is then controlled by the initial condition ^F^^ because the ad- 
mixture from ^F^" is suppressed, although only by Whether ^F{Q) is dominated by 
i^A»o Qj. depends on whether the suppressed factor in the lower row of (3.98) 
or the exponential e~^^^^y is smaller. In either case ®-F(C) is parametrically smaller than 
^F{(^). To which extent the large- limit gives a valid description of the physics for = 3 
depends on the relative size of Ky and K^i, as well as the relative size of ^F/*" and ^F^o. 
This can only be decided by a more detailed analysis, which we will not attempt here. 

There is, however, a region of phase space where (3.96) holds beyond the large- limit. 
From its definition (3.90) we see that K^ vanishes if z\ = Z2 = 0. In the double-scattering 
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process one has \zi\ ~ \z2\ ~ ^/qt, so that the hmit |2i|,|22| ^ \y\ is relevant in the 
region where \y\ is much larger than I/qt- Taylor expansion then gives 



vr 



2(l/2i)(^/Z2) Z1Z2 



(^2)2 

^ 2Cf log 



This implies 



< 



^1 ^2 



log' 



+ 



y 



y^ 



Zi 



y^ 



+ 



\y\ 



+ 



\Zi\ 



\y\ 



(3.100) 



(3.101) 



y V i^iii^2i 

for the factor in the off-diagonal elements of the matrix (3.98), whereas for the exponential 
factor we find 



\Z2\ 

y2 



Np 



with p 



vr 



as . \zi\\z2\xiX2C'^ 
^=2^^°^ C[ 



(3.102) 



where we have chosen ^uq as in (3.86). We thus find that the octet admixture in the ^ 
evolution of ^-F(C) is power suppressed by \zi\\z2\/y'^, and that ^-F(C) is power suppressed 
compared with ^F{Q by 



^Faj^^a2{Xi,Zi,yX,fJ^) ^Fai,a2{Xi,Zi,y;^) (\zi\\z2 



rmii(l, Np) 



(3.103) 



"'^-^ai,a2 (^i) -^ii Z/i C? /^) "'^-^oi ,02 (^i j -^^i > l/j /^) V Z/^ 

In straightforward generalization of single Drell-Yan production [48, 50], an adequate choice 
of Q in the cross section for double hard scattering is xiX2C^ ~ Q^- Together with \zi\ \z2\ ~ 
l/g|, this gives p ~ (os/vr) log((5/(/T)- Again, a more quantitative picture can only be 
obtained by a detailed analysis. 

The preceding results all rely on the validity of perturbation theory for the soft CS 
kernel and thus require not only zi and Z2 but also y to be perturbatively small. We 
cannot draw any strict conclusions about the case where zi and 22 are small, whereas y is 
in the nonperturbative region. However, we observe that the power suppression parameter 
\zi\ \z2\jy^ becomes smaller rather than larger in this case. One may thus speculate that 
the general features of our analysis, namely the autonomous C, evolution of and the 
suppression of '^F will continue in the nonperturbative regime. 

We conclude this section by noting that the CS equation and its solution for quark- 
antiquark distributions -^01,02 is readily obtained from the previous results by replacing 
Z2 — >■ —Z2 and that corresponding replacements are to be made for Fa^^a^ ™d -Fai,a2- 
This follows from the corresponding property of the soft factor Sgg discussed at the end of 
section 3.3.2. 



Interference distributions. The Collins-Soper equation for interference distributions 
^Iai,a2 and ^Iai,a2 has the same form as (3.71) with F replaced by I. The appropriate kernel 
M/ in the perturbative regime follows from S/ in (3.57) and reads 

/ Ky c{Ky+Kd) \ 



M.i{zi,Z2,y) 



(3.104) 
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with and Ky given in (3.90). This matrix has eigenvalues 



-(1 - c") {2Ky + Kd) ± Vl + c^ a/c2 {2Ky + Kd)^ + 



(3.105) 



and exponentiates to 



1 



-cirrf(e^'^+-e"- 



d' - d'_ 



') 



d'.e"^"'- 



Kyie 



Ld' 



d'-d'_ 



'1 c ' 
.c -1, 



For \K(i\ <C i.e. if \zi\, \z2\ <^ \y\, we can Taylor expand the eigenvalues as 



d' 



N 



The result simplifies if we use the orthogonal matrix 



d' 



N + 2 



Ky + 0{Kd 



U 



1 /\/]v^ -VivTi 
^^/ivTl 



(3.106) 



(3.107) 



(3.108) 



that implements the basis transformation from H to the combinations introduced 
in (2.117). We then have 



Ue 



LM.I 



exp 



iV- 1 ^ 



/ 

\AcKy 



Kd 
AcK, 



-LKy/Cp 



-LKy/Cp 



) 



^-LKy/Cf 



Repeating the argument that led to (3.103) we obtain 



/ 

(3.109) 
(3.110) 

^F{xi,z,,y;C,fi) ^F^'o{xi, Zi,y- fi) \ y'^ ) ' ^^"^^"^^ 

For |2;i|,|2;2| <C \y\ the sextet combination of I is hence suppressed compared with the 
antitriplet one, which in its turn is small compared with the singlet combination ^F. 

We finally note that, in contrast to the case of F, the limit \Kd\ <C Ky just discussed 
does not give the same result as the large- limit. In the latter one finds 



with p given in (3.102). Likewise, comparing (3.109) with (3.98) we find 

- - (iV+l)(Af-2) 

^I{xi,Zi,y;C,i_i) ^I^o{xi,Zi,y;n) f\zi\\z2\'^ ^ ^ 



-LKy/Cp 



) 



-LKy/Cp) 



-LKy/Cp 



(3.112) 

with relative corrections of order Ky/{NKd)- This expansion is obviously not useful if one 
has \Kd\ <C Ky. More generally, the large- A?" limit for the C, evolution of / will only be 
useful in kinematical regions where Ky/{NKd) is small enough for A'^ = 3. 
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3.5 CoUinear factorization 

The analysis in the previous sections was concerned with transverse-momentum dependent 
(TMD) factorization, i.e. with cross sections differential in transverse momenta that are 
small compared with the large scale. We now turn to collinear factorization, adequate 
for cross sections with integrated transverse momenta. We point out the main changes 
compared with TMD factorization but do not work out the formalism in detail. 

As in previous sections, we first recapitulate the situation for single Drell-Yan produc- 
tion. The main changes compared with TMD factorization are as follows. 

• We recall from section 3.2.2 that, after a complex contour deformation that avoids 
the Glauber region, the effect of soft gluon exchange is described by a soft factor 
Sq{z). Integration of the cross section (3.35) over q sets z equal to zero in this 
factor. Because the cancellation between real and virtual graphs gives ^^(O) = 1 as 
discussed in section 3.3.1, there is no net effect of soft gluons in the qt integrated 
cross section. 

With the elementary soft factor Sq reduced to unity, subtractions of soft-gluon con- 
tributions as discussed in section 3.2.3 are not required, neither for the parton dis- 
tributions nor for the hard-scattering subprocess. 

• Setting 2 = in the quark and antiquark distributions f{x,z), which is equivalent 
to integrating f{x,k) over k, gives rise to short distance singularities in addition to 
those that are removed by defining the distributions with renormalized quarks fields 
and Wilson lines. The dependence on the ultraviolet subtraction scale /x is described 
by the well-known DGLAP evolution equations. 

The rapidity divergences (from gluons with small i'^ and large i~) that prevent us 
from taking lightlike Wilson lines in the definition of f(x, z) cancel between real and 
virtual corrections when z = [116]. Indeed, the relevant one-loop graphs are those 
in figure 27b and c (without the derivative —d/dyy), and the approximation discussed 
in section 3.4 , which connects these graphs to the soft factor Sq{z), is valid for any 
i~ as long as is small. 

Collinear parton distributions can hence be defined with lightlike Wilson lines and 
do not depend on a parameter Correspondingly, the qt integrated cross section 
is free of Sudakov logarithms. In the operator definition of f{x,fj,), the Wilson lines 



in q{—^z)W^{—^z;v) and W{^z;v)q{^z) merge to a single Wilson line W[—^z, ^ 
given by 



Jo 



(3.113) 



The sections of the paths that go to infinity in {—^z;v) and W{^z;v) have can- 
celled, and a path of finite length between —^z and ^z remains. 

The hard-scattering subprocess now receives radiative corrections not only from vir- 
tual graphs but also from real ones, since emission of partons with large transverse 
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momenta in the final state is permitted once we do not fix the transverse momentum 
qt of the Drell-Yan photon. As already mentioned, a subtraction for the soft-gluon 
region is not required in this case, in contrast to the situation for TMD factorization 
discussed in section 3.2.3. Subtractions are however needed for the regions where 
momenta are collinear to one of the partons entering the hard subprocess. These 
subtractions must be performed in a way that matches the ultraviolet subtractions 
in the parton densities. In particular, the /j, dependence due to collinear subtrac- 
tions in the hard subprocess has to cancel against the fj, dependence of the parton 
distributions in the cross section. 

Let us now investigate the situation for double Drell-Yan production, limiting our- 
selves to a one-loop analysis as we have done throughout the preceding sections. A key 
to understanding the role of soft gluons is to set = in Sqg{zi,y), which results from 
integrating the cross section (3.33) over q^. Our discussion in section 3.3.1 implies that 
5fe(0,/x) = [Ub{0,^l)y'"' = = -Saifi) for the graphs in figure 23. Together 

with the relation Sc{z) = Sb{z), this turns our general one-loop result (3.52) into 

1% - 1 = 1% = s% = , - 1 = 2(1 + c2) S{y, fi) (3.114) 

at zi = Z2 = 0. The one-loop contributions to ^^Sgg cancel between the vertex corrections 
25a and the real graphs 25bi, in full analogy with the case of Sg discussed above. In ^^Sgg 
and ^^Sgg we have a cancellation between the real graphs 25b2 and the virtual graphs 25c. 

We see from (3.114) that in the qt integrated cross section for double Drell-Yan produc- 
tion the contributions from color singlet and color octet distributions decouple from each 
other, and that they have a different behavior concerning soft gluon exchange. In the term 
with color singlet distributions we have a cancellation of soft gluon effects, in full anal- 
ogy to single Drell-Yan production. Also, the graphs for the hard-scattering subprocess are 
exactly as in the single Drell-Yan process, with a cancellation of the soft-gluon region but 
with necessary subtractions for the regions of collinear parton momenta. From our above 
discussion it follows that collinear two-parton distributions can be defined with the same 
operators as their single-parton analogs, with lightlike Wilson lines W[y — ^zi,y + ^zi] and 
W[—^Z2, \z2\ between quark and antiquark fields, and that their contribution to the in- 
tegrated cross section is free of Sudakov logarithms. Like their single-parton counterparts, 
the distributions have ultraviolet divergences; the scale dependence that follows from their 
subtraction will be discussed in section 5.3.2. 

The contribution of collinear color-octet distributions ^F is quite different. Because 
real and virtual graphs have different color factors, soft gluon effects do not cancel between 
them, and their net effect is described by ^^Sgg{y). As a consequence, the different factors in 
the cross section formula require soft subtractions, as they do in the case of measured trans- 
verse momenta. Since the color indices of [q{—^Z2)W'^{—^Z2]v)]y and [VF(|z2; f )(?(|^2)]fc 
are not contracted, the two Wilson lines do not merge into a single one of finite length, and 
the same holds for their analogs with arguments y — ^zi and y + ^zi. The vector v in the 
Wilson lines cannot be taken lightlike, so that collinear octet distributions will depend on a 
parameter The resulting Collins-Soper equation gives rise to Sudakov logarithms, which 
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suppress the color octet contribution to the qt integrated cross section. This important 
resuh was aheady obtained in [121], based on the observation that in the hard-scattering 
subprocesses there is no cancehation of the soft-gluon region. An adequate scale /io for 
the initial condition of the CS equation will in this case be a hadronic scale, inverse to the 
typical distance \y\ between the two scattering partons. 

Let us finally take a look at the interference distributions Iai,a2- From (3.56) we find 
a soft factor 

S,= ('„ °)+2Sto,,)(; (3.115) 

at zi = Z2 = 0. There is hence no cancellation of soft-gluon effects, so that a formulation of 
collinear factorization will in this case be similar to the one for the color octet distributions 
just discussed, with the additional complication of mixing between the color singlet and 
octet channels. 



4 Some properties of two-quark and quark-antiquark distributions 
4.1 Spin structure 

Multiparton distributions have a nontrivial spin structure because the polarizations of 
different partons can be correlated among themselves, even in an unpolarized proton. In 
the following two sections we first investigate some general properties of spin correlations 
between two quarks and then show that they have observable consequences in multiple 
scattering cross sections. We will encounter several examples for parton spin correlations 
in section 5.2.2. 

4.1.1 Spin decomposition 

Let us first take a closer look at the spin dependence of the two-quark distributions i^ai,a2 
introduced in (2.86), making use of rotation and parity invariance. We always assume that 
the hadron is unpolarized, i.e. that the matrix element (2.86) is averaged over the hadron 
spin. The simplest cases are the distributions 

Fq,q{Xi, ki, y) = fq^q{xi, X2, ki, k^, ^1^2) ^lUi ^2 2/i V^) 1 

FAqAq{xuki,y) = fAq,Ag(x'i-^X2,kj,kl,k^k2,k^y,k2y,y'^) , (4.1) 

which are parity even, i.e. scalar functions. By contrast, the distributions Fq^j\q and F^q^q 
are parity odd, i.e. pseudoscalar functions. Their general form is 

Fq,Aq i^iiV) — 2/"' fq,Ai ~^ ^2 V'^ fq,AQ ~'~ ^2 fq,Ag ' 

F^^^^{x,, k,,y) = e^^'ki y^' /l,,, + e^^'k^y^' /I,,, + e^^'k\ k( /I,,, , (4.2) 

where and are scalar functions with the same arguments as in (4.1). The scalar 

functions are in general neither even nor odd in ki or y since their dependence on kiy, k2y 
and k^k2 is not constrained by symmetry. Note that the three two-dimensional vectors 
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y, ki and k2 are linearly dependent, so that the three cross products in (4.2) are linearly 
dependent as well. Expressing e.g. y as a linear combination of ki and ^2 one obtains 

jfU f _ fci {k2y)-{k^k2){k^y) ., ■ . 
1^ ~ (k k ^2 ^ 1 2> 



Jf u ,J' _ fc2 {k,y)-ik,k2){k2y) • . 



and can thus write Fq^/^q as e^^'k-^k-i^ times a single scalar function. However, that scalar 
function is singular when ki and ^2 become collinear, as is evident from the denominators 
in (4.3). To avoid such artificial singularities, one can use (4.2) if it is necessary to make 
the appearance of e^^' explicit. 

Using that ^(1 ±75) projects on quarks with definite helicity, one can readily identify 
the combinations of quark polarizations that are described by the above functions. In a 
schematic notation one has 

Fq,q ^qfqt + (1102 +Q1Q2 + Q1Q2 ' 
FaqAq ^ Qi^lt + Q1Q2 -Q1I2 - lilt ' 

FqAq ^ it 0.2 - 1112 -Q1I2 + € it ' 

FAq,q ^ qt(l2 - (I1Q2 + Qi(l2 - (I1Q2 ^ (4-4) 

where the superscript in denotes the quark helicity. The distribution FAq,Aq thus 
describes the degree to which the two quark helicities are aligned rather than antialigned, 
whereas Fq^^q and FAq,q describe the correlation between the helicity of one of the quarks 
and one of the cross products in (4.2). 

To illustrate that spin correlations between two partons need not be small, let us 
consider the simple case of a SU (6) symmetric three-quark wave function of the proton. 
Its spin-flavor part reads 



1 



V6 



-u-d+) + |ti-u+d+) - 2|u+M+(i-)) , (4.5) 



where + and — respectively indicate that the quark spin is aligned and antialigned with 
the proton spin. As is well known, this wave function gives Au/u = 2/3 and Ad/d = — 1/3 
for the longitudinal polarization of u and d quarks, which reproduces at least the trend of 
what is empirically found for the lowest x moments of the polarized quark densities. For 
two-quark distributions one finds 

Fau,Au/Fu,u = 1/3 , F^^^^d/Fu4 = -2/3 (4.6) 

and thus an appreciable correlation between the longitudinal polarizations of the quarks. 

Of course, the study of a three-quark wave function tells us little about partons with 
X ~ 10^^ or smaller, which are of particular relevance for LHC phenomenology. To the 
extent that they are known, polarized single-parton densities in this x range are small 
compared with their unpolarized counterparts, which means that there is only a weak 
spin correlation between a small-x quark and the proton as a whole. This is not too 
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surprising, given that a small-x quark and the proton are far apart in phase space. It 
does however not imply smah spin correlations between two quarks that have small but 
comparable momentum fractions xi ~ X2 and are thus closer in phase space. How large 
such correlations are is an important open question. 

The distributions defined with one or two tensor operators Og^ are associated with 
transverse quark polarization, since ^(7''' + sHa^~^j^) projects on quarks with a transverse 
spin vector . We now discuss their parametrization in terms of scalar or pseudoscalar 
functions. Let us begin with F^^g^ and which transform like two-dimensional 

vectors. They can hence be written as a sum of three scalar functions that are respectively 
multiplied by , k-[ and Only two of these functions are independent because of the 
linear dependence of the three vectors. A minimal parametrization is obtained by taking 
and 

as basis vectors. This gives 

PsqAqi^i^ ^^^y) = y'^fSqAq + V^^QSqAq > (4-8) 

where we have inserted the proton mass M on the r.h.s. so that / and g have the same 
mass dimension as F. Here and in the following we denote scalar functions by / and 
pseudoscalar ones by g. The latter can be represented in the same way as Fq^^q. If one 
wants to avoid pseudoscalar functions, one can replace the basis vector y^ by 

kl = ik, + k,y-^-^l±p^y^. (4.9) 

Since both y and ky are orthogonal to y, they must be proportional to each other, and 
explicitly one finds 

If one inserts this into (4.8) then ky is multiplied by scalar functions, which are however 
singular when ki + k2 and y become collinear. One could replace fei + ^2 in (4-9) by 
another linear combination of ki and but this would only move the singularities to a 
different part of phase space. 

The distributions F^^^ and F^^^ transform like axial vectors, so that one has 

^i&qi^i^ ^i^y)= y^Mfq^Sq + y^Mgq^Sq , 

Fiqjx,, ki,y) = y^Mfsq,q + yndgsq,q . (4.11) 

A decomposition in terms of scalar functions can be obtained by replacing y^ with e^^'ky. 
The tensor distribution Fg^^^ can finally be written as 

Fi[sqi^^, K y) = ^"'fsqM + {^y' y'' - ^"'y') M'f^Sq 

+ {y'y'' + y'y'') M^glq^gq + {y^y^' - y^y^') M^glq^gq . (4.12) 
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Notice that the four basis tensors t-p in this decomposition are orthogonal to each other, 
i.e. t-p t-q oc 6pq with p,q = 1,2,3,4. Contracting (4.12) with the transverse polarization 
vectors sjsg of the quarks, we see in particular that fsq^sq goes with S1S2 and thus describes 
the correlation between the two transverse quark spins. 

In summary, we can represent the spin structure of Fa^^a2{xi^^i^y) by eight scalar 
and eight pseudoscalar functions for each combination of quark flavors. The pseudoscalar 
functions can be traded for scalar ones, which have however artificial singularities for 
particular values of the vectors y, ki and fc2- 

If one integrates over transverse momenta to obtain collinear distributions, one finds 

F/^q,q{Xi,y) = Fg^Aq{Xi,y) = (4.13) 

because one cannot construct a pseudoscalar function with only one vector y. Likewise, 
the functions g in (4.8), (4.11) and (4.12) integrate to zero, so that we are left with 

Fq,q{Xi,y) = fq,q , FAq,Aq{Xi,y) = fAqAq ' 

Ksqi^i^y) = y'Mfq^Sq , Fi^,Sq^X,,y) = y^Mf^q^ , 

^Sq^M^y^ = y^MfSq,g ' HqAq^^'^y^ = y^^^fSlAq ' 

y) = S"'fsq,Sq + ( Vy^' - S"'y') M'fsq,Sq ■ (4-14) 

The eight functions / on the right-hand side now depend on xi, X2 and y^ and are obtained 
from their counterparts in (4.1), (4.8), (4.11) and (4.12) by integration over ki and fe2- 

The above decompositions are given for distributions in a right-moving proton. For a 
left-moving proton one has to change the sign of e^^ and hence of y and of all pseudoscalar 
functions, see our remark below (2.100). Analogous decompositions can be written down 
for distributions Fai,a2 5 -^ai,a2) -^ai,a2 that involve antiquarks, as well as for interference 
distributions /ai.oi and Iai,a2- 

Symmetry properties The terms appearing in the decompositions (4.1) to (4.14) are 
consistent with rotation and parity invariance. Let us now discuss their symmetry proper- 
ties. Using that the operators in (2.80) satisfy 0*{yi,Zi) = 0^{yi, —Zi) one finds that the 
distributions Fa^^a2{xi-,ki,y) are real valued, 

Ki,a2i^i,f^i,y) = Fai,a2iXi,ki,y)- (4.15) 

For distributions that are purely defined in momentum or position space, see (2.9) and 
(2.12), this implies 

Kua2(^i^^i^y) = Fai,a2iXi,-Zi,y) . (4.16) 

These functions are in general not real-valued, nor are the scalar or pseudoscalar functions 
one can introduce to parameterize them in analogy to (4.1) to (4.12). 

For the symmetry properties of parton distributions with respect to time reversal, the 
Wilson lines appearing in their definition are essential. As we argued in section 3.2.1, 
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multiparton distributions involve the past pointing Wilson lines W given in (3.13). Upon 
time reversal these turn into the future pointing Wilson lines W in (3.17). A distribution is 
called T even (odd) if it is even (odd) under time reversal without taking into account this 
change of Wilson lines. The time reversal invariance of strong interactions thus implies that 
T odd distributions are only nonzero thanks to Wilson line effects; a prominent example 
from spin physics is the Sivers distribution function [122]. However, time reversal does 
force distributions to vanish if they are T odd and have Wilson lines that are invariant 
under time reversal. This is the case for the Wilson lines along a finite lightlike path that 
appear in collinear single-parton densities and in the collinear two-parton distributions 
in the color singlet sector, as we discussed in section 3.5. By contrast, collinear color octet 
distributions ^F, as well as interference distributions and have Wilson lines that do 
change under time reversal. 

After these preliminaries we can now investigate the time reversal properties of two- 
quark distributions. We find 

PZa2(^i^^i^y) = "naiVai F^'aii^r, h, -y) (4.17) 

with sign factors r/g = +1 and r/Ag = VSq = where the superscripts indicate the type of 
Wilson line in the matrix element defining the distributions. The relations (4.15) to (4.17) 
also hold for the distributions -Fai,a2i -^ai,a2) -^ai,a2 with antiquarks and for the interference 
distributions /ai.aa, hi,a2- 

Since the scalar functions parameterizing F^ ^^{xi, ki,y) are in general neither even 
nor odd in y, they are not T even or odd either. The scalar functions that parameterize 
the collinear distributions in (4.14) are however even in y. As a consequence, F^^ ^^{xi,y) 
and ^^{xi,y) are T odd and all other distributions in (4.14) are T even. For the color 
singlet sector this implies 

'FA,M^^^>y) = H,A,(^^>1/) = 0, (4.18) 

whereas the corresponding color-octet distributions ^F^^ g^{xi,y) and ^F^^ ^^{xi,y) can 
be nonzero due to the Wilson lines appearing in their definitions. Analogous statements 
hold for the corresponding distributions with one or two antiquarks. Collinear interference 
distributions and are not restricted by time reversal invariance. 

4.1.2 Spin effects in gauge boson pair production 

In this section we show that the quark spin correlations discussed in the previous section 
have observable consequences in multiparton interactions. As we did earlier in this paper, 
we consider the production of a pair of gauge bosons 7, Z or W. We include the decay 
of each boson into a lepton pair, which carries information on the spin state of the gauge 
boson. While these processes have a rather small cross section, they may be suited for 
experimental studies due to their clean final-state signature. We do not present a full 
analysis here, but highlight the effects of selected parton spin correlations. 

For simplicity we limit our attention to those distributions that do not involve explicit 
vectors y or j/ on the r.h.s. of the decompositions in the previous section, i.e. to Fq^q, F^q^/\q 
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and the term S^^'f^^ in F^^ . For definiteness we analyze the graph of figure 6a, with 
two quarks emitted from the right-moving proton and two antiquarks from the left-moving 
one. We approximate the transverse momenta of the bosons by zero when calculating 
their production and decay, as deviations from this limit are suppressed by powers of qt/Q- 
The partonic cross section for the production of a lepton pair can be written as the product 
of a production tensor of the boson and a tensor for its decay, 

^a,a = «-ar''' ^mm' ' (4-19) 

where fi is associated with the boson in the amplitude and //' with the one in the complex 
conjugate amplitude. For unpolarized or longitudinally polarized quarks one easily finds 

{^g,gr' = -i^AgAgr' = -AqT' - BieT' , 
(^A,,,-)'^'^' = -i^gA^r^' = -BgT' - Ai^T' , (4.20) 

with coefficients A and B depending on kinematic variables and the vector and axial- 
vector couplings of the gauge boson to the quark q. In the case of a photon one has 
5 = 0. The transverse tensors and have as nonzero components g^ = g"^^ = — 1 
and e^^ = — e^^ = 1. From (4.20) it follows that the overall cross section depends on the 
combinations 

Fq,q Fg^g + FAg,Aq FAgAg ^ud Fq^qFAqAq + FAqAqFq,q (4.21) 

of multiparton distributions. The contraction of g^ and e^'^ with the vector boson decay 
matrices results in different angular distributions of the leptons. If one integrates over their 
angles then only the contribution from g^ remains. We thus find that nonzero values of 
FAqAg ^^"^ FAqAq modify both the total rate and the lepton angular distribution compared 
with the contribution from the unpolarized term Fq^q Fq^q. 

We now turn to transverse quark polarization. In this case the production tensor 
from each hard scattering depends also on the transverse indices associated with the quark 
polarization. We recall that the quark field bilinears qia'^^'j^q are chiral odd, so that in 
the helicity basis they correspond to quarks or antiquarks with opposite helicities in the 
scattering amplitude and its conjugate. As a result transverse quark polarization does 
not contribute to the production of a W boson. Keeping only the term d^^'fsq,Sq in the 
decomposition (4.12) of F^j^^^ and the corresponding term in the decomposition of Fg^^-, 
one finds for the neutral bosons 7 or Z 

^'^'S''' i^i'sq^qr' i<tsqr' ^ ^^^^ + 5^5^^ " ^f'^f > (4-22) 

where the polarization indices fi, fi' belong to the boson with momentum qi and the indices 
z/, u' to the boson with momentum q2. We observe that the polarization indices of the two 
bosons are entangled in (4.22). Contracting with the well-known boson decay matrices, 
one obtains an azimuthal dependence like cos(299), where (p is the relative azimuthal angle 
between the two leptons (as opposed to the antileptons).^^ We thus obtain the important 



Since we are working in the approximation Qi = q2 = 0, the azimuthal angles for the leptons of both 
boson decays are naturally defined w.r.t. the z axis in the pp center-of-mass. 
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a b c 



Figure 28. Parton-level graphs for qq annihilation into four leptons. Further graphs are obtained 
by interchanging the leptons £i and £2 together with the corresponding antileptons in a and c, and 
by reversing the charge flow of the central fermion line in c. Graph b requires at least one of the 
produced bosons to be a W . 

result that a correlation between transverse quark and antiquark spins, as expressed by 
the distribution fsq,5q in (4.12), leads to a correlation between the decay planes of the two 
produced bosons. 

It is instructive to rewrite the production tensors in terms of boson polarization vectors 
e+ = —(0, 1, i, 0) / \/2 and e_ = (0, 1, —i, 0) / \/2, which respectively correspond to angular 
momentum +1 or —1 along the z axis. One finds 

-qT' = 4 + e'^ e-''' , -i^T' = 4 - e- e*- ' (4-23) 

and 

\ 5^ ^' + g'/gi" - gT'aT') = < e*/ + e*/ e\ 6*/ . (4.24) 

We can easily understand why each boson is transversely polarized. Recall that a massless 
quark and antiquark can only annihilate into a vector boson if their helicities are coupled 
to ±1. Since we neglect the transverse momentum of the bosons, their angular momentum 
along z must also be ±1. The tensors and ie^ in (4.23) correspond to the same 
boson polarization in amplitude and conjugate amplitude and thus do not give rise to an 
azimuthal dependence in the leptonic decays, but they do give different distributions in the 
polar angles of the leptons, or equivalently in their rapidities. By contrast, (4.24) involves 
the interference between = 1 and —1 for each of the bosons, which readily translates 
into the cos(29j) dependence already mentioned. 

It is natural to expect that spin correlations between partons also lead to angular 
correlations in the final state for other double-scattering processes, such as the production 
of two dijets. In this case two-parton distributions involving linear gluon polarization can 
contribute as well. We note that in the analysis of [9] uncorrelated dijet planes were taken 
as a characteristic feature of the double-scattering mechanism. This is only adequate if 
parton spin correlations in the proton are negligible. 

Returning to four-lepton production, let us compare our results for double hard scat- 
tering with the contribution from a single qq annihilation, remaining in kinematics where 
Qi and ^2 can be neglected compared with Q. The corresponding Feynman graphs involve 
either quark exchange in the t ox u channel, or an intermediate boson in the s channel 
in case one or both final-state bosons are charged, see figure 28a and b. In addition, the 
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four-lepton final state can be produced by graphs as in figure 28c, wliere only one lepton 
pair comes from the decay of a vector boson. Such graphs were recently discussed and 
termed "single-resonance graphs" in [13]. They should be taken into account unless each 
lepton pair has an invariant mass inside the Z or W mass peaks. 

The dependence of the single-scattering cross section on the azimuthal angle (p between 
the leptons can be deduced from symmetry arguments since we set Qi and ^2 to zero in the 
parton-level process. Let us assume that the initial qq pair has total angular momentum 
= 1 along the z axis. This must then hold for both the scattering amplitude and its 
conjugate, since a single quark or antiquark in an unpolarized proton is unpolarized. One 
of the lepton pairs originates from the decay of a gauge boson and is thus in a partial 
wave with J = 1. The possible combinations of {Jf, Jl ) for the two lepton pairs are thus 
(2,-1), (1,0), (0,1) and (—1,2), where the first and last possibilities are only possible 
for single-resonance graphs. In the cross section this gives a ip independent term and a 
modulation with cos(p from all graphs, as well as modulations with cos{2ip) and cos(3(^) 
from single-resonance graphs. The same angular terms are obtained when the initial qq 
pair has = —1, whereas the configuration where the pair has = decouples in 
the hard-scattering graphs of figure 28. In summary, the cos(2(/3) modulation we found in 
the double-scattering mechanism competes with single-scattering contributions involving 
single-resonance graphs. One may hope that the two sources of cos(2(/3) dependence 
can be separated by a more detailed analysis — making for instance use of the fact that 
the single-resonance graphs also give a cos{3ip) term — but this issue must be left to a 
dedicated study. 

4.2 Mellin moments and lattice calculations 

If one takes Mellin moments in xi and X2 of the color singlet distributions ^F, then the 
light-cone operators Oa in their definition turn into local operators. The corresponding 
moments of single-parton densities can be calculated in lattice QCD, and we will now 
investigate to which extent the same can be done for two-parton distributions. From our 
discussion in section 3.5 it follows that the Mellin moments of color octet distributions 
do not involve local operators because of their Wilson line structure. We therefore limit 
ourselves to the color singlet sector. 

Using the relation (2.88) and its analogs for Fa^^a^ and Fa^^a2 5 one obtains double 
Mellin moments 

M2,\Z'iy') = Td^i^r"' [' dx,x^'-'\'F,,,,,{x,,X2,y) + {-iraa,'F^,,a,{xi,X2,y) 
Jo Jo L 

+ i-lT^ aa,^Fa,,a,{xuX2,y) + i-lT'-^''^ aa,aa,^Fa,-a,{xi,X2,y) 

= 2(p+)^-"-"^ I dy^{p\ O+-+(0) <■■+(?/) |p),+=o (4-25) 

with (Tq = asq = +1 and cTAg = — 1. For each label dq the moments Mal'a2 and the corre- 
sponding operator on the r.h.s. carry an additional transverse Lorentz index not displayed 



*This was not realized in [47], where single-resonance graphs were not taken into account. 
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in (4.25). On the r.h.s. we have the twist-two quark operators famihar from the operator 
product expansion, 



0m- 




= T 






= T 






T 



S q{y)^^'HD^'^{y)■■■iD^'-{y)q{y), 
S ^ qiy) 7^*^75 ^B^- (y) • • • iB'^" {y) q{y) , 
A S «-(y)i(T^^i75iB'^^(y)---iB'^"(y)(?(y), (4.26) 

where S'^(y) = ^{d'^ — d^) + igA^{y) is the antisymmetrized covariant derivative. S, A 
and T respectively denotes symmetrization, antisymmetrization and subtraction of traces 
in the indicated indices. It is understood that the composite operators Oai (y) and (0) 
in (4.25) are each renormahzed in a standard manner, e.g. in the MS scheme. As long as 
they are taken at a finite spacelike distance y, no further ultraviolet divergences appear. 
The operators Oq and Oaq mix of course with their gluon counterparts Og and O^g 
under renormalization. We shall return to the renormalization of collinear multiparton 
distributions in section 5.3.2. 

Let us rewrite the r.h.s. of (4.25) in a manifestly covariant form. We first introduce 
the covariant decomposition 

{p\ O,""-""^(0) O^q'-^"^ {y)\p) = 2p'^^ • ■■p''-2 pMi . ..pf^^, ipn, o^^-){py, y2) + . . . , (4.27) 

where the ellipsis represents terms with uncontracted vectors y^ and terms involving the 
metric tensor g^" . The reduced matrix element (Og^ Og^) can only depend on the invariants 
py and y^ . We then choose a frame where p = and y"*" = 0, so that py = p~^y~ and 
= —y^. This allows us to write (4.25) in the desired form 



M™(-y^) = J d{py) {O^q^ 0^^)ipy,y') . (4.28) 

A corresponding representation is readily obtained for M^^^^. For one or two polarization 
labels 5q the analogs of (4.27) involve the tensor structures in the decomposition (4.14) of 
the two-quark distributions. 

The matrix element in (4.27) and its counterparts with polarized quarks can be eval- 
uated on a lattice in Euclidean spacetime if one chooses y^ = 0. This is rather similar 
to lattice studies of transverse-momentum dependent single-quark distributions [123, 124], 
with the main difference that the operators taken at different spacetime points are single 
quark fields in that case, whereas they are gauge invariant bilinear operators here. The 
restriction y^ = entails 

^ = ^<P-, (4.29) 

_y2 y2 

where p and y denote the spacelike three- vectors. Thus, the integral over all py at fixed 
y^ on the r.h.s. of (4.28) can unfortunately not be evaluated from results on a discrete 
Euclidean lattice, where the maximal momentum is fixed by the lattice spacing. This is 
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completely analogous to the single-quark discussed in [124]. Despite this limita- 

tion of principle, we hope that lattice data in a certain range of py and will in the 
future provide genuinely nonperturbative information about the behavior of multi-parton 
distributions. 

We note that a lattice calculation has been reported in [125] for the correlation function 
of two vector currents at equal time in a proton at rest. This corresponds to setting 
ni = n2 = 1 and = in (4.27). The reduced matrix element (O"^ ) at py = is 
directly related to an integral of the two-quark correlation function $(fci, k2,r) defined in 
(2.75). The relative plus-momentum r"*" is integrated over in that case, rather than being set 
to zero as required for the distributions F{xi,ki,r) that appear in double hard-scattering 
cross sections. 

4.3 Relation with generalized parton distributions 

In section 2.1.5 we derived an approximate relation between multi- and single-parton distri- 
butions in a model theory with scalar partons. We now extend this relation to distributions 
for two quarks or antiquarks, taking into account the different combinations of fermion 
number and color. For the time being we neglect aspects related to the proton spin, which 
will be discussed in section 4.3.1. The distributions we will deal with are 

^Fai,a2 = {{{q3'^a2q2) {Qi'^^ai qi))) , ^Fai,a2 = {{ {Qi'^' ^2) (gsTaiQl))) , 

^Fa-, ,a2 = (( (92 Taa ^s) (^4 ^ai 9l ) )) , ^^ai ,a2 = {{ (^4 Taj Qs) (^2 ^ai Ql))) , 

%i,a2 = {{{q2^a2qi) (^sTai^l))) , ^/ai,a2 = {{{q'i^a2qA) {q2^aiqi))) ■ (4-30) 

The general result (2.69) for n scalar partons readily carries over to the two-quark 
distributions ^F: 

^Fa^^a2{xi,Zi,y) J cP'b fa2{x2, Z2;b + \xiZi) fa^{xi, Zi^b + y - \x2Z2) , (4.31) 

where the impact parameter dependent single-quark distributions /a(x, z; b) are defined in 
analogy to the scalar case in (2.64) and (2.65). Setting zi = Z2 = in (4.31), one obtains 
collinear distributions on both sides and has the probability interpretation represented in 
figure 4. As a counterpart to (2.72) one can transform the relation (4.31) into transverse- 
momentum space, where it reads 

^Fai,a2{xi,ki,r) ^ fa2{x2,k2 - \x2r; -r) fa^{xi,ki - ixir;r) (4.32) 

with distributions fa{x,k;^) defined in analogy to (2.66). Integrating over fcj we obtain 
the relation recently given in [59]. Relations analogous to (4.31) and (4.32) are obtained 
for ^Fa^^a2 by replacing the label 02 with 02 on both sides. 

To reduce ^F to single-particle distributions we could repeat our earlier derivation that 
started with (2.59). We find it more convenient to work in the transverse-momentum rather 
than impact parameter representation. We insert a complete set \X) of intermediate states 



-93- 



between the color singlet operators (94 92) and {qsTaj qi) and assume that single-proton 
intermediate states dominate. We then have 



dz- d „•„ y ^ -\- „• 



(2vr) 



■ z- p'^—iz k- 



X " '-4=1 

X {p+ ,p\Oa,{y,z2)\x) {x\OaAy.zi)\p'^ ,p) 



p=0 



dp'+ d'^p' 
2p'+ (27r)2 



2p+/ dy^d^y e'^ 



j-j- I dz- d^Z^ p+^^z^k^ 
i=l 



(2^)2 



x{p+,p\Oa,{y,Z2)\p'+,p'){p'^,p'\Oa,{y,Zl)\P~',P) ■ (4.33) 

p=0 

After an appropriate shift of the position arguments in the bilinear field operators and a 
change of integration variables from y, zi, Z2 to uq = ^{zi — Z2), ui = ^{zi + Z2) + y and 
U2 = ^{zi + Z2) - y this gives 



^Fa„aM,h,r) « ^ / dp'+dy [n / 



du- d^ 



(27r)3 



Jy-o (p'- [1-^1+2:2] p)^-«'"o(p'-*'i+*'2) 



X e 



{p+,p\Oa2{0,U2)\p'^,p') {p'+,p' I Oai(0,ni) |p+,p> 



p=0 



(4.34) 



The integrations over Uq and Uq fix the momentum of the intermediate state. In partic- 
ular, its plus-momentum is p'~^ = {1 — xi +X2)p'^, which reflects that the operator qsTa^ qi 
describes the emission of a quark with plus-momentum xip'^ and the reabsorption of a 
quark with plus-momentum X2P'^. The matrix elements in the approximation (4.34) are 
thus given by generalized parton distributions (GPDs), which play a prominent role in the 
description of hard exclusive processes, see [126-128] and the reviews [129-131]. To evalu- 
ate the unapproximated form in (4.33) one would need the corresponding matrix elements 
for all transitions p ^ X. This is obviously impractical, although for selected transitions to 
single baryons, e.g. for p — t- A(1232), some information can be obtained [129, 131]. GPDs 
are defined by 



fa{x,i,k-p,p') 



dz 



Axz 



dh 



(2vr) 



e-'^'' {p'+,p'\Oa{0,z)\p-^,p), (4.35) 



where P = ^{p + p') and Tq is one of the matrices in (2.81). The parameter 



p-i" — p 
p+ _|_ p' 



/+ 



/+ 



(4.36) 



is often called skewness. One finds that k is the average transverse momentum of the two 
quark legs and x their average plus-momentum divided by the average plus-momentum 
of the proton states. For ease of notation we do not indicate the polarization states of the 
protons, which are in general different. A parameterization of the matrix elements (4.35) 
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for spin 1/2 hadrons in terms of scalar functions can be found in [132]. Invariance under 
the transverse boost specified hy v ^ v — {v~^/P~^) P gives the relation fa{x,i, k;p,p') = 
fa{x, i,k- xP- A, i A) with A = (1 + i)p' - (1 - i)p. Abbreviating 

fa{x, fc; A) = Ux, e, k- -i A, i A) (4.37) 

we can thus rewrite the relation (4.34) as 

1 + 

'^Fa^^a2{xi,ki,y) ^ ^— ^ /^^(a;,-^, |(fc+ + r + xfc_);(l+^)fc_) 

X fa, (x, e, \{k+ -r + xk_)- -(1 + e)fc-) , (4.38) 

where k± = fci it ^2 and 

. = ^i±^. S = ^f^i^. (4.39) 

2- XI+X2 2 - Xi + X2 

In complete analogy one derives 

1 + C 

\,,a2{3^i,ki,r) s» Y—^fa2{x,-S,,^{k+ + r + xk_)-{l + ^)k^) 

X /,,(x,e,i(fc+-r + xfc„);-(l + e)fc-), (4.40) 

where the generalized parton distributions fa for antiquarks are defined by replacing 
Oa{0,z) with OaiO,z) in (4.35). One finds 

fa{x,^,k;A) = aafa{-x,^,-k;A) (4.41) 

with the same sign factors aq = asq = +1 and fiAg = —1 that appeared in (2.88). We also 
note that generalized parton distributions with positive and negative skewness parameter 
are easily related to each other by taking the complex conjugate of (4.35). 
For the distributions ^Fai,a2 and ^Iai,a2 we obtain 

1 + f 

^Fa,,a2{xi,ki,r) « ^— ^ /a2(-2;,-C,-i(fc- +r + xfc+);(l + ^)/c+) 



X 



fa,{x,^,^{k_ - r + xk+y,-{l + ^)k+) 



1 + c 

%,,a2{xi,ki,r) ^— ^ /^^(x,-^, i(fc_+r + xfe+);(l+OA;+) 

X fa, {x, ^, i(fc_ - r + xfc+); -(1 + k+) , (4.42) 
where we have again k± = ki ziz k2 but now 



2 — Xi — X2 2 — Xi — X2 

In this case we have < ^, which describes the emission of a quark-antiquark pair. Again, 
this could be anticipated from figure 5 since now the parton lines combined to color singlets 
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are {12} and {34}, with each pair being on the same side of the final-state cut in the double 
parton distribution. 

An important difference between the approximations for ^F, ^I, I and the one for 
given in (4.31) is that the generalized distributions on the r.h.s. of (4.38), (4.40) and (4.42) 
do not reduce to collinear functions if we integrate over ki and ^2, because these momenta 
appear in their fourth arguments. 

The preceding derivations can easily be extended to distributions describing inter- 
ference between different quark flavors. The distributions corresponding to figure 7 are 
defined with bilinear operators dVu or uTd. For a proton target, the ground state in 
the sum over intermediate states inserted between the two operators is then a neutron. 
Isospin symmetry relates the resulting matrix elements to matrix elements in the proton: 
{n\dTu\p) = {p\uTd\n) = {p\uTu\p) — {p\dTd\p). Under the assumption of 5?7(3) fla- 
vor symmetry one can derive similar relations for distributions involving strange quarks 
[129, 131]. 

Although the relation between multiparton distributions and GPDs is an approxima- 
tion whose accuracy is not easy to estimate (and although our current knowledge of GPDs 
is far less advanced than that of ordinary parton densities) this relation provides opportu- 
nities to obtain information about multiple interactions that is hard to get by other means. 
One example are the different interference distributions discussed above, which are so far 
entirely unknown. Perhaps even more important is that GPDs give rather direct informa- 
tion about the distribution of single partons in the impact parameter b, which is Fourier 
conjugate to a transverse momentum transfer A that can be measured in physical pro- 
cesses. This is in stark contrast to the interparton distance y in two-parton distributions, 
which appears as an integration variable in cross section formulae like (2.91) and is not 
directly related to observable kinematic quantities. We already mentioned in section 2.6 
that studies of GPDs give evidence for a correlation between the longitudinal momentum 
and the impact parameter of partons in the proton. For values of the momentum fraction 
where such a correlation is strong, it is hardly plausible that there should be no correlation 
between xi,X2 and y in two-parton distributions, even if there were important corrections 
to approximations like (4.31). 

4.3.1 Spin correlations 

We now take a closer look at the role of the proton spin in the approximate relation 
between double and single parton distributions, which we have glossed over up to now. For 
our purpose, a suitable choice to describe the spin state of a proton is the light-cone helicity 
A = zb^, which is equal to the usual helicity in a frame where the proton plus-momentum 
p"*" tends to infinity (see e.g. [133] or [130, section 3.5.1]). We denote the corresponding 
momentum eigenstates by A). 

When inserting intermediate proton states between the two color singlet operators in 
a two-parton distribution for an unpolarized proton, we schematically have 

l^{p+,P,X\Oa, Oa, b+,P,A> 
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or an analogous relation with states \p^, b, A) of definite transverse position. The sum over 
states on the r.h.s. thus includes single-parton matrix elements where the proton helicity 
differs in the bra and the ket state. Corresponding sums over polarization states should 
hence be inserted in the relations (4.31), (4.38), (4.40) and (4.42). 

To discuss the implications of this observation, let us focus on the collinear distribution 
^F(xi,r). At the level of matrix elements we have a relation 



'Fa„a,{x^,r) ^ ^J2f^f{x2;-r)f^['\xv,r), (4.45) 



A,A' 



where the superscripts A and A' of fa denote the proton helicities as in (4.44) and an 
average over proton helicities is understood in ^F. A standard decomposition of the spin 
dependence of generalized parton distributions involves two distributions H and E for 
unpolarized quarks and two distributions H and E for longitudinally polarized quarks. 
The distribution E does not contribute in the case of zero skewness ^ = we are dealing 
with here. Using the conventions and the matrix elements for definite proton light-cone 
helicity in eq. (54) of [130], we have 

f++{x,r) = H''{x,0,-r^), /""(x, r) = F«(x, 0, -r^) , 

1 ■ 2 1 ■ 2 

fq^ix,r) = L±^E^{x,0,-r'), f^-{x,r) = E\x,Q,-r') (4.46) 

and 

/++(x, r) = H^{x, 0, -r^) , J^" (x, r) = -H^x, 0, -r^) , 

f-+{x,r) = 0, /+-(x,r)=0, (4.47) 

where M is the proton mass and H'^{x,^,t), E'^{x,^,t) and H'^{x,^,t) are the usual GPDs 
defined in [130]. H'^ and H'^ are the respective generalizations of the unpolarized and 
longitudinally polarized quark densities q and Aq. Changing the basis of the proton spin 
states, one can see that E*^ is related to unpolarized quarks in a transversely polarized 
proton [65]. Inserting (4.46) into (4.45), we get 



2 



%,g{x,,r) « H''{X2,0, -r') H'^{xi,0, -r') + ^ i?'?(x2, 0, -r') E^xi,0, -r') . (4.48) 

The term with H corresponds to the simplest approximation of the two-parton distribution 
as a product of single-parton distributions, whereas the one with E appears in addition. 
E describes a correlation between the position of a single quark and the proton spin, 
and (4.48) shows how such a correlation may lead to a correlation between two quarks 
in an unpolarized proton. It is difficult to say whether this correction term alone already 
provides an improved approximation of ^Fq^q{xi,r), but one may take it as an indicator 
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for the possible departure from a simple factorized ansatz that neglects all correlations. 
In cases where the correction from E is large, it is plausible to expect that the factorized 
ansatz involving only H will fail. 

Applying the same method to the distribution ^-FAg,A<j) we obtain from (4.47) 



Again one should be cautious regarding the validity of this approximation. For the region 
of small but similar xi and X2 we already argued in section 4.1.1 that one may well have 
sizeable correlations between the longitudinal polarization of two quarks, even if there is 
little correlation between the longitudinal polarizations of one quark and the proton as a 
whole. Conversely, in kinematics where the product of quark-proton spin correlations in 
(4.49) is sizeable it seems natural to assume that quark-quark spin correlations are sizeable 
as well. 

5 Perturbatively large transverse momentum 

So far we have treated multiple interactions as a two-scale problem, in which the virtualities 
11112 ~ define a large scale whereas the transverse momenta \qi\, 1^2 1 scale A 

of nonperturbative interactions are treated as small. We now make a distinction between 
the different scales previously treated as small, requiring |q^| ~ 1^21 ~ qr to he large 
compared with the hadronic scale A. We thus have a three-scale problem characterized by 
the hierarchy 



Large Qi implies that at least some of the transverse parton momenta fcj and fcj must 
be large. The occurrence of partons with large transverse momentum kx can be thought 
of as resulting from the perturbative splitting of partons with low kx, which leads to a 
factorization formula for transverse-momentum dependent parton distributions in terms 
of a hard-scattering kernel and collinear distributions. This significantly adds predictive 
power since collinear distributions depend on fewer variables than kx dependent ones. 

Example graphs for the case of a single-quark distribution are shown in figure 29. 
The description based on such graphs was extensively used for spin effects and azimuthal 
correlations in Drell-Yan production [134-136] and semi-inclusive DIS in [117, 136, 137], 
building on the seminal work in [48, 118]. 

This description carries over to the case of two-parton distributions and is discussed 
in section 5.1. In the subsequent sections we investigate a competing mechanism for the 
generation of high transverse momentum, in which the two partons with momentum frac- 
tions xi and X2 originate from the perturbative splitting of a single parton. We will see 
that this mechanism has profound consequences for the theoretical description of multiple 
interactions. 




(4.49) 



A < < Q ■ 
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a 



b 



Figure 29. Graphs for a single-parton distribution /(x, fc) at perturbatively large fe. Here and 
in tlie following it is understood that lines emerging from the lower blob have virtualities of order 
A <C The cikonal line in graph b results from the Wilson lines in the definition of f{x, fc), see 
section 3.2.1. 



5.1 Ladder graphs at large y 

The single and double ladder graphs in figure 30 are natural generalizations of the ladder 
graph for a single parton density in figure 29a. In the following we concentrate on these 
ladder graphs, bearing in mind that in a covariant gauge there are further graphs with 
eikonal lines as in figure 29b. Those do not change the conclusions we obtain for the ladder 
graphs. Indeed, they are absent in the axial gauge vA = 0, where the Wilson lines in the 
definition of parton distributions reduce to unity (apart from pieces at infinity, as discussed 
at the end of section 3.2.1). 

When interpreting the graphs in figures 29 and 30 it is important to bear in mind that 
they represent a separation of dynamics at different scales, with lines attached to the lower 
blob having virtualities of order A, whereas propagators in the upper part of the graphs are 
for virtualities of the order of the large transverse momentum qT- An important feature of 
figure 30 is that no hard gluons are exchanged between the parton lines that have different 
momentum fractions xi and X2 at the top of the graphs. The requirement that both I — 
and l+^r have small virtualities forces \r\ to be of order A, which translates into interparton 
distances \y\ of hadronic size in the Fourier transformed distributions F{xi,ki,y). 




Figure 30. Ladder graphs for the region of small r with large fci (a) or large fci and fe2 (b). The 
power behavior refers to F(xi, ki, r) and is discussed in the text. Here and in the following we omit 
the dashed line that indicates the final-state cut as in figure 29. 
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5.1.1 Power behavior 

Let us first investigate the general power behavior associated with ladder graphs. In the 
following we refer to qt as the hard scale (compared with A), keeping in mind that qx is 
still much smaller than Q. 

We proceed in a similar way as in section 2.4. In particular, we use the modified 
parton correlation functions which contain a factor for each quark or antiquark 

of momentum / and in which pairs of quark fields are contracted with a Dirac matrix Va 
from (2.81). For the transition from k to m partons in the t channel, we correspondingly 
use hard-scattering coefficients V^_>^ that include a factor \//+ for each incoming quark or 
antiquark and a factor l/v7+ for each outgoing one. Spinor indices in V' are contracted 
with an appropriate matrix ^T"*", ^'J~^J5, ^io'~^-'"f5 for outgoing lines and with ^7", ^757", 
^i<J~^j5 for incoming ones. V is invariant under a boost along z and thus can only depend 
on the scale qt but not on Q (cf. the corresponding argument for <!>' in section 2.4). One 
thus obtains 

~ q'^-'-'^ , (5.2) 

as one can easily check for the example graphs below. Note that compared with the 
hard-scattering amplitudes in (2.127) we now have 3m instead of m because V includes 
the propagators of the outgoing partons, as is appropriate for the calculation of parton 
distributions. 

The power behavior for the single-ladder graph in figure 30a can be obtained from 

F(xi,fci,r)|j.g .^^^=p+kfk^ j dr" dk^ dk^ d\ ^4 

p+A;+A;+ j dk^ dlf j dr" dk^ dq d^l^ $'4 . (5.3) 

The factor p+ and the integrations over minus-momenta come from the definition of F, 
whereas kfk2 compensates the corresponding factors in V' and It is understood that 
V' includes a 6 function for each parton line going across the final-state cut. This does 
not affect the power counting, since one could first consider the hard-scattering amplitude 
without cut and then take the appropriate discontinuity in the s channel. The momenta 
k2 ± ^r and li ± ^r attach to the parton distribution at the bottom of the graph and hence 
have virtualities of order A, whereas ki it emerges from the hard scattering and hence 
has virtuality of order qt- As a result, the momentum components fc^,Z^,r~ ~ hP'/p^ 
and |Zi| ~ A are small and can be neglected in the hard-scattering kernel V . We used 
this when rearranging the order of integrations in the second step. By contrast, the large 
components k^ ~ Qt/p^ ^^"i ^i"' ^2 ^ ~ kept in V' . For the power behavior 

we obtain 

F{xi, ki, r) I ~ Q, X X x (A^)^ x A'^" = q, x (5.4) 

A qj, 

with (2.126) and (5.2). We recognize the l/q^ behavior that is characteristic of the splitting 
of one parton (the incoming quark) into two partons (the outgoing quark and the gluon). 
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The power behavior for the double-ladder graph in figure 30b is obtained by the same type 
of analysis: 



F{xi,ki,r) 



fig. 30b 



fk^ j dk^ dlf V2^2 j dk2 dl^ ^2-^2 j dr~ dq dq d^h d^h $4 



1 



ai X {q't X q^y X (A")^ X A-^" = x ^ . (5.5) 

(It 



If all transverse momenta are small, the distribution F{xi,ki,r) scales of course like A~^. 

Let us now see how the power behavior of the two-parton distributions translates into 
the power behavior of the cross section 



s^dcF 



Wi=i dxi dxi d^Qi 



oc 



figs. 31a, b 



2 -| r 2 „ 

sai{xiXis) JJ / d^ki d^ki 5^'^\qi - ki - ki 

i=l i=l 



X d^r F{xi,ki,r) F{xi,ki,-r) , 



(5.6) 



where we have omitted numerical factors as well as labels for parton species, spin and 
color. We have multiplied the cross section with for convenience, since this gives factors 
sa{xiXis) of order 1 on the r.h.s. 

To have both large Qi and ^2 requires at lowest order in either a single-ladder 
graph in the distribution for each colliding proton, or a double-ladder graph in one of the 
distributions with no hard gluons in the other, as shown in figure 31. In both cases one has 
F{xi,ki,r) F{xi, ki, —r) ~ a1/{h!^qj^) for the product of distributions, and the integration 
volume d'^ki d^ki 6^'^\qi — ki — ki) is of order A^ since k^ and thus ki — are restricted 
to be of size A. Similarly, one finds d^k2d^k2 5(2) (92 - ^2 - ^2) ~ in both cases, so 
that the overall power behavior is 



s^da 



Y\i=i dxi dxi d'^qi 



ai X (A2)3 X 



figs. 31a, b 
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(5.7) 



By similar arguments one finds that the power behavior remains the same at higher order 
in as, when one can have more than two ladders in the graphs for the cross section. Any 
decrease by a factor A'^/q^ in the product F{xi, ki,r) F{xi, ki, —r) is compensated by an 
increase from A^ to q^^ in the integration volume over the transverse parton momenta fcj 
or ki. 

5.1.2 Factorization formulae 

Let us now investigate the structure of the factorization formulae for ladder graphs. Still 
omitting spin and color indices for the moment, we write our result (5.3) as 

I. ^| _ /"^P"^")^ dl+ 



ki J dk^ ^2->2 
P'^^iH j dr^ dq dk2 d^li^'^{li,k2,r) 



(5.8) 
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Figure 31. Ladder graphs contributing to tlic cross section for boson pair production at large qx- 

The lower limit on the if integration reflects that the gluon with momentum li — ki 
crosses the final-state cut and hence cannot have negative plus-momentum. Up to a nu- 
merical factor, the second expression in square brackets is just the transverse-momentum 
dependent double-parton distribution F{ui,X2,li,k2,r) with = if /p^ ■ The first factor 
in square brackets is invariant under boosts along z and can thus only depend on fci and 
if /kf = Ui/xi- We can write it as a numerical factor times k^"^ P(^Ui/xi,kf) , where P 
is dimensionless. If P is a scalar it depends on fc^ only via kf and only because there are 
other dimensionful variables C and /i, which we have not displayed for ease of writing. For 
certain parton polarizations, P is a tensor with transverse indices and can hence depend 
on the components of fc^, as discussed below. 

We can finally Fourier transform (5.8) from r to y and then have 

F(xi,X2,fci,fc2,y)|fi 30a = -r2 / P[—,ki] (fliF{u,X2,h,k2,y), (5.9) 

TTK^ Jxi ^1 / J 

where the factor l/vr has been chosen for convenience. Inserting this and its analog for 
F(xi, X2,ki,k2,y) in the cross section formula (2.35), we have 



da 



Y\i=i dxi dxi d'^Qi 



1 f - 1 

= o-i{xiXis) a2{x2X2s) / d^kid'^k2 

fig. 31a ^ J 



vr(qi - k^Y 7r(q2 - k^f 



^--dni^/xi X ['-^^du2p(X2 



XI Ul \Ui J U2 \U2 



^ j d^y j d%FiuuX2ji,k2,y) j d%F{ xi,U2,ki,l2,y) , (5.10) 

where we have used the 6 function constraints in (2.35) to eliminate fci and ^2- We can 
now approximate — fc^ ~ and ^2 ~ ^2 ~ 92' after which the integrations over ki 
and k2 only concern the double-parton distributions, which are then integrated over both 
transverse-momentum arguments. The result is 

da 1 11 /-^-^^ 



ni=i dxi dxi d^Qi 



-Cri(xiXiS)cJ2(X2X2S) it ^ / P\ ' 9l 

fig. 31a C Trqi vrq^ ui \ui 



x/ ^P(^,q2) d^yF{ui,X2,y)F{xi,U2,y). (5.11) 

Jx2 ""2 ^U2 ) J 
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Only collinear two-parton distributions appear on the r.h.s., so that the relation (5.9) is 
only needed in the form 



d^fc2i^(a;i,X2,fci,fc2,2/)|fi = -72 / P(—,ki\F{u,X2,y). (5.12) 



Repeating the preceding arguments for the double-ladder graph, one obtains 

'■1—X2 



F{xi,X2,ki,k2,y)\^^^^ = -^-^ [ —p( — ,ki 

X ^p(^,k2)Fiu,,U2,y) (5.13) 

A, U2 VU2 / 



from the result (5.5) and 



da 



Y\i=i dxi dxi d'^Qi 
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x/ ^P^^,q^ d^yF{ui,U2,y)F{xi,X2,y) (5.14) 

for the contribution of figure 31b to the cross section, where again only collinear two-parton 
distributions appear. 

The analog of (5.12) for single-parton distribution reads 

/(-.'=)l..,« = i./'t^(;''')^<''*- <^-'^' 

At order Og the kernel P{x/u) is just the usual DGLAP splitting function, up to terms 
proportional to 6{1 — x/u) which will be discussed shortly. To see this, let us consider the 
kx integrated parton density defined with a naive cutoff, 

f{x;f,)''^^' I d^ke{fi'-k^)fix,k^). (5.16) 

Using that the transverse-momentum dependent density depends only on the square of k 
we then have 

/(a;;/^) ''=™ TTi? f{x,k'^ = ix^) , (5.17) 

and comparing with the DGLAP equation for the l.h.s. we can identify the kernel in (5.15) 
as the familiar splitting function. 

The preceding argument is oversimplified in two respects. Firstly, the calculation of 
f{x,k) for large k only involves real graphs like those in figure 29 at leading order in a^, 
because to obtain a parton with large k one needs a recoiling parton in the final state. 
(Higher-order graphs can include virtual loops, so that our argument cannot be applied 
any more.) By contrast, the evolution equation for the collinear parton density f{x), which 
is integrated over all k, involves both real and virtual graphs at order as- The latter give 
contributions proportional to 5(1 — x/u) to the DGLAP splitting kernels, which are absent 



- 103 - 



from P{x/u) in (5.15). Secondly, the distribution /(x, k) must be defined with non-hghthke 
Wilson lines as discussed in section 3, which leads to a dependence on the parameter 
defined in (3.58). Since the collinear distribution f{x) has no such dependence, it is the 
kernel in (5.15) that must depend on The explicit calculation in [117] shows that the C 
dependence of P[x/u) comes with a factor 5{1 — x/u), which in the light of our discussion 
in section 3.2.3 is plausible if one observes that the point x = u m. (5.15) corresponds to 
infinite negative gluon rapidity in figure 29. 

Let us mention that there is an analog of (5.15) for the distribution f{x,z) at small 
transverse distance z. The corresponding hard-scattering kernel differs again from P{x/u) 
by terms proportional to (5(1 — x/u). This is because f{x,z) is given by an integral 
J d^ke^^^ f{x,k) over all k, so that already at order virtual graphs appear in addi- 
tion to real ones. 

We now turn our attention to the role of color in two-parton distributions at high 
transverse momentum, which we have glossed over so far. Since the graphs we are discussing 
do not connect parton lines with different momentum fractions xi and X2, the color coupling 
of the distributions on the left and on the right of the factorization formulae (5.12) and 
(5.13) are the same. For distributions ^F{xi, ki,y) in the color singlet channel, the kernels 
P in (5.12) and (5.13) coincide with the one in (5.15) at least at leading order in Us, 
since the relevant hard-scattering graphs to be calculated are identical. The leading-order 
kernels for the color octet distributions ^F{xi,ki,y) differ by an overall color factor from 
those for ^F{xi, ki, y), which is the subject of the next section. Note that, unlike their color 
singlet analogs, the collinear color octet distributions appearing on the r.h.s. of (5.12) and 
(5.13) depend on as discussed in section 3.5. Given our discussion of the cancellation or 
non-cancellation of soft contributions between real and virtual graphs in section 3.5, we 
expect that the kernels for the position space distributions ^F{xi, Zi,y) and ^F{xi, Zi,y) 
at small Zi will differ by more than an overall color factor already at leading order. A 
systematic investigation of this is left to future work. 

The power counting in section 5.1.1 and the discussion in the present section do not 
depend on whether the parton lines in the ladder graphs are quarks or gluons. As is 
well-known for single parton distributions, a quark with high transverse momentum can 
originate from a gluon with low transverse momentum and vice versa. The corresponding 
elementary ladder graphs are shown in figure 32 below. 

We now discuss the spin structure of the ladder graphs. As we have seen in sec- 
tions 2.2.1 and 2.2.2, there are three polarization combinations for each quark or gluon in 
a two-parton distribution, which we can choose as unpolarized {q,g), longitudinally po- 
larized {Aq, Ag) and transversely {5q) or linearly {6g) polarized. The possible transitions 
between these combinations in the factorization formulae (5.12) and (5.13) are restricted 
by symmetries. Transverse quark polarization 5q is described by a chiral-odd operator 
qia~^^^^ q and the chirality of light quarks is conserved in hard-scattering subprocesses, so 
that the only transitions for transversely polarized quarks are of the form 6q — t- 5q. In 
the longitudinally polarized sector one has all possible transitions between Aq and on 
the left and on the right of (5.12) and (5.13), with transitions to other polarizations being 
forbidden by parity invariance. Likewise, one has all possible transitions between unpo- 
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larized quarks and gluons. In addition, ladder graphs allow the transitions g ^ 6g and 
q ^ 5g from unpolarized collinear distributions to linearly polarized gluons at high k, as 
has been observed in the study of single Higgs production in [83, 84]. Since 5g corresponds 
to a helicity difference of two units between the gluon on the left and the gluon the right 
of the final-state cut, there is no collinear distribution for a single linearly polarized gluon 
in a proton, so that transitions 6g ^ g and 5g ^ q played no role in [83, 84]. However, 
one finds that the corresponding hard-scattering kernels are nonzero and hence allow these 
transitions for two-parton distributions. One thus has all possible transitions between q,g 
and 5g. 

For distributions that involve polarizations 6q or 6g the kernels P in (5.12) and (5.13) 
are tensors with transverse Lorentz indices, constructed from 5-^' and from the large trans- 
verse momentum ki or k2 in the ladder. Explicit calculation at order ag shows that the 
kernel Pg^^^ for the transition from F^^^ to F^^ ^ (with arbitrary a) is proportional to 
As a result, ladder graphs do not generate the distributions and g^^ in the decom- 

position (4.12) of F^^g^{xi, ki, y), given that they come with tensors that are absent in the 
collinear distributions Fg^^ g^{xi,y) according to (4.14). This adds to the predictive power of 
the perturbative mechanism at large transverse momenta. For transitions involving linear 
gluon polarization we find kernels 

P&g' « r^'''' , PsU « PsU « P^i « P^Sg ^ 2fclfcf - , (5.18) 

where the transition from Fg^ ^ to F^^ ^ is described by P^g/g , the transition from Fg ^ to 
Psg a by P^lg etc., and where fcj is fci or fc2- The second tensor in (5.18) is symmetric and 
traceless and describes two units of orbital angular momentum, which compensates the 
mismatch of helicities in the transitions g — t- 5g, q — t- 6g, 6g ^ g and 6g — t- q. For later use 
we note that terms involving this tensor vanish by rotation invariance when (5.9) or (5.13) 
is integrated over fci and fc2. 

The representation of the cross section derived in this section is based on a two-step 
procedure. In the first step we have used factorization to separate the annihilation processes 
qq ^ V into vector bosons V with mass or virtuality of order Q from transverse-momentum 
dependent two-parton distributions, in which the largest scale is qx- In a second step, we 
have used factorization to compute these distributions in terms of hard-scattering processes 
at scale qT and distributions that reflect physics at a hadronic scale A, where it turned out 
that in the cross section we only need the latter distributions integrated over fci and fc2. 

An alternative procedure is to first use factorization to represent the graphs in figure 31 
as the product of collinear two-parton distributions and inclusive hard-scattering processes 
qq — )■ V + X, where at lowest order in Og the unobserved system X consists of just one 
gluon. In a second step one can then simplify the corresponding hard-scattering kernels by 
taking the limit ^ Q we are interested in. The relation between these two procedures 
has been studied in detail for single Drell-Yan production or for semi-inclusive deep inelastic 
scattering in [117] and [134-137]. An important property of the procedure using transverse- 
momentum dependent distributions in a first step is that it permits the resummation of 
Sudakov logarithms with the method of Collins, Soper and Sterman [118]. 
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Figure 32. Elementary ladder graphs for transitions between quarks and gluons. 
5.1.3 Color factors and quark-gluon transitions 

Let us now compute the color factors for the transitions from collinear to transverse- 
momentum dependent two-parton distributions, restricting ourselves to the leading order 
in as- For this it is sufficient to consider ladder graphs for the hard scattering, because 
graphs with eikonal lines as in figure 29b can be eliminated by choosing the gauge vA = 0. 
Multiplying the color structure of the ladder graphs in figure 32 with the appropriate color 
matrix for the incoming partons, one obtains the following color transitions: 
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^ ^g, 'g ■■ 4, tt^ i]^, = -i ir'^' + i d-'^' . (5.19) 

Reversing the fermion lines in the ladder graphs changes the order of multiplication for the 
t matrices on the l.h.s. of the above relations. This leads to a change of sign on the r.h.s. 
for the transitions "^g — )■ \ and \ — )• ^g, whereas transitions between \ and ^g or between 
singlets are unchanged. This implies that the difference of distributions for q and q does 
not mix with gluons in the singlet or the symmetric octet channel but does mix with gluons 
coupled to an antisymmetric octet. 

In the cases where there is mixing, the relation (5.12) for a double-parton distribution 
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at large ki becomes a matrix equation, which can be written as 



/ 



(fk2 F\xi,X2,ki,k2,y) 



1 



E 



'""^p"(^i)f''{«i....«) 



(5.20) 



where the arguments of F and P are as in (5.12). For polarizations 5q and 5g the distribu- 
tions and kernels carry tensor indices, which were discussed in the previous section and will 
be omitted here. The matrix structure in (5.20) can be generalized to the relation (5.13) 
for double- ladder graphs, but the resulting expressions are rather cluttered with indices 
and will not be given here. 

In the color singlet sector one has 



IrpJ 
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^F&n,a 
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IpJJ' 
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CpPg, NPgg NPgSg 

yCpP&g, NPSgg NPSgSgJ 



(5.21) 



where ni;' is the number of quark flavors and where the second parton a may be an un- 
polarized or polarized quark, antiquark or gluon.^^ To obtain the color factors for the 
off-diagonal elements one must take into account the prefactors in the definitions (2.103), 
(2.118) and (2.123). In the upper left 2x2 submatrix of ^P we recognize the structure of 
the mixing matrix in the usual DGLAP equations. Mixing in the symmetric octet sector 
involves the vectors 



[Eg ['F,,a + %A 



9,0. 



\ 



5g,a 
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9,0. 



^Fx 



(5.22) 



If a indicates a gluon, one should replace ^Fq ^ + ^Fq^a by Fq^a + %,a and by 

Ap Aj 



q^a — Fq^a- The Splitting matrices now read 
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(5.23) 



and 
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(5.24) 



We note that in [47] the possibility of transitions between q, g and 5g was overlooked, and only the 
mixing between q and g was considered. 
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There is no mixing for the combinations [^Pq,a — ^Pg,g] > iior for the difference of distri- 
butions for two quark flavors, nor for matrix elements where the quark flavors differ on 
the two sides of the final-state cut as in figure 7. In these cases the behavior at large ki 
is described by (5.12) with P replaced by CpPqq for color singlet and by —^Pqq for color 
octet combinations. 

Experience with the usual parton densities tells us that gluons quickly dominate over 
sea quarks as one goes to momentum fractions below 0.1 (except possibly if one considers 
very low factorization scales). One therefore expects that for typical values of xi and X2 
at LHC, two-parton distributions at high transverse momentum are dominated by those 
combinations that can originate from gluons in (5.20). 

Comparing (5.21) with (5.23) and (5.24) we find that the color factors are always 
smaller in the octet channels than in the singlet channel, with the biggest suppression 
occurring for Pqq. In the large- limit the singlet matrix ^P has one eigenvalue y^qq ^^'^ 
two eigenvalues with color factors N for the submatrix in the g, 6g sector. Both ^P and 
'^P have the same two eigenvalues, but with color factors ^ instead of N, and another 
eigenvalue of order 1. One can hence expect a dominance of color singlet distributions for 
sufficiently large transverse momentum, which would significantly simplify the theoretical 
description and the phenomenology of multiple interactions. How strong the suppression 
of color octet channels is in given kinematics should, however, be studied quantitatively 
before drawing strong conclusions. 

We have also calculated the color factors for higher color representations of gluon 
distributions, restricting ourselves to = 3 as we did in (2.121). Mixing with quark 
distributions is of course not possible in this case. We find that the color factors for 
decuplet and antidecuplet distributions are zero, so that ladder graphs do not admit these 
color combinations, at least not at leading order in Og- For the 27 representation we obtain 

where a = g or 6g. The color factor is equal to —1 and thus smaller in magnitude than the 
factors N or ^A^ we have in the singlet and octet sectors, respectively. 

The color factors we have obtained agree with those given in [138], provided that one 
restores a missing factor ^JljN in the expression of Pg/ in eq. (54b) of that paper. Our 
color factors for transitions in the gluon sector are also in agreement with eq. (A. 6) in [140]. 

The splitting matrices for longitudinally polarized quarks and gluons in (5.20) are 
obtained from the upper left submatrices for unpolarized quarks and gluons in (5.21) to 
(5.24) by changing the kernels but keeping the color factors. Likewise, the splitting kernel 
P^qiq for transverse quark polarization comes with the same color factors as Pqq. In both 
cases, it is again the color singlet sector that has the largest color factors and will therefore 
be enhanced at high transverse parton momentum. 

Let us finally consider ladder graphs for the quark-antiquark interference distributions 
/ represented in figure 5c. The color independent part of the splitting kernel is different 



^The projector Pq in [138] appears only for SUiN) with > 3 [139]. 
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Figure 33. Graphs for a quark- antiquark distribution tliat involve the perturbative splitting of 
one into two partons on both sides of the final-state cut. The power behavior refers to F{xi, ki, r) 
and is discussed in the text. 

for distributions / and F because of their different spin structure (in one case a gluon is 
exchanged between a quark and an antiquark line and in the other case between two quark 
lines), but we shall not pursue this issue further here. We can, however, easily determine 
the color structure of the graphs. For definiteness, consider the exchange of a gluon between 
the two lines with momentum fraction xi and color indices j, j' . The color decomposition 
that remains invariant under this exchange is the one in (2.116), since 

fji i'j'v i^ik' 5vk + 5ik Si'k') = f^k' ifk + *ifc *i'fc' = 2N ^^'^ ^ ^^'^'^ 

tji t°jiii {Sik' Si'k — ^ik^i'k') = t'jk' t'j'k ~ ^°jk t'j'k' = 2W ^-^'^ ~ ^^^^i'^'^ (5.26) 

For = 3 we thus find a color factor ^ if the quarks are coupled to a sextet and — § if they 
are coupled to an antitriplet. Both factors are smaller than Ci? = | in the color singlet 
channel. 

5.2 Parton splitting at high transverse parton momenta 

We now turn to another mechanism that generates large transverse momenta in multiparton 
distributions: the perturbative splitting of one parton into two partons, both of which 
subsequently take part in hard-scattering processes. This mechanism turns out to be 
enhanced by powers of qt/^ in the cross section. In addition, it leads to conceptual issues 
concerning the very notion of multiparton interactions. In the following sections we derive 
several results about the splitting mechanism, but we will be left with a number of open 
questions for future research. 

5.2.1 Power behavior 

There are a multitude of graphs involving the splitting of one parton into two partons, and 
in order to assess their importance we use power counting as our first guiding principle. 
Simple examples for parton splitting are shown in figure 33. They allow all transverse 
momenta fci, ^2 and r to be large. This leads to large virtualities for ki, ^2 and r, so that 
their minus components are large as well, 

fci ~ ^2 ~ ~ (Zt ) fe^' ~ /cg^ ~ r~ ~ Qt/p^ ■ (5-27) 
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At lowest order in as one has disconnected graphs as in figure 33a, which leads to the 
kinematic constraint 

fci + fcs ~ A , /cf + A;^ ~ AVp"^ • (5.28) 

We will see shortly that this constraint plays a special role when two-parton distribu- 
tions are combined to calculate a cross section. The power behavior of graph 33a can be 
determined using the same method as in section 5.1.1, and we have 

F(x„fci,r)|gg 33^=p+A;+/c+ j dr^ dk^ M^2\l~=_k- ,k,=-k, j dk^ <^'^{k^ + k^) 

~ Os X gf, X {q:^^f X X A~^ = x ^ , (5.29) 



Qt 



where the power counting of integration volumes is dr~ ~ dk^ ~ Qt/p^ '^^2' ~ A^/p"*" 
in order to fulfill the constraint (5.28). Notice that in V' one should set A;^ = —k^ because 
the difference of fc^ and —k^ is negligible compared with k^ . Likewise, one should set 
^2 equal to —k^ in V' . The factor / dki^ ^2{ki + /C2) in (5.29) is proportional to the 
transverse-momentum dependent distribution of a gluon in the proton. 

Starting at order one has connected graphs as in figure 33b. The restrictions (5.28) 
are then lifted, and we obtain a power behavior 

F{xi, ki,r)\^^ = p^kfk^ j dl^ dr' dk^ dk2 Vi^^ j dr d^l ^>2 

~ X 4 X qZ,^^ X A^ X A"^ = X 4- • (5.30) 

Qt 

The graphs just discussed describe the transition from two to four partons in the t 
channel. Let us compare them with transitions starting from three or four partons. The 
corresponding graphs admit a variety of topologies, and we shall not give a comprehensive 
treatment here. Important examples are shown in figures 34 and 35. They are subject to 
different kinematic restrictions: 

• graphs 34a and 35a require -|- fc2| ~ A and are thus analogous to 33a, 

• in analogy to 33b, graphs 34b and 35b produce partons with unconstrained transverse 
momenta, 

• in graphs 34d and 35d one must have |fci -|- ~ A since the rightmost quark line is 
disconnected, 

• graphs 34c and 35c are subject to both constraints \ki + k2\ ~ A and |fci -|- \r\ ~ A. 

The power behavior of the resulting two-parton distributions can be obtained by the same 
methods as previously and is given in the figures. We see that within a given kinematic 
group, the graph with the smallest number of partons initiating the hard scattering is 
dominant by power counting, i.e. two partons in cases a and b, and three partons in cases 
c and d. The graph with the leading power behavior also has the lowest power of Ug- 
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c 



d 



Figure 34. Example graphs for the transition of three to four partons in the t channel and the 
corresponding power behavior of F[xi,ki,r). They involve the splitting of one into two partons 
only on one side of the final-state cut. 




Figure 35. As figure 34, but for the transition of four to four partons in the t channel. 

We note that for graphs starting with four partons there are topologies leading to different 
kinematic constraints, such as those in figure 36. We shall not discuss these in the following. 

So far we have assumed that the transverse momenta fci, are all large. How- 

ever, the graphs we have discussed remain under perturbative control in more restricted 
kinematics as well. In graph 33a for instance, we need large transverse momenta for the 
four upper parton lines, i.e. large ki it (recall that k2 ~ — fci for this graph). This 
allows either r or k\ to be small, as long as the other is large. Both configurations will be 
important in our further discussion. The power behavior we have derived above remains 
unchanged in those kinematic regions, as we shall see explicitly in section 5.2.2. 

Cross section. Let us now see how the different contributions to multiparton distribu- 
tions enter the cross section for large Qi and Taking the lowest-order parton splitting 
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a 



b 



Figure 36. Ladder graphs with tlie kinematic constraint |fci — A;2| ~ A and the resulting power 
behavior of F{xi,ki,r). 



contribution of figure 33a for both protons, we arrive at the graph in figure 37a. Both 
ki + k2 and fci + k2 are restricted to be of order A in this case, which imphes 

9l + 92 = ^1 + ^1 + ^2 + ^2 ~ ^ • (5.31) 

In other words, the produced bosons must be almost back to back in transverse momen- 
tum. To determine the power behavior of the cross section, we note that the integration 
element cPki cPki 6^'^\qi — ki — ki) scales like since ki can be freely chosen of size qt- 
Once this choice is made, k2 can only differ from —ki by an amount of order A, so that 
(Pk2 (Pk2 S^"^^ {q2 — k2 — ^2) scales like A^. The corresponding constraint on ki + ^2 is then 
automatically fulfilled by virtue of (5.31). With cPr ~ and the scaling behavior (5.29) 
the cross section formula (5.6) then gives 



„2 



da 



l<Ji+'J2l~A / 1 \2 1 



X (4)2 X A2 X (-_) = X — . (5.32) 

fig. 37a ^■'^ 



Going one order higher in a^, one has graphs as in figure 37b with a connected two-to-four 
parton transition on one side. The constraint (5.31) is then lifted, and Qi and ^2 can be 
chosen independently. The integration elements in the cross section formula scale as in the 
previous case, but due to the stronger falloff in qj- in (5.30), one now has 

' ~ X {q'j.f X A^ X ^ X = x \ . (5.33) 

fig. 37b Qt ^ It It 



Yli^i dxi dxi d^qi 



At yet higher order in as one obtains the same power behavior if both two-parton distri- 
butions contain a connected two-to-four parton transition: the extra factor h? /q^ from 
the two-parton distribution is compensated by an increase from A^ to q^ in the loop phase 
space, since both ki and k2 can then be chosen independently of order qT- 

We note that both (5.32) and (5.33) contribute at the same power of ^/q^ if one 
integrates the cross section over Qi and ^2 ™ ^ region of size qx- This is because the con- 
tribution (5.32) has a restricted phase space of order d'^Qi d'^q2 ~ h?q^. In the differential 
cross section, however, the contribution (5.32) gives a peak in the distribution of -|- q2, 
which is enhanced not only by a power of as but also by q^j . 

There are more contributions to the cross section with the same power behavior as 
the one we have just encountered. We recall from our discussion in section 2.4 that in 
the differential cross section, double parton scattering has the same power behavior as 
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a 



b 



Figure 37. Graphs for the cross section with two-to-four parton transitions in the t channel for 
both colliding protons. Graph a only contributes when \q-^ + ^ 




Figure 38. Graphs for the cross section with two-to-three parton transitions in the t channel for 
both protons. 




a b 



Figure 39. Graphs for the cross section with a single hard scattering. 

the interference of two hard scatters in the amplitude with a single hard scatter in the 
conjugate amplitude (see graph 9a). If the two partons initiating the two hard scatters in 
the amplitude come from the splitting of a single parton, we have graphs like in figure 38.^'' 
For graph 38a one finds the same scaling behavior (5.32) as for graph 37a, and for graph 

^^The single hard scattering to the right of the final-state cut proceeds through a loop in our example, 
because gluons have no direct coupling with electroweak gauge bosons. Other processes, like the production 
of two dijets, can proceed already at tree level. The powers of as in our example are thus not representative 
of the generic case, whereas powers of A/qt are. 



- 113 - 



38b one finds the same scaling behavior (5.33) as for graph 37b. 

The same power behavior is again found for the case where the two gauge bosons 
are produced in a single hard scatter, both in the amplitude and its conjugate. The 
corresponding cross section formula can be found in (2.42). For incoming gluons the hard 
scattering proceeds through a loop as on the r.h.s. of graphs 38 a and b, whereas for 
incoming quarks or antiquarks one has graphs as those in figure 39. If |q^ + ~ A 
then one needs no parton exchange of virtuality qx in any of the parton densities, and one 
immediately finds the same power behavior as in (5.32). If Iq^^ + ~ 9t then at lowest 
order in Og one parton distribution has large transverse momentum generated by a ladder 
graph as shown in fig. 39b. According to (5.15) one has f{x, k) ~ Os/Qj^ for |fc| ~ qx, and 
for the cross section one obtains the same power behavior as in (5.33). 

Adapting our discussion at the end of section 5.1.2 we see that the graph in figure 39b 
can be calculated either in terms of transverse-momentum dependent parton densities, one 
of which involves a ladder graph, or in terms of collinear parton distributions and the 
parton-level process qq — )• V1V2 + g, where Vi and V2 denote the produced vector bosons. 
The result is the same in both cases. 

In a similar way, figures 37 and 38 can be interpreted as graphs for two-boson produc- 
tion by a single hard scattering process at one-loop level, namely by gg — ?• Vi V2 for graphs a 
and gg — t- ViV2 + g for graphs b. The quark lines in each loop are then typically off-shell by 
order Q, which is the hard scale set by the final state. Note that this differs from the case 
when one interprets the same graphs as representing double hard scattering (figure 37) or 
the interference of double and single hard scattering (figure 38). In that case the quark 
lines in the loops (except for those on the r.h.s. of graphs 38a and b) are understood to 
have typical virtualities of order qx, which allows one to treat them as incoming on-shell 
partons in the tree- level subprocesses qq ^ V, whose large scale is Q. Detailed inspection 
of the quark loops in figure 37 shows that they receive contributions with the same scaling 
behavior in qr/Q from the two regions where all quark virtualities are either of order q-r 
or of order Q}^ One thus obtains the same power behavior for the graphs in each of the 
two interpretations just discussed. The interpretation in terms of a single hard-scattering 
process producing two gauge bosons for graphs 37a, 38a and 39a and two gauge bosons 
plus a gluon for graphs 37b, 38b and 39b makes it clear that each group of graphs has the 
same scaling behavior, respectively given by (5.32) and (5.33). 

In section 5.1.1 we found that ladder graphs as in figure 31 contribute to the scaled 
cross section with a power A^/g^, with no distinction between the cases where qi + ^2 is 
of order A or qx- This means that the contributions of figures 37, 38 and 39 are enhanced 
over the ladder graphs by q^/ ^ for Iqi + ^2! ~ ?T and by gf^/A^ for |q^ + 92! ~ ^- As 
we already discussed in section 2.5, one can however expect that at small xi and X2 this 
enhancement is counteracted by the stronger rise of the ladder graphs with decreasing par- 
ton momentum fractions, since the ladder contributions involve two-parton distributions, 
whereas the graphs in figures 37, 38 and 39 depend on single-parton distributions. Whether 

^*By contrast, there is no kinematic region where aU quark hnes in the loops on the r.h.s. of graphs 38a 
or b are oflF shell by order qr- This is easily seen by analyzing the flow of large plus and minus momenta. 
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a 



b 



c 



Figure 40. Contribution of tlie perturbative transition from two to four quarks in one proton, for 
the region of small r. 




a b c 



Figure 41. As figure 40 but for the region of perturbatively large r. 

this small-x enhancement is more important than powers of q^/A'^ cannot be determined 
on generic grounds, so that one will want to keep both types of contribution in a flexible 
theoretical description. 

With this in mind, we now turn our attention to the graphs in figures 40 and 41, which 
involve parton splitting and thus single-parton distributions for one proton but a two- 
parton distribution for the other. In the graphs of figure 40 the two-parton distributions 
in one proton force r to be of order A, whereas in figure 41 an additional gluon exchanged 
between partons with momentum fraction xi and X2 allows r to be of order qt- 

The corresponding power behavior is readily obtained from our results for the relevant 
parton distributions (given in figures 30, 33 and 35) and the available loop phase space in 
each graph. We find 



c2 



da 



|<3l+<32l~A 



Y\i=idxi dxi d^Qi 
s^da 



~ a., X — 



fig. 40a '?r 

A2 



~ "2 X :V (5-34) 



figs. 40b,c 



and the analogous scaling behavior with an extra power of ag for the graphs in figure 41. 
For + 92! ~ Qt, we thus find the same behavior as for the ladder graphs in figure 31, 
which involve however two two-parton distribution in the cross section and therefore have 
a stronger small-x enhancement. In the region Iqi + 92! ~ A we have an extra power of 
q^/A"^, as in the other parton splitting contributions discussed so far. 
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a 



b 



Figure 42. Graphs with three partons in the t channel for one proton (a) or for both (b). 



In figure 34 we have graphs initiated by proton matrix elements with three partons in 
the t channel. Examples for their contribution to the cross sections are given in figure 42. 
They behave as 

191+921 ~ A 



„2 



da 



Y\i=i dxi dxi d'^Qi 
s^da 



fig. 42a 
l'Jl+'J2l~^ I 

~ as X — . (5.35) 

fig. 42b It 



Wi^i dxi dxi d'^Qi 

The contribution from graph 42a is thus suppressed compared with the one from graph 37a, 
although only by a power oi A/qt, which corresponds to the loss of one power A/qt between 
the splitting graphs initiated by two or three partons in the t channel (cf. figures 33a and 
34c). Likewise, graph 42b is suppressed by A^/g|, compared with graph 37a, thus having 
the same power behavior as graph 40a, but lacking the small- j; enhancement of the latter. 
An analogous situation is found for graphs that are like in figure 42 but have an extra gluon 
across the final-state cut (constructed e.g. from graphs 33b or 34d) and thus contribute 
to the region + ~ It- To the extent that one power of A/qt is a small enough 
suppression parameter and that the small-x enhancement of four-parton matrix elements 
is important, one can hence neglect contributions involving three-parton matrix elements. 
Our results for the power behavior of the different contributions are collected in table 1. 

The contributions in the first three rows of the table were recently investigated in 
[98]. It was pointed out in that work that the 2x4 contribution in our table has a further 
enhancement compared with the 4x4 term, which is due to the fact that the latter involves 
the product F{xi,r) F{xi, —r) of two distributions that decrease with r, whereas the former 
involves only one factor F{xi,r) multiplied by the perturbative splitting contribution that 
is approximately r independent for \r\ ~ A. 

To compare our results with those in [98] we note that | Qi + ^2 1 ~ required to 

be in the perturbative domain in that work, whereas we treat it as comparable to a soft 
scale. Our power counting results apply to this case as well as far as the qt behavior is 
concerned, if one understands A as either Iq^^ + ^2! or a generically soft scale, without the 
ability to distinguish between them. What is important for our results is the hierarchy 
A <C gr ^ Q, which in the notation of [98] reads 6' <^ S <^ Q. The fact that in [98] four-jet 
production rather than the double Drell-Yan process was studied does not prevent us from a 
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partons 


example graphs 


power 


behavior 


(t channel) 








4x4 


31 


1 


1 "J 


4x2 


40, 41 


1/4 




2x2 


37, 38, 39 


1/A2 


1/4 


2x3 


42a 


l/(AgT) 


A/4 


3x3 


42b 


1/4 


AV4 



Table 1. Power behavior to the scaled cross section da j (dx\dx2dx\dx2(P' q\<P q^) from various 
contributions. An entry m x n in the first column means that the cross section involves matrix 
elements with m and n partons in the t channel for the first and the second proton, respectively. 



comparison since, as we pointed out earlier, our power counting results hold independently 
of the particular hard-scattering processes. We agree with [98] that the 2x4 and the 
4x4 contributions to the cross section respectively behave like 1/4 ^-i^d 1/4) ^'^d that 
the 2x2 contribution does not have a 1/4 behavior but depends logarithmically on q^- 
However, the authors of [98] write that the 2x2 term is comparable to the 2x4 and 4x4 
terms. We emphasize that the 2x2 contribution goes like 1/A^ in the scaled cross section 
and is therefore power enhanced compared with the other two contributions for q^ ^ A. 
This comes out of our power counting analysis and is confirmed by explicit calculation, see 
(5.74) below. What can potentially make the 2x4 and 4x4 terms more important is their 
small-x behavior, as we noted above. 



5.2.2 Splitting in two-parton distributions 

After the general analysis in section 5.2.1 we now investigate splitting contributions to 
two-parton distributions in more detail. We begin with the graph in figure 33a, which 
describes the splitting process g — qq. 

From the color structure of the graph we readily find that it gives rise to color octet 
distributions that are suppressed compared with the color singlet ones by a factor 
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IF - 

0.1,0.2 

The color singlet distributions are given by 

47ra, 1 



(5.36) 



01,02 



(Xj, fcj, T") 



g^gq (Itt)^ 2 



2p+j dr~dk^ dk^ ^ifsih + k2) 



X tr 



(fcl-H7 a 

Ol 1 NO , . / 



(/ci - ir)2 ie {k2 + lry + ie 
ih - ^r)-f g jki + ^r)7 

X 1 r, ^' 



(fc2 



ir)2 - ie ' (fci + ir)2 



(5.37) 
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where a and /3 are polarization indices of the gluon potentials in the correlation function <I>^, 
whose definition follows from (2.96). is already summed over the gluon color indices, 
and the corresponding trace over color matrices has given a factor |. As discussed in 
section 2.2.2, a and /3 are restricted to be transverse at leading-power accuracy. The 
second and third line of (5.37) represent the hard part of the process, where we can neglect 
the difference between the transverse and minus-components of k\ and ^2, see (5.29). We 
introduce 



ki + k2 



(5.38) 



and change integration variables from k-^ and ^2 to k and k . The integration over k 
only concerns the gluon correlation function <I>^, which can be decomposed as [82] 



-5^-> f{{x,K) + - 



4M2 



where M denotes the proton mass. In terms of the operators introduced in (2. 



(5.39) 
we have 



4M2 



h^^{x, k) 



1 



1 

xp'^ 



dz-d^ 



(27r)3 

dz~d^z 
(27r)3 



{p\OgiO,z)\p), 



{p\Oll (0,z)|p> 



(5.40) 



/f is the usual transverse- momentum dependent density of gluons, whereas the gluon 
Boer-Mulders function h^^ describes linearly polarized gluons and is essentially unknown 
at present (see [85-87] for processes where this distribution could be studied). Writing the 
product of propagator denominators in (5.37) as 

1 1 



2xip+ (k 



k 



\r) + ie 2x2P'^(k 



1 



1 \2 

r) 



It 



2x2P+{k + \r) +{k + \r) + ie 2xip+ {k + \r) - {k + ^r)' 



(5.41) 



we see that the integrations over r~ and k~ can conveniently be performed using the 
theorem of residues, after a change of variables to (A; — \r)~ and {k + \t)^ . Performing 
the fermion trace, we finally obtain 



11,12 



9-^qq 

h^^{xi + X2, k) 
X\ + X2 

With the abbreviation u 



as 
47r2 

ll'mm' 



ff{xi+X2,K) If, 

ai,a2 



Xi + X2 

Xl 



Xl 



Xl + X2 



+ 



2K™fv" 



K 



2AP 



Xl + X2 



[k + ir) 



{k + \rY{k-\r) 



2 • 



(5.42) 



1 — u the kernels read 



"q,q 



{u) 



n 



Aq,Aq 



(u) 



^qAq 



-26^^'6"'uu 



(5.43) 
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and 



[U^IJ^' {u)Y^' = iT^^''^'"^' d^'^u + lT^^''^"'' P'^'u-lT^^'^"'"'' 6"'uu, (5.44) 

where j and j' are the indices of the Dirac matrices ia^^'j^ in the definition of the distribu- 
tions. The kernels T and U not hsted in (5.43) or (5.44) are zero. Note that fc — is half 
the difference between the transverse parton momenta ki — and k2 + on the left of 
figure 33a, whereas k + ^r is half the corresponding momentum difference on the right. To 
ensure that the quark lines with momenta ki it and k2 ± in figure 33a are far off-shell 
it is sufficient that one of the transverse momenta r and k is large, as already mentioned 
earlier. This implies that (5.42) describes the large r behavior of Fai^a2{xij ^ii *") foi' small 
ki and k2, as well as its behavior for large k at small r. 

We see that the short-distance splitting process gives rise to a rich spin structure, with 
all chiral even two-parton distributions being nonzero. The relations Fq^q = —F/\q^/\q and 
FAq,q = —Fq^j\q reficct that the perturbative gluon splitting leads to a 100% correlation 
between the helicities of the quark and antiquark: if the quark has positive helicity the 
antiquark has negative one, and vice versa. For values of u around |, the transverse spin 
correlation encoded in Fsg^sg is as large as the unpolarized distribution Fq^q. 

The splitting contributions to other two-parton distributions are obtained in close 
analogy to the case we have just discussed, and in the following we only give the relevant 
starting expressions and results. A reader not interested in the details may skip forward 
to the paragraph after equation (5.60). 

The graph in figure 33a also contributes to the interference distributions /ai,a2; with 
the same ratio ^Iai,a2/ ^^aifii = —l/VN'^ — 1 of octet and singlet distributions as in (5.36). 
The expression for ^Iai,a2 can be obtained from the one in (5.37) by interchanging {k2 — -^r)j 
and {ki + ^'')7 in the fermion trace. The result has the same structure as in (5.42), with 
the kernels Ta^^aj replaced by 

V^!^{u) = -V^i^qiu) = -26'''uu, 
[V5g,5g{u)]"' = S^^' S^^' (u^ + u^) + {S^^6^'^' - 6^^' 6^'^) (n - u) (5.45) 
and the kernels Ua^^a2 replaced by 

Wq,q W = -W^q^^q (u) = T ' (u + U j , 

All other kernels are zero. We see that the splitting contribution to the interference distri- 
butions la^ ^0.2 is generically of the same size as for the distributions Fa^ ^0.2 ■ 

We now turn to the analog of figure 33a for the splitting process q ^ gq- This graph 
(not shown here for brevity) involves propagators for the outgoing gluons and requires a 
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choice of gauge. If we work in the light-cone gauge = ^4+ = with n = (1, 0, 0, — l)/-v/2, 
then the gluon propagator has a numerator 



(5.47) 



and the q ^ gq sphtting contribution to quark-gluon distributions reads 

47ras 

q^gq (27 



Fai,a2{xi,ki,r) = l Cf (xip+) 2p+ / dr'dk^dk^ 

q^gq [iTTy J 



{ki - ^rY + ie 



X tr 



(^2 



(^2 + 



7° + 



(5.48) 



Since j is a transverse index, the numerator factor of the first gluon propagator simplifies 
to —Qaj + f^a{ki — ^r)j/ [ki — ^r)~^. If we work in covariant gauge instead, these two terms 
correspond to the first two terms of the gluon field strength G^^ = d~^A^ — + 0{g) in 
the operator definition of the quark-gluon distribution. An analogous statement holds for 
the second gluon propagator. The expression (5.48) involves the quark correlation function 
for an unpolarized proton, for which one has 



1 _ 1 _ e^^ 

-7-/f(x,^) + -ia^-75^^ 



(5.49) 



to leading-twist accuracy, or equivalently 

dz~d'^z 



M 



(27r)3 

dz~d^z 
(27r)3 



^tXZ p^ — lZK, 



[p\O,i0,z)\p), 
[p\oi^{0,z)\p). 



(5.50) 



The q ^ gq splitting process gives rise to all possible color couplings in the quark-gluon 
distribution in (2.123), with color factors 
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0.2 



Ar2 



^ 0.1,0 

ai,a 
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71 



(5.51) 



Contrary to the case of g ^ qq analyzed above, the splitting mechanism now favors color 
octet distributions over color singlet ones. Evaluating (5.48) we obtain 



^ai ,02 {■^i i^ii'^') 



as 
27r2 



Cf 



11,12 



Xl + X2 



Xi 



Xl + X2 



^mm'^m' hj" {x^ + X^, k) / X^ 

M X1 + X2 ''I'^nxi+xa 



{k+lr)\k-lry' 



(5.52) 
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with 



Ti''{u) = 5''\l + u^)/u., 



(u) = -ie" (1 + 



Tl''^,{u) = -ie''' {l + u) 



and 



[ul^'6^{u)Y = 26^^'6^'^ u/i 



(5.53) 



(5.54) 



All other kernels are zero, in particular F^g^^q is not generated by the splitting mechanism 
at leading order in as- Analogous results can be derived for the splitting q — ?• gq. 

The splitting of one gluon into two gives a contribution to two-gluon distributions, 
which reads 



Ana,, 



{2nY 



I N {xip+) (X2P+) 2p+ / dr- dk^ dk^ ^a^(^i + h) 



fig. 33a 



i-K-f 1 ^ ttn„ ^ 
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.a afi 



(5.55) 



fc2=— fcl 



in the gauge = 0. Evaluating this expression, we obtain a result with the same structure 
as for g ^ qq in (5.42), 
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(5.56) 



with 



T^^' (li) = 2(5''' (n/u + u/n + nn) ^ 



rilA,(^) = 25"- (2 - uu) 



IV , 
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(5.57) 
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and 
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(5.58) 



The kernels T^g,g, Tsg^g, Usg^g and Usg^Ag are respectively obtained from Tg^^g, Tg^gg, Ug^sg 
and U/\g^Sg by interchanging u ■v^ u and the appropriate indices. The remaining kernels are 
zero. For the different color combinations we find 
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(5.59) 



where as in the case q ^ gq color octet distributions are enhanced over color singlet ones. 
The factors for the higher color representations in the case = 3 are 
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(5.60) 



g^aa 



The 27 representation is hence even more strongly enhanced than the two color octet 
combinations. Decuplet and antidecuplet distributions are not generated by perturbative 
splitting at lowest order. We recall that this was also the case for the ladder graphs 
discussed in section 5.1.3. 

We see that the perturbative splitting mechanism gives rise to a multitude of two- 
parton distributions at high transverse momentum, which we have collected in table 2. As 
the comparison of (5.42), (5.52) and (5.56) shows, a common feature of all channels is the 
dependence on the transverse momenta k and r. 

Position space. The Feynman graphs for the splitting contributions are naturally eval- 
uated in momentum representation. We now transform our results to position space. We 
restrict our attention to the splitting g ^ qq since the other distributions can be treated 
in close analogy. Using the relation 
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(5.61) 



one can easily transform (5.42) to impact parameter space, 
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(5.62) 
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Fg,q FAg,Aq ^Ag.q FgAq Fsg,5q P&g,q ^Sg^Aq ^gM ^^g-^q 



ft 



Fg,g F^gAg F^g,g ^gAg Fig,Sg F&g,g ^SgAg Fg,Sg FAg,5g 

ff X X X X X X X 

hf^ XX X X X X X 

Table 2. Overview of the two-parton distributions that receive nonzero contributions from pertur- 
bative sphtting of a single quark or gluon. A cross indicates a nonzero contribution at order a^. 
Not shown are entries for Fq^g and its analogs with polarization, which are like those for Fg^q and 
its polarized counterparts. 



where 



zi — 2:2 



(5.63) 



fi{x,C) and hi^{x,(^) are the Fourier transforms of the transverse- momentum dependent 
parton densities in (5.42). The distribution in (5.62) is singular at y — ^zi = —\z2 and at 
y + i^zi = ^Z2, i.e. at the points where in the operator definition of F{xi, z.i, y) either the 
two fields on the left or the two fields on the right of the final-state cut are taken at equal 
transverse positions (see (2.85)). 

The form of the distributions -Fa^^a2(xj, fcj, y) in the mixed representation is slightly 
more involved. For convenience we introduce the function 



D{k,r) 



(fc+ir)(fc-^r) _ 1 
(k + W)'(k-W)' ~ 2 



4' 



1 {k'^ + \r^ + tkrY 



(5.64) 



which appears in Fq^q, F/\q^/^q and Fgq Sq- Its Fourier transform can be written as 



D{k,y) 



(27r)2 



D{k,r) 



where the Gaussian integral over r is straightforward to perform. Using the representation 



da oP ^ 



exp 



-aa 



z 

Aa 



1 



Kp{az) 



(5.66) 
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of the MacDonald functions we obtain 



D{k,y) = - 



dt e 



~2itky 



(2itky - l) Ko(2a/i - t"^ \k\\y\ 



+ 2\/l - t2 \k\\y\ Ki [2\/l - i2 
For the factor appearing in Fq^/^q and F/\q^q one finds in a similar fashion 



(5.67) 



D{k,y) 



d'r ,_^6'''(fe + ir)'(fc-ir)'' 



(2vr) 



{k+\r)\k-\rf 



W 1 I /' 2 
e k^y - 

TT 



dte^'^'^^y KQ(2^l-t'^\k\\y^ 



(5.68) 



For small k^y^ we can approximate the MacDonald functions and perform the integral 
over t, which gives 



D{k,y) 



1 



TT L 



log(4fc2y2) + 27 + 0(fc^y 



D{k, y) = e'''fc'/ - \og{Ak'y') + 27 - 2 + 0(fc^y 

vr L 



'2 r. 



2„.2n 



(5.69) 



where 7 is the Euler number. In the short-distance limit y"^ ^ ^we thus have a logarithmic 
divergence in D{k,y) and hence in the distributions Fq^q{xi, ki,y), F/\q^^q{xi,ki,y) and 

Fiq^Sq{Xi^ fcj, y). 

5.2.3 Contribution to the cross section 

We now investigate how the splitting contribution to quark-antiquark distributions in fig- 
ure 33a enters in the cross section for double hard scattering, as shown in figure 37a. 
Concentrating on the factors that depend on transverse momenta, we have with the cross 
section formula (2.33) and the distributions from (5.42) 

2 ^ " /" (P'T 

1^ j d ki d ki 5^ ^ [qi — ki — ki) j -^11:^2 (-^j' ^) ^ai,a2 (^j) ~ 

j d'^K.d'^R.6^'^\K + K - - [ffixi + X2, k) /f (Xl + X2,K,) + ■ ■ ■] 

d'rd'kd'k6('Hk + k-q) (^-^^)'(^ + ^^)' (5.70) 

{k-lr)\k + lrf [k + \r)\k-\rf 

with q = \{qi — Q2)) where the ellipsis stands for terms involving the gluon Boer-Mulders 
functions /i^^. Changing integration variables to fc+ = k -\- and k- = k — ^r, we can 
rewrite the last line as 



oc 



d'k^ 



d'k^ ^ 



(5.71) 



- 124 - 



Each integral is infrared finite but has a logarithmic divergence at large k± . This logarith- 
mic divergence also appears if we use the impact parameter space representation (5.62). 
According to (2.36) the cross section is then proportional to 



j d^Zid'z2d'y e "^^^^ "^^^^ Fa,-a^{Xi,Zi,y) Fa,,a2{Xi,Zi,y) 

oc j d'C e-'^'^^^"^^^ [/f (Xi + X2X) /f (Si + S2, C) + . . . ] 

X / d^z d^y e-- iy-h-yjy-yr {y + h-y'{y + hr' 

The last line diverges logarithmically for y = \z and y = —\z. At these points one 
respectively has y — \z\ = —\z2 and y + ^zi = \z2, so that the singularities correspond 
to configurations where partons are at the same transverse position, either to the right or 
to the left of the final-state cut. 

To understand the origin of this ultraviolet divergence, we go back to the graph in 
figure 37a. As mentioned in section 5.2.1 this graph receives leading contributions from 
two kinematic regions. In the first region, the virtualities and transverse momenta of the 
quarks are of order qx and thus much smaller than Q, whereas in the second region they 
are of order Q. The approximations that are necessary to derive factorization for double 
hard scattering are only valid in the first region. However, the integrand in (5.71) does 
not decrease fast enough with k± = \{k,i — k2 ^ r) to suppress the second region, so 
that the factorization formula (2.33) requires a suitable regularization in order to remove 
contributions from that region. A corresponding statement holds in the position space 
formulation. 

A simple way to regularize the cross section formula in impact parameter space is to 
impose a lower cutoff l//-i^ on (y + ^2)^ and [y — ^z)"^. The integral in the last line of 
(5.72) then becomes 



TT 



5'" / — Jo(u;) + U''" 



(5.73) 



times the same expression with I — )• /' and m — )• m' . The integral in (5.73) behaves like 
log(/i/|q|) for 1^1 S> \q\. The fi dependence of the cross section obtained in this way must 
cancel when one adds the contribution from figure 37a in the region of transverse loop 
momenta of order Q. That region is naturally associated with single hard scattering as 
discussed in section 5.2.1. At this point, one must obviously be careful to avoid double 
counting between the parts of the graph that one associates with single or with double 
hard scattering. The analogous double counting problem in multijet production has been 
pointed out in [141]. 

To use a cutoff in (5.72) is of course rather ad hoc, and there should be better ways 
to construct a consistent factorization scheme in which the formula for double-parton scat- 
tering has a controlled ultraviolet behavior and in which the double counting problem is 
properly taken care of. One may for instance think of subtracting the perturbative splitting 
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contribution of figure 33 at large momenta or small transverse distances in the definition 
of the two-parton distributions, so that graphs like in figure 37 are not included in double 
hard scattering at all. To solve this issue is a nontrivial task and must be left to future 
work. 

We already remarked that the integrals in (5.71) are finite in the infrared. This is due 
to the numerator factors and can be understood in simple physical terms, as noted in the 
detailed analysis given in [142]. The points where one of the four momenta fc+ — q, 
k- or fc_ — q vanishes correspond to configurations where one of the four g ^ qq splitting 
processes in figure 37a proceeds in strictly collinear kinematics. The amplitude for the 
collinear splitting g ^ qq is zero because an on-shell gluon has helicity ±1, whereas the 
helicities of q and q add up to zero due to chirality conservation for massless quarks. 

Referring to the end of section 5.2.1 we finally determine the dependence of (5.70) on 
qT and on A. With \qi + q2\ the second line scales like 1/A^, and with the behavior 
of the third line just discussed we find 



s^da 



fig. 37a It 



Hj^i dxi dxi d'^q. 
where /i is an ultraviolet cutoff much larger than q^- 
5.3 Parton splitting in collinear distributions 

The results of the previous section are relevant not only for transverse-momentum de- 
pendent two-parton distributions but also for collinear ones. As we have seen, collinear 
two-parton distributions appear in transverse-momentum integrated cross sections and in 
cross sections at perturbatively large qr via the ladder graphs discussed in section 5.1. 
Since ki and /c2 are not fixed in collinear distributions, the splitting contributions we com- 
puted in section 5.2.2 are relevant for F{xi,r) at large r and, after Fourier transform, for 
F{xi,y) at small y. 

5.3.1 Ultraviolet behavior 

Integrating (5.42) over ki and k2, i.e. over k and k, one formally obtains 

j2 



n / (2vr)2 _ Kr) \^^^_ - ^ /f (xi + X2) Ti[-, 

d^k {k+lrfik-^rY 
(2^)2 (fc+ i^)2(;,_ 1^)2 



Xi+ X2_ 

(5.75) 



where the integration over k gives the collinear gluon distribution ff^xi +X2), whereas the 
term with h^^ disappears due to rotation invariance. In the case where r^f^j 

oc 5^^', i.e. 

for Fg^g, F/\g^/^g aud Fsg^sq, the integral over k is ultraviolet divergent. The corresponding 
integrals of Fg^^g and -^A<j,g are proportional to e^^ rV and hence vanish, as they must 
according to the constraint (4.13) from parity invariance. An analogous discussion can 
be given for the interference distributions Iai,a2 and for the distributions resulting from 
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the splitting processes q ^ gq or g ^ gg. In all cases, contributions going with the 
Boer-Mulders functions h^'^ or h^^ vanish after integration over k and one is left with 
contributions from the unpolarized distributions f1 or ff. 

With the kernels T or V given in (5.43), (5.45), (5.53) and (5.57) we find that the 
splitting mechanism generates nonzero collinear two-parton distributions 

^q,q ) ^Ag,Ag ) PSq,5q ) ; lAq,Aq ; ISq,&q i -^9,9 ' -^AjjAg ; ^&g,q i 

Fg,g ,FAg,Ag ,Fg^Sg ,Fsg,Sg , (5.76) 

as well as the distributions obtained by interchanging the first and second subscripts in 
(5.76) or by replacing quarks with antiquarks in Fg^g and its polarized counterparts. With 
the exception of Fsg^Aq (and F^g^sg^ F^g /^q, F^g^Sg) these are indeed all collinear distribu- 
tions that are allowed by parity invariance and that are chiral even. For the distributions 
depending on polarization indices we have 

pjj' Tjj' xjf pjf,kk' jj',kk' 

^5q,5q ^Sq,5q " ' Sg,Sg ^ ^ ' 

^Sg,q Pg^Sg 2t^'J''''W' = 2r^ T^' - b"' . (5.77) 

To further investigate the ultraviolet divergence mentioned below (5.75) we focus on 
Fg^q for definiteness. Since D{k,r) in (5.64) falls off as l/k^ for fixed r and as for 
fixed fc, one obtains logarithmic divergences if one integrates over one or both of these 
variables. To regulate these divergences one may work in 4 — 2e dimensions. The result for 
Fg,g{xi, ki, r) is then the same as in (5.42) with a modified kernel 

r^;:(n; e) = j''' [u" + (1 - uf - e] . (5.78) 

Integrating over both transverse momenta and changing integration variables to = 
k + and k^ = k — ^r, one obtains 

which is zero due to rotation invariance. We note that integrating over ki, k2 and r puts 
all four fields in the matrix element defining Fg^g at the same transverse position, so that 
one obtains a twist-four operator. If (5.79) were not zero but finite after subtraction of 
the logarithmically divergent pieces, the graph in figure 33a would contribute to the scale 
evolution of a twist-four distribution. The vanishing of (5.79) thus reflects the fact that 
distributions of twist four and of twist two (the collinear gluon distribution in (5.75)) do 
not mix under evolution. The same zero result is obtained in any regularization scheme 
that respects rotational invariance. 

Integrating over k at fixed nonzero r and using the integral representation in (5.64), 
one obtains 

^2n ^\ /„2- 



/.-.m.)..'-.Iilf^r,.)(r,)- 



TT 



1 

- + \og^ + const. + 0{€) 



(5.80) 
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This contains an ultraviolet pole and an associated logarithm of the renormahzation scale fj?. 
The value of the constant is not of relevance for our discussion. If r = then the integral 
on the l.h.s. is scaleless and therefore vanishes. To isolate the ultraviolet singularity in that 
case, one can for instance give a small mass to the quarks. Up to corrections of order jtt,^ 
this leads to 

The ultraviolet pole and the associated logarithm are hence the same as for nonzero r. 

If one defines the collinear distribution Fq^q{xi,r) = J d?ki d?k2 Fq^q{xi, ki,r) in the 
MS scheme, the above 1/e pole is subtracted, together with a constant. To leading order 
in Os one finds for the scale dependence^^ 



d 



d log /i^ 

where 



FgA^^,r) = ffixi + X2) Pq,g ( , (5-82) 

g^qq Xl + X2 \Xl + X2 J 



is the familiar DGLAP splitting function (now including a color factor Tr = 1/2, unlike 
the function Pqg we used in section 5.1.3). We come back to this in the next section. Let us 
note that with the results in (5.52), (5.53) and (5.56), (5.57) we obtain relations analogous 
to (5.82) for Fg^q and Fg^g. On the r.h.s. of these relations we respectively find the DGLAP 
splitting functions Pg^q{u) and Pg^g{u), except for terms proportional to 5{1 — u) in Pg^g{u). 

Quite interestingly, the situation changes if we consider Fq^q{xi, y) instead of Fq^q[xi,r). 
The Fourier transform of (r^)^*^ in 2 — 2e transverse dimensions is 

I d^~^'r e-'^y (r^)-^ = A^-^^tt^-^ ^^T{^ , (5-84) 

which can be seen by writing (r^)"*^ = r~^(e) Jq°° da q*^"^ e""*"^, performing the integral 
over r and then the one over a. The factor r(e) responsible for the ultraviolet divergence 
in (5.80) is thus canceled if one Fourier transforms from r to y, and the result is finite for 
e = 0, 

^e"'''Dik,r) = ±. (5.85) 

The l/y^ behavior can be obtained directly in 4 dimensions by setting z = in (5.62). 

We thus find that Fq^q{xi,r) requires an ultraviolet subtraction for the graph in fig- 
ure 33a, whereas Fq^q{xi,y) does not. Let us see what we obtain if we define a modified y 



^^We note that both Qs and the gluon distribution in (5.75) also have a scale dependence, which 
becomes relevant at order in the evolution equation. 
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dependent distribution as the Fourier transform of the ultraviolet subtracted distribution 



[jimod 

q,q 



(2vr)2 



We have 



Fq,q{Xi-, y) 



f{xi,X2)TTlog^ + g{xi,X2) , 



g^qq 



g^qq 



1 



f{xi,X2) , 



(5.86) 

(5.87) 
(5.88) 



where the explicit expressions of / and g are easily obtained but not relevant for our 
discussion. As shown in appendix A the Fourier transform of (5.87) gives 



g^qq 



fixi,X2) 



1 



y' 



+ gixi,X2)6^'^\y), 



where 



" 1 " 






= lim 


y2 


+ (M) 



y^ 



0{y' 



S^'Hy)ldy^9{y''-e)9{bl-^,y') 



(5.89) 



(5.90) 



with bo = 2e~^ . We thus find the same y dependence in (5.88) and (5.89), up to terms 
concentrated at the singular point y = 0. This is not surprising since the ultraviolet 
divergent term in (5.80) is independent of r. 

In section 5.1.1 we have shown that the contribution of ladder graphs at large y to the 
cross section involves an integral 



/ 



(fyF{ui,y)F{ui,y) 



(5.91) 



according to (5.11) and (5.14). Since the collinear two-parton distributions behave like 
1/y^ at small y, the above integral has a linear divergence for small y^ and is hence not 
defined as it stands. A corresponding linear divergence is found if one Fourier transforms 
from y to r, 



I 



(fr F{ui, r)F{ui, —r) 



(5.92) 



where according to (5.87) the distributions behave like log(r^/;U^) for large r. The ultra- 
violet subtraction already included in the definition of F{ui,r) is hence not sufficient to 
render the integral in (5.92) finite. 

The reason for the unphysical divergences in (5.91) and (5.92) is that the cross section 
formulae containing these integrals have been derived for the region where 1/y'^ or is 
much smaller than q^. We thus encounter a similar problem as in section 5.2.3, with the 
difference that the divergence to be regulated is now linear instead of logarithmic. To 
make the cross section formulae (5.11) and (5.14) well-defined, one must either remove or 
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suppress the y integral in the region where \y\ is not large compared with l/qx, or one 
must define F{xi,y) such that in this region the contribution from perturbative sphtting 
as in figure 33 is subtracted. Along with such a procedure, one must provide a prescription 
for evaluating the splitting contribution at small |^/| in such a way that there is no double 
counting, as discussed in section 5.2.3. 

Integrating the cross section (2.91) over Qi and q2, we readily obtain the integral in 
(5.91) with Ui = Xi and Ui = Xi. The discussion of the previous paragraph carries over 
to that case, with the difference that the requirement for the validity of the cross section 
formula is then \y\ ^ 1/Q instead of \y\ ^ ^/qt- For the corresponding momentum 
integral (5.92) with Ui = Xi and Ui = xi one must require |r| ^ Q instead of |r| ^ qx- 
We note that in [143] it was proposed to regulate this integral by imposing an upper cutoff 
< min(g^,g|). By itself this is clearly insufficient to obtain a reliable result, since the 
contribution from outside that region is large and needs to be evaluated as well. 

5.3.2 Scale evolution 

Let us now investigate the scale evolution of collinear two-parton distributions. We focus 
on the color-singlet combinations ^F, which are most closely related with single-parton 
densities as we already saw in section 3.5. For definiteness we consider the quark-antiquark 
distribution ^Fq^q, which we studied extensively in the previous section. The generalization 
to other parton and polarization combinations is straightforward. 

The dependence on the scale ^ of collinear parton distributions arises from the reg- 
ularization and subtraction of ultraviolet divergences in their definition. This involves 
divergences from self-energy graphs (which also occur in transverse-momentum dependent 
distributions and can be expressed in terms of suitable Z factors) and divergences from re- 
gions of large transverse parton momenta. For a single-parton distribution the contribution 
from the high-transverse- momentum tail was already discussed in section 5.1.2. 

The scale dependence in the collinear distributions ^F(xi,y) arises from self-energy 
graphs and in addition from the ladder graphs in figure 30, which according to (5.12) 
and (5.13) give rise to ultraviolet divergent integrals f d?kid'^k2 ^F{xi, ki,y) unless one 
performs suitable subtractions. Since the ladder and self-energy graphs have exactly the 
same structure as for single-parton distributions, the corresponding evolution equation 
reads 

+ [ ' —Pg,bJ—)%,b2{xuU2,y) (5.93) 

for a quark-antiquark distribution. The splitting functions Pa,fe now include the contribu- 
tions from virtual corrections, unlike the corresponding kernels in section 5.1.3. Note that 
the labels bi and 62 do not take the value 5g here, because the corresponding kernels Pq^Sg 
and Pq^Sg vanish due to rotation invariance, see our remark below (5.18). 
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The evolution equation (5.93) has the structure of a usual DGLAP equation for each 
parton. The corresponding operator appearing in the distribution Fq q{xi,y) is 



{g(-iz2)M^[-iz2,iz2]r,-g(iz2)} 
{q{y-hi)W[y-lzi,y + lzi]rgq{y + ^zi) \ ^'^ 



(5.94) 



where the Wilson line VF[i^',^] is defined in (3.113) and where { . . . Y^^'^^ indicates that 
each bilinear operator qWq is renormalized at scale fi in the same way as for single-parton 
distributions. As long as the transverse distance y between the two bilinear operators 
remains finite, no further ultraviolet divergences appear, and one has the product of two 
renormalized twist-two operators. As remarked earlier in the literature, one may choose 
different renormalization scales //i and for the two operators, which appears useful 
when one has two hard-scattering processes with rather different hard scales. The separate 
evolution equations in fii and fi2 are then simply the usual ones with a single DGLAP 
kernel. 

For the collinear distributions ^Faj^aj (^i) that depend on the relative momentum r 
the situation is different, as we have seen in the previous section. The splitting graph in 
figure 33a and higher-order corrections as in figure 33b give rise to additional ultraviolet 
divergences. Their subtraction leads to an inhomogeneous term in the evolution equation. 
At leading order in Os one has 

Fg^g{xi,X2,r) = > / Pq,bA—] Fb^,gi^i,X2,r) 



d log ji 



bl=q,g •''^l 
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1-^1 dm 



+ yZ I ~ ^9,62 ( — ) %,b2 {xi, U2, r) 

^ 1x2 U2 \U2j 



"1, 

1, 



P. 



X1+X2 ''\Xi + X2 



/f(xi + X2) 



(5.95) 



where the extra term follows from (5.82). At higher orders in as, the inhomogeneous term 
will involve a convolution integral, as can be anticipated from the graph in figure 33b. 
The appearance of the extra term in the evolution equation can also be understood in the 
impact parameter representation by writing Fq^q{ Fourier transform 



d^ye '''y Fq^q{xi,y;n) 



- rcn, fi 



(5.96) 



Since Fq^q{xi,y) has a 1/y'^ singularity at small y, the integral over this variable is log- 
arithmically divergent and requires a subtraction in addition to those already made in 
Fq,q{xi,y). We have indicated this extra subtraction by [ . . . Y'^'^'^. 



For the distribution 



F, 



g,q\Xi, fl) — Fq q{Xi,r — 0, /i) 
dcf 



d^y Fq^q{xi,y;fi) 



(5.97) 
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the evolution equation (5.95) has long been known in the literature, see [144-146] and the 
recent detailed study [147]. We wish to comment in this context on an ansatz that is often 
made in phenomenological studies, in which the y dependent two-parton distributions are 
written as 

F{xi,y;fi) = f{y)Fixf,fi), (5.98) 

where f{y) is a smooth function normalized as J d?y f{y) = 1. A typical choice for f{y) 
is e.g. a Gaussian or a sum of Gaussians. This type of ansatz is obviously inconsistent if 
F{xi,y; fj,) is defined from the product (5.94) of twist-two operators, since the /i dependence 
of the l.h.s. is then given by the homogeneous evolution equation whereas the /x dependence 
on the r.h.s. is governed by the inhomogeneous evolution equation (5.95). If one instead 
defines the y dependent distribution as the Fourier transform of F{xi, r) as in (5.86) then 
the ansatz (5.98) is consistent regarding evolution since by construction F'^°'^{xi,y; fi) 
evolves as in (5.95). We do however not think that this procedure is satisfactory. As 
we have seen in (5.89), F{xi,y)'^°'^ differs from F{xi,y) only by terms proportional to 
S^'^\y), and such terms do not appear in the ansatz (5.98), which is smooth and finite at 
y = 0. In more physical terms, we recall that the inhomogeneous term in the evolution 
equation (5.95) has its origin in the 1/y'^ behavior of F{xi,y) at short distances, which is 
not described by (5.98). 

We have seen in section 5.2.3 that this short-distance behavior prevents us from using 
either F[xi,y) or F'^°'^{xi,y) in the double-scattering factorization formula as it stands. 
An ansatz like (5.98) with a smooth function f{y) does not have this problem and may be 
regarded as modeling a y distribution where the perturbative splitting contribution that 
gives rise to the l/y^ singularity has been removed. Since the ansatz is ad hoc, one cannot 
say which evolution equation should then be used on both sides of (5.98). Our discussion 
suggests that the homogeneous form (5.93) may be more appropriate, at least for values 
y of typical hadronic size, which are of course most important when the ansatz is used in 
the factorization formula. With this choice, one also retains consistency with respect to 
evolution if one makes the additional ansatz F{xi,fi) = /(xi, ^u) /(x2, /u), as is often done. 
To find a systematic solution that treats both splitting and non-splitting contributions in 
a consistent manner remains a task for future work. 

For the reasons discussed in section 3.5, the evolution of color octet distributions ^F 
differs from the one of ^F, and these differences have not yet been worked out in detail. 
However, the issues discussed in the present section affect ^F in the same way as ^F, given 
that both the ladder graphs in figure 30 and the splitting graph in figure 33a differ only by 
overall factors between the singlet and octet channels. They hence give rise to the same 
logarithmic divergences when the relevant transverse momenta are integrated over. One 
may therefore expect that, once a solution of the above problems for singlet distributions 
is found, it will be possible to adapt it to the octet sector. 
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6 Conclusions 



We have investigated several aspects of multiparton interactions in QCD. Such interactions 
can contribute to hadron-hadron cohisions whenever one has a final state with several 
groups of particles for which the vector sum of transverse momenta is small compared with 
the large scale Q that characterizes the process. As we have shown in sections 2.1.3 and 2.4, 
multiple interactions are then not power suppressed in 1/Q compared with the mechanism 
where these groups of particles are produced in a single hard scattering. Examples are 
the production of two lepton pairs originating from the decay of two vector bosons with 
low transverse momenta, or the production of two dijet pairs that are approximately back- 
to-back. For small parton momentum fractions x, which are typical of collisions at the 
LHC, multiple hard scattering can even be enhanced because one expects multiparton 
distributions to rise faster with decreasing x than single-parton densities, as we argued in 
section 2.5. 

Given the importance of transverse momenta in the final state, we have given a fac- 
torization formula for multiple hard scattering in terms of multiparton distributions that 
depend on the transverse momenta of the partons. Such a formula can be fully derived 
for lowest-order Feynman graphs and generalizes the more familiar description in terms 
of collinear (i.e. transverse-momentum integrated) multiparton distributions given in the 
literature [52, 53, 66]. A physically intuitive interpretation is obtained if one expresses the 
cross section in a mixed representation, in which the multiparton distributions depend on 
the average transverse momentum of the partons and on their average transverse distance 
from each other, where the "average" refers to the scattering amplitude of the process and 
its complex conjugate. These distributions have the structure of Wigner functions. 

The simple picture just sketched is however complicated by the presence of correlation 
and interference effects, some of which have been pointed out earlier in the literature 
[89]. The spin and color of the partons described by a multiparton distribution can be 
correlated, and such correlations change the overall rate of multiple interactions. Two- 
quark distributions allow two color couplings, which we classified as color singlet and color 
octet, whereas for gluons a number of color couplings appear in addition to the color singlet 
one, see section 2.3. In section 4.1.2 we have shown that spin correlations can also affect 
the distribution of particles in the final state, using four-lepton production as an example. 
Further contributions to the cross section can come from interference effects in fermion 
number or in quark flavor (figures 6c and 7) and from the interference between single and 
multiple hard scattering (figure 9a). One can however expect that these interference effects 
will not benefit from the small-x enhancement of multiple interactions mentioned above 
(although in the case of interference between single and multiple scattering the situation 
is not entirely settled as explained in section 2.5). Regarding "rescattering contributions" 
of the type shown in figure 12a, we have shown that their evaluation in terms of two 
sequential scattering processes with on-shell external partons is inappropriate and that, 
when calculated properly, such contributions are suppressed by powers of 1/Q. 

How large the above correlations and interference effects are remains an important 
open question, both for the phenomenology of multiple interactions and from the point of 
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view of hadron structure. The possibility to study moments of multiparton distributions on 
the lattice as explained in section 4.2, as well as the approximate relations with generalized 
parton distributions we derived in sections 2.1.5 and 4.3 provide two possible avenues to 
investigate these issues further. 

A proper factorization formula in QCD requires much more than an analysis of the 
lowest-order Feynman graphs contributing to the process in question. In section 3 we 
have taken first steps towards a factorization proof for double hard scattering in terms of 
transverse-momentum dependent distributions. Our investigation only applies to processes 
where each hard scatter produces color-singlet particles, given the limitation of our current 
understanding for single hard-scattering processes [111]. For definiteness we have restricted 
our analysis to the double Drell-Yan process. We have shown how collinear and soft gluon 
exchange at order can be arranged into Wilson lines, which are basic building blocks 
in the construction of an all-order factorization formula. We have also seen that at this 
order soft-gluon effects cancel in factorization formulae that involve collinear two-parton 
distributions in the color singlet sector, whereas they do not cancel in the color octet sector. 
In section 3.2 we have listed the many issues that remain to be clarified and worked out for 
a full factorization proof. The most critical questions are probably whether one can show 
that the effect of soft gluons in the Glauber region cancels in the cross section and whether 
the double counting problem mentioned below can be solved in a satisfactory way. 

Our calculation of soft-gluon effects at leading order in Os also allows us to investigate 
the structure of Sudakov logarithms in the double Drell-Yan process, extending the method 
of Collins, Soper and Sterman [118]. We find that the leading double logarithms are given 
by the product of the corresponding Sudakov factors for each single scattering process, 
whereas beyond this approximation soft gluon effects connect the two hard scatters in 
a nontrivial way. In the region where all transverse parton momenta are large and the 
transverse distance y between the two partons is small compared to a hadronic scale, we 
find that Sudakov effects favor the color singlet coupling in two-quark distributions. If this 
result could be generalized to large y, it would provide a valuable simplification. 

In generic kinematics, the description of multiple interactions involves a multitude of 
terms, with many unknown distributions that describe correlation effects already in the 
case of double hard scattering (not to speak of the case with three or more scatters). 
The predictive power of the theory is increased in the region where the net transverse 
momentum qt for each final state produced by a hard scattering is large compared with 
the scale A of nonperturbative interactions (while still being small compared with the scale 
Q characterizing the hard-scattering processes). Apart from the possible simplification 
due to Sudakov effects just mentioned, the transverse-momentum dependent multiparton 
distributions can then be computed in terms of collinear distributions and a hard scattering 
at scale qt- The generation of high transverse momenta can proceed by ladder graphs as in 
figure 30, and we find that the color factors of these graphs favor the color singlet coupling 
in two-parton distributions. 

A different mechanism is shown in figure 33, where one parton splits into two partons 
that subsequently take part in a hard scatter. By explicit calculation at order Ug we find 
that this splitting mechanism generates a multitude of spin correlations between the two 
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emerging partons. For qq distributions the color singlet coupling is preferred, whereas for 
qg and gg distributions the opposite is the case. Contributions from ladder graphs and 
from parton splitting graphs compete with each other in the double scattering cross sec- 
tion. An overview is given in table 1, where we see that compared with splitting graphs 
the contribution of ladder graphs is suppressed by powers of A/gj-. On the other hand, the 
splitting graphs lack the small-x enhancement discussed earlier, so that one cannot decide 
on generic grounds which mechanism is more important in given kinematics. Finally, we 
find that splitting contributions require a modification of the formalism outlined so far, 
because they increase so strongly for decreasing interparton distance y that one obtains 
divergent integrals when inserting them into the factorization formulae. This is closely 
related with the problem that graphs like in figure 37a can either be interpreted as repre- 
senting double hard scattering with parton splitting in each two-parton distribution, or as 
representing a single hard-scattering process at two-loop level. A consistent factorization 
scheme must ensure that there is no double counting of this graph in different kinematic 
regions. A satisfactory solution of these problems remains to be found, and as we argued 
in section 5.3.2 such a solution will also have consequences on the evolution equation for 
collinear multiparton distributions. 

In summary, we find that a systematic description of multiparton interactions in QCD 
involves a considerable degree of complexity, but that there are several elements that hint 
at possible simplifications. More work is required to work out these simplifications and to 
put the theory on firmer ground. 



A Two-dimensional Fourier transform of the logarithm 

In this appendix we prove the relation 
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with bo = 2e , which we used in (5.89). To this end we integrate the relation over a test 
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which leads to 
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The integration region has been split in such a way that integration by parts does not give 
any boundary term, so that one has 



For an alternative derivation of the result in this form (with test functions depending on y^) 
we refer to eqs. (129), (133) and (141) in [148]. Further rewriting 



we can separate the terms with f{y) and with /(O) and thus obtain (A.l). 
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Changes in arXiv version 2 compared with version 1 
abstract: removed reference to lattice calculations 
p.l added Ref. [39] 
p. 2 added Ref. [40] 

p. 5 2 lines above (2.11): specified transverse positions and momenta 

p. 8 in (2.21) and everywhere else in the paper: changed notation for the symmetry factor 
from S to C, so as to avoid confusion with the cm. energy. Added a footnote 
comparing with a common notation in the literature. 

p. 15 second line of (2.54) corrected order of arguments p and p' in F 

p. 37 second line: added a forward reference to a later section 

p. 39 6 lines below (2.137): corrected "graph 9" to "graph 9d" 

p. 50 figure 15: corrected assignment of open and filled circles 

p. 89 omitted unnecessary (and confusing) subscripts on a and a in (4.19) 

p. 91 extended title of section 4.2 

p. 94 2 lines above (4.35): changed formulation to avoid confusion between A as a baryon 
and as a momentum transfer 

p. 102 equation (5.10) and text two lines below: corrected ^2 ~ to q2 — ^2 

p. 108 3 lines below (5.25): updated comparision of our results with [138] and added 
comparision with (new) Ref. [140] 

p. 117 9 lines before section 5.2.2: updated comparision with Ref. [98], which has changed 
from arXiv version 1 to 2 

p. 119 9 lines above (5.45): added a small paragraph to help orient the reader 

p. 123 corrected exponent on the l.h.s. of (5.66) 

p. 126 corrected power of logarithm in (5.74) 

p. 130 5 lines above section 5.3.2: corrected \y\ ^ qt to \y\ ^ 1/q't 

In addition, a small number of grammatical mistakes have been corrected. 
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